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Abstract 

We physicists are on a journey towards the ultimate theory which describes everything in 
our Universe. A great milestone achieved in 201 2 is the discovery of the Higgs boson by the 
CMS and the ATLAS collaborations at the Large Hadron Collider (LHC); it completes the 
Standard Model of particle physics, which was developed through the mid to late twenties 
century. 

The long twentieth century was over. Happily making a steppingstone of the Standard 
Model, we are now heading to more fundamental theories. Nature has many unsolved 
features: the Dark Matter, the Dark Energy, and the mechanism which produced current 
baryon asymmetry of our Universe, etc. Also we have to build a unified explanation of 
the three forces embedded in the Standard Model, and to develop a description of the 
gravitational force in the language of the quantum theories. 

The supersymmetry is one of the most promising candidates for theories beyond the Stan- 
dard Model, and its tail was expected to be caught at the early stage of the LHC. However, 
our expectation was not fulfilled, and we have no footprints observed yet. What does this 
mean? 



The current status of the supersymmetric theories is described in this dissertation. First, 
the simplest model of supersymmetric Standard Model is introduced, which is called the 
MSSM. Under this framework, the discrepancy in the muon anomalous magnetic moment 
between the prediction from the Standard Model and the experimental result suggests the 
supersymmetric particles are of order "lOOGeV, which is also supported by discussions on 
the little hierarchy problem. However, the LHC experiments have found no scalar-quarks 
or gluinos in such mass range, and moreover, the Higgs boson mass of 126GeV requires, 
within the MSSM framework, the scalar-top mass of order 1-10TeV. This current status 
forces us to abandon the simplest supersymmetry-breaking frameworks of the CMSSM 
and the GMSB scenarios. 

Two promising possibilities remain there: the first is that the scalar-quarks and the gluino 
are much heavier than of order "lOOGeV while the other SUSY particles remain near the 
order, and the second is to extend the MSSM with extra fields. The second scenario is 
investigated in this dissertation; the V-MSSM is proposed as an extension of the MSSM 
with a (10 + 10) pair of the SU(5) decuplets. In the framework the Higgs mass is increased 
by effect from the extra matters, and thus the 1 26 GeV is achieved with the scalar-top having 
a lighter mass. This fact resurrects the CMSSM and the GMSB scenarios. This dissertation 
examines the GMSB scenario under the V-MSSM; it is called V-GMSB scenario. 

It is shown that the V-GMSB has a potential to realize the 1 26 GeV mass of the Higgs boson 
with holding the explanation of the muon magnetic moment discrepancy, if the masses of 
the extra quarks are approximately less than 1.2TeV. Constraints on the V-GMSB from 
the LHC experiments are discussed then; it is concluded that the gluino mass must be 
approximately heavier than 1 .1 TeV, and that the extra quarks be heavier than 300-650 GeV 
depending on the decay branches of them. 

LHC prospects are briefly discussed. As the extra quarks are expected to be approximately 
less than 1.2 TeV, searches for the particles are of great interest at the 14 TeV LHC; con- 
straints from the supersymmetry search, especially on the gluino mass, are expected to be 
much improved there. Therefore, it is expected that the fate of the V-GMSB is adjudicated 
at the court of the 14 TeV LHC. 
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Chapter 1 



Overture 



A new boson is observed in Higgs hunting. The magisterial thesis of Author, submitted in January 2010, 
begins with the following paragraph: 

We have the Standard Model, which describes almost all physics below the energy scale 
100 GeV. Although it is still under verification, especially the existence of the Higgs boson, 
the experiments held in the Large Hadron Collider (LHC) will work out the answer soon, 
which will be a declaration of the triumph of our philosophy. 

This prognostication became reality. On 4th July 2012, the ATLAS and the CMS collaborations claimed that 
they respectively observed a new boson with a mass approximately 126 GeV in searches for the Standard 
Model Higgs boson Q]|2]. It is not confirmed that the new particle is the Higgs boson; we have to measure 
its property precisely, especially the Yukawa couplings with the Standard Model particles. Nonetheless, it 
allows us to expect that the Standard Model with the Higgs mechanism GJH) will be validated with the LHC 
and the International Linear Collider (ILC) in the near future. 

In this dissertation, we assume that the new particle is the Higgs boson, and discuss the Standard Model 
with the Higgs mechanism, simply we call it the Standard Model hereafter, and higher-energy theories above 
it. 

Standard Model has been completed, but. . . The discovery of the Higgs boson has completed the Stan- 
dard Model. This model has the electroweak symmetry breaking as its heart, which is governed by the Higgs 
mechanism, and explains almost all of Nature. 

We physicists, however, expect that more fundamental theories are hidden at higher energy scale Mh; g h, 
and hope that they solve the following prospects or problems: 

(a) unification of the electroweak force and the strong force ("grand unification"), 

(b) description of the gravitational force in harmony with the electroweak and strong forces, 

(c) the reason why neutrinos have such extremely tiny masses, 

(d) the mechanism which generated current baryon asymmetry of our Universe, 

(e) the mechanism which caused the inflation in the early universe, 

(f) identification of the Dark Matter and the Dark Energy, 

but these are so fundamental that cannot be solved in a night. Towards these problems physicists have 
examined the Standard Model, and have found the following problems, discrepancy, or unnaturalness, which 
can be used as keys towards more fundamental theories: 

((problem)) Dark Matter problem We know that our familiar matters, e.g., electrons, protons, and 
neutrons, account for approximately 4% of the substance of our Universe ||5]|6). We consider that 
approximately 20% is some other matter, called "Dark Matter," and that the rest is not even matter, 
which we call "Dark Energy." We do not know nature of them; as for the Dark Matter, since it is still 
considered to be "matter," expected to be ascertained more easily. 
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((problem)) current baryon asymmetry of our Universe Although the baryon number B is slightly 
violated with the sphaleron effect 0SI9), the violation is too feeble to explain current baryon asym- 
metry of our Universe. 

((unnaturalness)) fine-tuning problem (hierarchy problem) Based on our current understanding of 
quantum field theory, the mass of the Higgs boson, nth, should naturally be much above (9(100) GeV, 
but we found it is just around 126 GeV. This means that Nature employs unbelievably unnatural 



tuning. This topic is examined in Sec. 2.1.2 



((discrepancy)) muon g - 2 problem The anomalous magnetic moment of muons, or the muon g-2, 
is precisely measured in experiments, but the measured value has 3<x-level discrepancy compared with 



the theoretical prediction under the Standard Model, as will be discussed in Sec. 2.2 



((discrepancy)) slight mismatch on gauge coupling unification Physicists expect the three forces, 
the electromagnetic force, the weak force, and the strong force, would be unified in more fundamental 
theories. If so, the strength of these forces, expressed in terms of gauge coupling constants, should be 
common at the energy scale of such theories. Actually the coupling constants are dependent on energy 
scale, and the matching does roughly occur at ~ 10 16 GeV. However, there lies a slight mismatch, and 
we have to modify the energy dependence slightly in order to realize complete unification. We will 



discuss this topic in Sec. 2.3 



The SUSY and the MSSM 

Amazingly, we have a silver bullet. With examining those hints, physicists invented the supersymmetry 
(SUSY) IfTUl several decades ago. 

The SUSY is a symmetry between bosons and fermions: it transforms bosons into fermions, and vice 
versa. The SUSY extends the Standard Model. The minimal version of the supersymmetric Standard Model, 
the minimal supersymmetric standard model (MSSM) |[TTl[T2l[T3l has, therefore, scalar quarks (squarks) and 
scalar leptons (sleptons) as the bosonic partners of the Standard Model fermions. Also the Higgs boson and 
the gauge bosons meet their fermionic partners called Higgsino and gauginos. Chapter|4]of this dissertation 
is devoted to topics around the MSSM. 

The MSSM completely solves the fine-tuning problem. It also has capability to provide a candidate for 
the Dark Matter, and to solve the muon g-2 problem. Also the slight mismatch of the gauge couplings 
is resolved. Moreover, the SUSY is considered as a key to build string theories, which are considered as 
promising candidates for the ultimate theory. 1 

However, there ain't no such thing as a free lunch. The great MSSM has still many problems. The 
principal one is the fact that the SUSY is not realized in Nature. We do not have scalar electrons with 
a mass of 0.5 11 Me V. We thus consider the SUSY is violated so that the masses of the superpartners 
become much heavier. This is realized with appending extra terms that do not respect supersymmetry to the 
(supersymmetric) Lagrangian of the MSSM. Then this patch causes CP- and flavor problems, and another 
patch, called the /^-parity, is required to make protons stable. Another weak point of this model is that no 
evidence has been found at experiments, although the LHC is expected to observe such signals. It is nothing 
but a hypothesis, or a daydream, at this stage. These problems, and solution candidates, are discussed in 
Chapter g] 

*'A too philosophical (and not scientific) note: 

It is instructive that a symmetry between bosons and fermions solves the hierarchy problem. The hierarchy problem ultimately 
originates the fact that a scalar boson, the Higgs field, is appended to the Standard Model, where all matters are fermionic and all forces 
are governed by vector bosons, in order to realize the electroweak symmetry breaking. It is somewhat expedient that the electroweak 
symmetry breaking, the heart of the Standard Model, is realized by a strange, and lastly appended, particle, and we had to worry why 
we have a sole scalar particle, a muggle, in the theory of fermions. 

The history of physics is characterized as the cycle of unifications, from that of motions of apples and the moon to that of the 
SU(2) wi;il | t and the U(l)y gauge symmetry (in other words, that of the W ± bosons, the Z boson, and the photon). Therefore, the 
principal problem of the Standard Model, in this historical and philosophical viewpoint, is the co-existence of bosonic and fermionic 
matters. The Standard Model should be, "historically," supersymmetric, and it is one important reason that I, Author of this dissertation, 
prefer the SUSY as the model beyond the Standard Model. 
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About this dissertation 

We in this dissertation will see how the status of this SUSY daydream is altered by the discovery of the 
Higgs boson. 

After discussing the problems related to the Standard Model in the next chapter, we will review the 
MSSM and our current understanding on it in Chapter|4] The Higgs mechanism and the mass of the Higgs 
boson under the MSSM, and the mechanism to solve the muon g-2 problem, are examined there, and current 
status of LHC SUSY searches is reviewed. We will then notice that the mass of the Higgs boson, 126 GeV, 
is a bit heavier than natural expectations, and that the Higgs boson now requires the squark to have masses 
of (9(1-10) TeV. This fact forces us to abandon the SUSY solution to the muon g-2 problem as long as we 
adopt the most simple and neat set-up, called the gauge-mediated SUSY-breaking (GMSB) framework. 

To revitalize the GMSB framework, in Chapter [5] we extend the MSSM with extra particles. This 
model, called the V-MSSM, is the main topic of this dissertation. In the model the extra quarks yield extra 
contributions to the Higgs mass, and the 126 GeV is realized with the squarks having a mass of < 1 TeV. 
Such lighter masses of the SUSY particles allow the SUSY contribution to the muon g-2 to be large enough 
to explain the discrepancy even under the GMSB scenario. We call this the V-GMSB scenario.* 2 

Chapter [5] is the main chapter of this dissertation. There phenomenology of the V-MSSM is discussed. 
We start from the discussion on the Higgs boson mass, on the gauge coupling unification, and on the muon 
g — 2, under the V-MSSM. Then the V-GMSB scenario is introduced in Sec. |5.3| The characteristics of the 
model are examined, and the vacuum stability bound enters the discussion as a very severe constraint on 



the V-GMSB scenario. However, even under the bound, as can be seen quantitatively in Sec. 5.5 the Higgs 
mass of 126 GeV can be explained with holding the explanation of the muon g-2 discrepancy. Here the 
GMSB is revitalized. Finally, interpreting reports from the LHC experiments, we will obtain current collider 



constraints on the model in Sec. 5.6 and Sec. 5.7 Especially, Fig. 5.7 is the conclusive figure of our tour. 



Chapter[6]is the coda, where the main theme of this dissertation is recapitulated. 



2 The same discussion can be applied to the CMSSM framework; that is, in the CMSSM framework the muon g-2 anomaly cannot 
be solved under the 126 GeV Higgs constraint 1141 . and extending to the V-MSSM resolves this conflict 1 1 511 1 61 - 
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Chapter 2 

Foundation 



This chapter is a review of the problems stated in Chapter [T] of the Standard Model. We will in Sec. |2.1 
start from the Standard Model and review the electroweak symmetry breaking. The Standard Model is 
now completed with the discovery of the Higgs boson, but we will see it has a strange problem, "hierarchy 
problem," once we regard the Standard Model as a low-energy effective theory of more fundamental theories. 



This is the main topic of this chapter. Then in Sec. 2.2 we review the muon g - 2 problem, a discrepancy 



in the muon anomalous magnetic moment between the Standard Model prediction and the measured value. 
Sec. |2.3| is a discussion towards the grand unification theories (GUTs), where we will see a slight mismatch in 
unification of the gauge coupling constants. These problems are all solved with the supersymmetry (SUSY) 
in Chapter|4] 

Section 2.1 The Standard Model and Hierarchy Problem 
2.1.1 The Standard Model 

The Standard Model [3| is a model of the gauge theory. It has a gauge symmetry of Gsm = U(l)y X 
SU(2) wea k x SU(3) st rong- The symmetries appear in Nature as "forces," which in language of the quantum 
field theory are governed by gauge bosons: a B-boson for U(l)y, three W -bosons for SU(2) wea k, eight gluons 
for SU(3) strong . 

However, Nature does not have U(l)y x SU(2) wea k; this symmetry, called the electroweak symmetry, is 
spontaneously broken with the Higgs mechanism [4|, discussed below, and falls into a U(l) electromagnetic 
symmetry, which governs the electromagnetic force with the photon y. The remnant yields the weak force 
with the W ± -bosons and the Z-boson. This is the kernel of the Standard Model, the electroweak symmetry 
breaking. 

Let us start from the Standard Model Lagrangian of our triumph. 

-£ — -£gauge -^Higgs "t" -^-matter -^Yukawaj (2.1) 

Aauge = - \ B ^ B ^ ~ \ Tr(W v WV) - \ TtiCTG^), (2.2) 

Xffiggs = \{d M - mW» - ImB^ Hf - V(H), (2.3) 
Xmatter = Q.l*f - mG^ - xgiW^ - gitfyfip) PlQi 

+ Uiiy» (d M - - fitfrflj PrUi + Dpf (d M - igsG^ + \ig Y B^ P K D t (2.4) 
+ Ljiy (d M - igzWp + \ig r B,) P h L t + Epf (d„ + ig Y B M ) P R £/, 

^Yukawa = U i{y u ) u HP h Qj - Diiy^jH 1 P L Qj - Efo, ),,// Lj + H.c; (2.5) 

V(H)= ^A\H\ 4 -fx 2 \H\ 2 . (2.6) 
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Bfi, W/j and are the gauge fields, and their "field strengths" d^Xy - d v X^ - xgxiX^Xy] are denoted as X^. 
H is the Standard Model Higgs field, and the Higgs potential V(H) is characterized by the quartic coupling 
A and the quadratic coupling fi, where /i is the only parameter with mass dimension in the Standard Model. 
y„4 >e are matrices of the Yukawa couplings in the Standard Model* 1 ; i and j are generation indices of quarks 
(Q, U and D) and leptons (L and E). P L and P R are the well-known projection operators. This Lagrangian 
describes Nature very well, and exceptions are limited to be, e.g., the tiny neutrino masses and the Dark 
Matter. 

Here one should notice that the baryon number B and the lepton number L are respectively 
conserved in this Lagrangian; this is an accidental symmetry originating in the Standard Model 
gauge symmetry. Actually these symmetries are anomalous, and quantum effect causes especially 
in the early universe the so-called "sphaleron process," which violates B + L significantly (7][8][9). 
Meanwhile, the conservation of the number B - L is not anomalous, and thus it is kept even under 
quantum effect. 

Our main concern is the Higgs potential. For a successful electroweak symmetry breaking A > and 
/j 2 > must hold. Then the minimum of the Higgs potential falls in \H\ = /j/ ^J~A =: v, which results in the 
vacuum expectation value v of the Higgs boson; the value is well-known to be u « 174 GeV from the masses 
of the W- and Z-bosons. Finally, the Higgs field is parameterized as 

H( x )-( H+(x) )-(°\ + 1 (Mx)+i<t>2(x)\ (21 s 

H(x) ~ \h\ x )) ~ [vj + \ h(x) + ifc&c) j " (Z7) 

h(x) is the (Standard Model) Higgs boson, which we first observed in 2012. <pi(x) are would-be Nambu- 
Goldstone bosons, ignored here for simplicity. The potential is now 

VQi) = ^h 4 + + v 2 h 2 + constant. (2.8) 

What we observed in this year 2012 is nothing but = V2yU ss 126 GeV, and this results in A as 0.26. Now 
we know all the parameters in the Standard Model. The Standard Model is completed. 

The electroweak symmetry breaking makes the three W-bosons massive, which are observed as W ± - 
bosons and a Z-boson. This effect emerges from -Cmggs as , with ignoring the Nambu-Goldstone bosons, 

X H iggs = \{d, - ig 2 W M - Ifc r 2?J Hf d l -{d,hf + ^ (g 2 2 W +tl W; + ^Z%J . (2.9) 



Here redefinition of the gauge field is employed as 

V ± ■= — (W l + iW 2 ) , - , 

" V2 " " Wl' \ sin0 w cos0 w /\B /i 



Wt := ^(W^iW 2 ), l Z A := fe S f» -J^U?!. (2.10) 



where 



gz 



■ =^gr z + g2, gy = — =gzsm6 w , g 2 - . - -g z cosg w . (2.11) 

' inns H... sin c7 w 



y/gy 2 + gi 2 v cos0 ' 

We can find the mass of the gauge bosons as 



gi gz 

m-w — v, niz — v, (2.12) 

V2 V2 



where we can obtain the value of v ~ 174 GeV. 



'We will later use Y ul i e for the Yukawa couplings in the MSSM. 
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h 



k s 

(A) (B) 



Figure 2.1: The one-loop level diagrams which give quadratic quantum corrections to the mass of the 
Higgs boson, provided by (A) a Dirac fermion if/ and (B) a scalar boson (p. The particles are assumed to be 
coupled to the Higgs boson h with interactions kluf/iff and fc s |/i| 2 |0| 2 , respectively, where k and k s are coupling 
constants. 



2.1.2 Hierarchy problem 



The Standard Model with the Higgs mechanism is a complete beautiful theory for its own sake. However, 
we physicists do not consider the Standard Model as the ultimate theory, but a low-energy effective theory 
of more fundamental theories hidden at higher energy scales, which we label Mhi g h- Once we take this 
standpoint, which is the current paradigm of our science, the Higgs particle brings an unnaturalness into the 
Standard Model, which has been discussed for several decades, the hierarchy problem lTP7l . Let us briefly 
review the unnaturalness through the electroweak symmetry breaking in the Standard Model. 

The mass of the Standard Model Higgs boson receives quadratic quantum corrections from the diagrams 



in Fig. 2. 1 The contributions are respectively evaluated as, at one-loop level, 



Aw ^"^ Mh2 ^ +0(1 ° gMhi s h) 



from a fermion (Fig. 2.1 A), and 



Am h 



\6n 2 



T M h ia h +< 3(logM high ) 



(2.13) 



(2.14) 



from a scalar boson (Fig. 2.1 B). The Standard Model has a top quark with k ~ 1, but no scalars other than 
the Higgs boson itself. lSm L h is thus of order M 2 . h , and we need a finely tuned cancellation in the right hand 
side of the formula 



nr h (physical) = m\ (bare) + Am h 



(2.15) 



to realize m\ ~ 10 4 GeV 2 . 

How should we settle this unnaturalness? One possibility is to assume that Mhi g h is near 100 GeV, and 
that our quantum field theory cannot be applied above this scale. Then the tuning is not so fine, and moreover, 
we need not care quantum corrections at any higher scales. 

A more beautiful solution is provided by the supersymmetry (SUSY) flDl - The SUSY, a symmetry 
between fermions and bosons, supplies two scalar bosons for a respective Dirac fermion. In other words, 
a Weyl fermion is accompanied by a complex scalar field under the SUSY, and vise versa. Moreover, the 
SUSY guarantees that the scalar partners <f>\ 2 for a Dirac fermion iff that has a coupling k to the Higgs boson 
do couple to the Higgs boson with the exact coupling k s — \k\ 2 . Therefore, the "superpartners" cancel the 
quadratic divergence as 



\k\ 2 



L 

y^-j M high + <5(log Mhi g h) ~» <5(log Mhigh). 



(2.16) 



We will discuss this powerful hypothesis in Chapter|4] 
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Section 2.2 The Muon g - 2 
2.2.1 Foundation 

The g-value of the magnetic moment is one of the most famous and long investigated quantities in frontier 
of physics. It is defined as a ratio between the spin magnetic dipole moment fi spin and the spin vector S of a 
particle as 

\e\ 

/'spin = qgmS where w ■= ^— > ( 2 - 17 ) 

and expressed as a term in the Hamiltonian H as 

H?-ix spiD -B = -(qg t x B )S-B, (2.18) 

where m is the mass, and q\e\ is the electric charge, of the particle. 

The g-value for electrons and muons, g e and g M , are predicted as g e — — 2 under the quantum 
mechanics. However, radiative corrections shift the values slightly. The shift is known as the anomalous 
magnetic moment g - 2, or a := (g - 2)12. 

Historically, the discrepancy was known through the observations of Lande ^-factor gj of atoms. In 

1947, Schwinger calculated the one-loop level QED contribution to the electron magnetic moment as a e = 
a/2n = 0.00116 lfT8l . where a w 1/137 is the fine structure constant at the low-energy scale. Then, in 

1948, Kusch and Foley precisely measured the difference as a e = 0.00119 with comparing the gj values of 
gallium, indium and sodium fl9l . 

Currently the most precise measurements are achieved by the Harvard group for the electron at the sub- 
ppb level [20 1, and by the muon g — 2 collaboration at Brookhaven National Laboratory for the muon at the 
sub-ppm level. Theoretical calculations at similar precision are achieved, as exemplified by the five-loop 
level calculation for the QED contributions. The values are summarized as 

o e (exp) = (11 596521.8073 ±0.0028) xl0~ 10 (0.24ppb) (20), (2.19) 

a e (SM) = (11 596521.8178 ±0.0077) xl0~ 10 (0.67ppb) ED, (2.20) 

a^(exp) = (11 659208.9 ±6.3) xl0~ 10 (0.54ppm) |22||23), (2.21) 

a^SM) = (11 659 182.8 ±4.9) xl0~ 10 (0.42ppm) |24||25]|26l|27l|28l, (2.22) 

and the discrepancies between experiment and theory are 

a e (exp - SM) =(-1.06 ± 0.82) xl0~ 12 (1.3cr-level), (2.23) 
a^exp - SM) = (26.1 ± 8.0) xl0~ 10 (3.3<r-level). (2.24) 

The Standard Model predicts the experimental values very well, but there exists 3<r-level discrepancy in the 
muon g — 2. This is the muon g - 2 problem. 

Now let us see the muon g - 2 prediction under the Standard Model in detail. The g-factor of the muon 
magnetic moment is expressed in the QFT language as 

g = 2[F 1 (0) + F 2 (0)], (2.25) 
where Ffy 2 ) are the form factors of the p-p-y vertex function P 1 , which is described in Fig. : 



2.2 



A: 



PV, p) = yKFdq 2 ) + -^F 2 (q 2 ). (2.26) 
Lm 

(See, e.g., Ref. l29l for a review.) 

At the tree level Fi(0) = 1 and /*2(0) = 0, and thus g = 2. Quantum corrections modify these factors, but 
since the correction to F\(q 2 ) is sunk into the renormalization of e, the muon anomalous magnetic moment 
can be expressed as 

o„ = F 2 (Q). (2.27) 
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(A) (Bl) (B2) 

Figure 2.2: (A) The Feynman diagram which contributes to the g-factor of the muon magnetic moment. The 
tree level diagram corresponds to the classical result g M = 2, and quantum corrections, i.e., loop diagrams, 
yield deviation from 2. p, p' and q — p' — p are momenta of the particles. (Bl) The hadronic vacuum- 
polarization contributions to the muon g - 2. (B2) The hadronic light-by-light contributions to the muon 
g -2. The gray circles denote hadronic loop diagrams. 



2.2.2 Standard Model evaluation on the muon g - 2 

The Standard Model prediction of the muon g - 2 is summarized as: 



a^(QED) = (11658471.8951 


± 0.0080) x 10" 


10 


El, 


(2.28) 


<V(EW) = ( 


15.4 


±0.2 


)xl0" 


10 


nana, 


(2.29) 


^(HVP-LO) = ( 


692.3 


±4.2 


)Xl0" 


10 


eh, 


(2.30) 


( 


694.91 


±4.27 


)Xl0" 


10 


ED, 


(2.31) 


^(HVP-HO) = - ( 


9.84 


±0.07 


)Xl0" 


10 


ED, 


(2.32) 


^(HLbL) = ( 


10.5 


±2.6 


)Xl0" 


10 


ESI. 


(2.33) 


( 


11.6 


±4.0 


)xl0" 


10 


MM- 


(2.34) 



For the "HVP-LO" and the "HLbL" contributions two values are cited as a reference. 
The respective categories are defined as follows: 

QED contribution comes from the diagrams only with leptons and photons. It is calculated analytically 
up to the three-loop level (of order a 3 ), and recently an automated computation of the Feynman di- 
agrams finished its five-loop level (of order a 5 ) calculation I12T1 l24l . The uncertainty of the QED 
contribution, dominated by the uncertainty in measurement of a, is much smaller than those of the 
other contributions. 

Electroweak contribution is from the diagrams with Higgs, W and/or Z boson but without gluons. The 
contribution is evaluated at the two-loop level with including leading log three-loop effects [25 , 26 1. A 
hadronic loop uncertainty of ±0. 1 xl0~ 10 , and an uncertainty of ±0.2 xl0~ 10 from the "then-unknown" 
Higgs boson mass are considered, but the latter corresponds to the mass range of 114GeV < m/, < 
250 GeV allowed in those days, and thus is considered to be improved. 

Hadronic contribution is that including QCD interaction. This contribution is separated into two types: 



the diagrams of Fig. 2.2 Bl, called the hadronic vacuum polarization (HVP) contribution, and those 
of Figs. |2.2| B2, the hadronic light-by-light (HLbL) contribution. The lowest order contribution of the 
HVP, or HVP-LO, enters at the order a 2 , while higher order contributions (HVP-HO) are of order a 3 
as well as the HLbL. 

The HVP contributions cannot be calculated directly, and are evaluated through the dispersion relation 
using experimental data of the cross sections cr{L (e + <?~ — > y* — > hadrons). Due to limited accuracy 
of those data, the hadronic contributions are the dominant source of the uncertainty of the muon 
g - 2 prediction; especially there lies a disagreement among experimental data for the 2^-channel 
region. Several results from different collaborations are currently available [27, 30, 33], two among 
which Ii27ll30l are quoted as a reference. 
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The evaluation via the dispersion relation is simply summarized as the following equation: 
a„(HVP-LO) = f ds of 01 (e + e~ -> y* -> hadrons) K(s), (2.35) 

where K(s) is a kernel function (See, e.g., Ref. |34|). Note that cr° M is evaluated with final 
state radiations but without initial state radiations and vacuum polarization corrections. 

The HLbL contributions cannot be calculated directly, nor be evaluated with experimental data. Lattice 
calculations l35ll are expected, but currently low-energy effective theories [28 1 3T|[32) are exploited. 

Summing up all the above contributions, we obtain the Standard Model expectation. For example, if we 
combine the values from Refs. E4ll25l 26 27, 28 1 with errors in quadrature, the muon g - 2 is predicted as 



fy(SM) = (11 659 182.8 ± 4.9) xl(T 10 , (2.36) 



and we face a 3<x-level discrepancy. 



Several models have been invented to solve the muon g - 2 anomaly. The SUSY can be a solution again, 
where the superpartners, heavy but having ordinary couplings to muons, yield loop level contributions to shift 
the predicted value. Another candidate is a hidden photon |36, 37 1, an extra U(l) gauge boson with a feeble 
mixing with the photon. In this case, the mixing is as small as 10~ 4 , but if the mass of the hidden photon is 
of order 100 MeV, the theoretical value of the muon g - 2 is shifted enough to match the experimental result. 
This model is out of scope of this dissertation; we will concentrate on the SUSY, and see how the SUSY 



solves the muon g - 2 problem in Sec. 4.3 



Section 2.3 Slight Mismatch of the Gauge Coupling Constants 

As is mentioned when we discussed the hierarchy problem in Sec. |2.1| the Standard Model is considered 
as a low-energy effective theory of more fundamental theories, where we expect that some features of the 
Standard Model are "unified" with beautiful and sophisticated manners. An example is the grand unification, 
the unification of the three forces. 

Several models are proposed to realize the grand unification, and the most famous ones are the SU(5) 
grand unification theories (SU(5)-GUTs) Il38ll . Since the SU(5) gauge group includes Gsm as a subgroup, 
this scenario is very promising, and has been studied for several decades. SU(5)-GUTs have only one force 
of SU(5). Let us call the gauge coupling gs- The gauge group is expected (and assumed) to confront 
spontaneous symmetry breaking and to break down into three forces at a scale higher than (9(100) GeV: 
SU(3) with a gauge coupling g sm) = g 5 , SU(2) with g sm) = g 5 , and U(l) with g m) = ^/3J5g 5 . 

Here we face one problem. Our three forces do not have the strengths of {g s , g2, gy) = (<7s, 05, V3/5#5), 
but instead ~ (1.2,0.65,0.36) at the electroweak symmetry breaking scale mz- Fortunately, however, the 
gauge couplings depend on the scale. The dependence is embedded in the renormalization group equations 
(RGEs) of the gauge coupling. In the MS scheme [39|, the RGEs are given as |40| 

\ 1 / 199 2 27 2 44 2 17 2 1 2 3 2 \ 3 mvn 
^ J + (iW I 50* + W 02 + 5" 03 " 10^ " 2 yb ~ 2^P' (23?) 

1 1 9 -A + + 12 03 - \vt - Id - \y 2 A & (2-38) 





1 


dlogg 


16tt 2 


&gi 


1 


dlogg 


16tt 2 




1 


dlogg 


16tt 2 



(167r 2 ) 2 \lO yi 6 * 2 ^ 2 J ' 2 T 



1 /ll 2 



g\ + -gi - 26gi - 2yf - 2y L b g>, (2.39) 



?3 / (16^ 2 ) 2 \lO tfl 2* 
at the two-loop level, where Q is an energy scale of evaluation and we have defined 

(g3,g2,gi) ■= g s ,gi, yl jgr 



(2.40) 
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Q [GeV] 

Figure 2.3: The renormalization group evolution of the gauge coupling constants g^, g^ and g\ (black lines, 
from top to bottom) under the Standard Model, together with the evolution of y, (red line). Evaluation is done 
at the two-loop level for the gauge couplings. The Yukawa couplings entering the evaluation are evolved 
with the one-loop level RGEs. Several remarks on this figure are documented in Appendix |4.B| 



Using these RGEs together with one-loop level RGEs for the Yukawa couplings, which are given as 

,9,3 





1 


! 17 


dlogg 


16k 2 


^ 20 9 


dy b 


1 




dlogg 


I6n 2 


di/r 


1 




dlogg 


16k 2 



->yl + ylW 



2 3 2 ^ t 2 



Vb, 



g\ + 3yj + 3y 2 b 



(2.41) 
(2.42) 
(2.43) 



we can evaluate the values of the gauge couplings at any energy scales. 

The result is shown in Fig. 2.3 The couplings approach to each other, and gather closely at Q ~ 
10 I5 GeV, but complete matching is not achieved; there is a slight mismatch of the gauge coupling con- 
stants. 

How can we settle this situation? First we must keep in mind that SU(5)-breaking effect might change 
the runnings near 10 15 GeV. Since all what we need is just a slight shift, we can optimistically hope that the 
gauge coupling unification is realized with help from the effect, called "threshold corrections." 

Another solution is, again, the SUSY. Since the SUSY introduces many fields, the renormalization group 



running is modified. We here postpone the discussion in Sec. 4.4 where we will see the SUSY performs 
fabulous miracle. 



Section 2.4 Concluding Remarks 



In this chapter, the Standard Model and related issues are introduced. 
Sec. 



We saw the hierarchy problem in 



2. 1 the muon g - 2 problem in Sec. 2.2 and the unification problem in Sec. 2.3 



Miraculously the SUSY can solve these three problems. We have now enough reasons to believe Nature 
adopts the SUSY. However, on the contrary, we have no signature of the SUSY although it was optimistically 
expected to be found at the Large Hadron Collider. We have indirect evidences but no direct evidences. How 
should we understand this, somewhat strange, situation? 

Jurists have the principle of "in dubio pro reo". We scientists respect a similar principle; experimental 
evidence is the most important for the theory, and any beautiful theories could not win their validity without 
signatures in experiments. The SUSY is, therefore, just a hypothesis yet. 



12 



Dissertation / Sho Iwamoto 



Then, what we can do now is to fully utilize the LHC to determine whether the SUSY exists or not, and 
therefore what is important is to examine LHC SUSY searches without optimism and to invent more efficient 
way to detect the SUSY. From this viewpoint, the next chapter is dedicated to an overview of the ATLAS 
detector, one of the experiments at the LHC. The SUSY comes on the discussion after this excursion. 
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Chapter 3 

The ATLAS Experiment 

The ATLAS experiment ATI l42ll is a general-purpose detector for the Large Hadron Collider (LHC). With 
utilizing the detector, the ATLAS collaboration recorded a great milestone of the Higgs discovery in July 
2012 [2 1 together with the CMS collaboration They have also searched for signatures from models 
beyond the Standard Model such as the supersymmetry (SUSY). 

In Chapter [5] we will apply a result of a SUSY search from the ATLAS experiment with detailed dis- 
cussion on the ATLAS detector, such as consideration of lepton detection efficiency. In preparation for the 
discussion we here briefly review the detector instrument and how physical objects (such as electrons and 
muons) are detected. 

The CMS detector 11431 is the other general-purpose detector for the LHC. Its apparatuses are slightly 
different, but basic ideas of design and general strategy for the object detection are similar to those of the 
ATLAS detector. Therefore, although we simply concentrate on the ATLAS detector in this dissertation, the 
discussions can to a considerable degree be applied to the CMS detector. 



Section 3.1 Overview 
3.1.1 Coordinates 

For later convenience the definition of coordinates in the ATLAS detector is documented at first. 

The ATLAS collaboration uses a right-handed coordinate system with the origin at the nominal collision 
point. From the collision point the x-axis directs the center of the LHC ring, y-axis points upwards, and 
z-axis is defined along the beam axis. Spherical coordinates (r, 6, <p) is defined as usual: r := yjx 2 + y z + z 2 , 
cosO :- z/r, and tan</> := y/x. The pseudorapidity 77 is defined as 77 := - log tan(#/2). The transverse 
direction corresponds 77 = 0, while 77 = ±00 point the beam direction. 

Momenta and energies in the detector are usually expressed as "transverse momenta" p T and "transverse 
energies" Et, which is defined as p T = (p x , p y , 0), pj = llpJI = p sin 8 and Ej — E sin 6. Angular distances 
are often expressed in terms of AR := ^j(JS.rf) 1 + (A(f>) 2 . 



3.1.2 Detector 

The ATLAS detector [41] consists of three apparatuses: trackers, calorimeters, and muon spectrometer (MS). 
Fig. |3.1| provides an overview of the detector system, where we can find the trackers at the very center, the 
calorimeters as the brown apparatus surrounding the trackers and the gray surrounding the brown, and the 
MS as the outermost part. In the figure magnets are also drawn, which provide strong magnetic field in the 
ATLAS detector. A simplified cartoon is shown in Fig. |3.2| for a better understanding. 

The trackers, or the inner detectors (ID), are the innermost apparatus surrounding the LHC beam pipe 
cylindrically. Charged particles flying through the detectors provide many hits, where actually they ionize 
detector materials, on the layers of the trackers, and the hits are reconstructed as charged tracks. The ID 
itself has also layered structure; from inner to outer, the pixel detector, the semiconductor tracker (SCT), and 
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Figure 3.1: An overview of the ATLAS detector and its subdetectors. The innermost layer is the tracker, 
which is composed by the pixel detector, the semiconductor tracker (SCT), and the transition radiation 
tracker (TRT). The calorimeters surround the tracker, and the muon chambers (MS) is installed in the outer- 
most sector. This figure is cited from Ref. ED. (ATLAS Experiment ©2008 CERN) 
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Figure 3.2: An extremely-simplified (and thus inaccurate) version of Fig. 3.1 The red line in each figure 



illustrates an imaginary trajectory of a particle produced at the collision point. The upper figure shows the 
approximate scale of the detector, and the lower is a cartoon describing the detector system. The detectors 
are roughly separated to the "barrel" region (not shaded) and the "end-cap" region (shaded). 
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Figure 3.3: An overview of the ATLAS inner detectors (ID), which can be found at the center of Fig. |3.1| 
The pixel detector, the semiconductor tracker (SCT), the transition radiation tracker (TRT) are respectively 
indicated. See Fig. |3.4| for the actual installation of the detectors. This figure is cited from Ref. |41|. 
(ATLAS Experiment ©2008 CERN) 
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Figure 3.4: A cut-away view of the ATLAS inner detectors (ID). The collision point locates in the left- 
bottom corner of this figure. We can see that each detector is separated to the barrel part and the end-cap 
part. The trackers are installed within I77I < 2.5, but the TRT coverage is restricted to < 1.9. This figure 
is cited from Ref. E). (ATLAS Experiment ©2008 CERN) 



16 



Dissertation / Sho Iwamoto 



Tile barrel Tile extended barrel 




Figure 3.5: A cut-away view of the ATLAS calorimeters. Detailed explanation is found in the text. This 
figure is cited from Ref. glj. (ATLAS Experiment ©2008 CERN) 



the transition radiation tracker (TRT) are installed. In Fig. 3.1 the pixel is installed inside the very central 
yellow tube, and the SCT surrounds it. They as well as the TRT are clearly observed in Fig. 3.3 a detailed 
view of the ID. In addition, to see the exact size of each inner detector, a cut-away view of the ID is provided 
in Fig. |3.4| 

The calorimeters, surrounding the ID, are composed by the electromagnetic calorimeter (ECAL), the 
hadronic calorimeter (HCAL), and the forward calorimeter (FCAL). A detailed cut-away view is provided 
in Fig. 3.5 where we can clearly find the ECAL as "LAr electromagnetic barrel" and "LAr electromagnetic 
end-cap (EMEC)," the HCAL as "Tile (extended) barrel" and "LAr hadronic end-cap (HEC)", and the FCAL 
as "LAr forward (FCal)." 

The ECAL is the inner part of the calorimeter, which consists of a barrel part ("LAr electromagnetic 
barrel") for \rj\ < 1.475 and two end-cap parts ("LAr electromagnetic end-caps") for 1.375 < I77I < 3.2. They 
all are sampling calorimeter with liquid argon (LAr) and lead. Electrons and photons mainly interact with the 
ECAL, and lose almost all of their energy to be observed as electromagnetic jets in the ECAL. The energy 
loss of electrons is through Bremsstrahlung, while photons have various interactions such as e + e~ pair- 
production, photo-electric effect, and Compton scattering. Note that the muon does not stop here because 
its energy is too small to cause significant Bremsstrahlung; the energy loss AE/dx with Bremsstrahlung of a 
particle whose mass is m is approximately proportional to E/m 2 , and thus, E > 1 TeV for muon is required 
to trigger Bremsstrahlung. 

The HCAL of the ATLAS detector is a sampling calorimeter for hadronic particles. It is with large dense, 
and hadronic particles interact with the detector material to deposit their all energy and stops; finally they 
are observed as hadronic jets. Iron-scintillating-tile system is used to cover the barrel region: the apparatus 
for \t]\ < 1.0 is called "tile barrel", and for 0.8 < I77I < 1.7 are "tile extended barrel" as shown in Fig. 3.5 For 
the end-cap region copper-LAr system is used, which is displayed as "LAr hadronic end-cap" in the figure. 

The very forward regions of 3.1 < |?7| < 4.9 are covered with the FCAL. Copper-LAr system for 
electromagnetic objects and Wolfram-LAr system for hadrons are utilized. 

The MS, the outermost apparatus of the ATLAS detector, is another tracker for muons. All the other 
(known) particles but muons and neutrinos are already trapped in the previous detectors, and here momenta 
of muons are measured. 

In addition to the above detectors, the ATLAS experiment has three smaller detector systems in a very 
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forward (and backward) region: LUCID (luminosity measurement using Cerenkov integrating detector), 
ZDC (zero-degree calorimeter), and ALFA (absolute luminosity for ATLAS). They lie + 17m, +140 m, and 
+240 m from the collision point, respectively. The main purposes of the LUCID and the ALFA are to mea- 
sure the integrated luminosity and to provide online monitoring of the instantaneous luminosity and beam 
conditions [411 . The ZDC is mainly utilized to detect very forward neutrons in heavy-ion collisions BP . 

* * * 

The performance of the ATLAS detector is fully exploited within \tj\ < 2.5, or 9° < 8 < 171°. The 
coverage of the detectors is approximately summarized as: I77I < 2.5 for the trackers, \rj\ < 3.2 for the ECAL 
and the HCAL, and \rj\ < 2.7 for the MS. The FCAL is installed at 3.1 < \rj\ < 4.9. Only the central region is 
utilized for usual searches for models beyond the Standard Model. 



3.1.3 Triggering 

The nominal bunch-crossing rate of the ATLAS detector is 40 MHz, which means that enormously vast 
numbers of collisions occur at the center of the ATLAS detector. Not all of them can be recorded to the 
storage due to limitations of processing speed and storage capacity. Therefore, a system of event selections, 
called triggering system, is installed in the data-taking system. The selections are designed to pickup events 
of our interest to reduce the output rate to ~ 200 Hz (=* 300 MB /s), which is acceptable by the data-taking 
system [42|. 

Several selections are installed inclusively; i.e., an event passes the trigger and is stored for the physics 
analysis if it fulfills at least one of the trigger requirements. A famous selection is that requiring hard jets 



plus missing energy; for example, in the analysis which we will examine in Sec. 5.6.2 events which passed 
the trigger requiring missing energy of^j > 100 GeV and at least one jet with pj > 80GeV are utilized. 

Note that the criteria depend on the instantaneous luminosity of the LHC; in the run of 2012, the bunch- 
crossing rate was < 20 MHz, and thus the criteria were looser than those designed for the nominal rate. 



Section 3.2 Object Identification and Reconstruction 

The detectors are installed to detect particles produced in collisions, that is, in fact, to identify their species 
and to reconstruct their momenta and energy. Detection is always accompanied by efficiency, resolution, and 
fake rate. These are studied at first in silico, i.e., in Monte Carlo simulations, and now are measured in situ. 

We will briefly review those issues on identifications, especially for jets, electrons, and muons in detail, 
which are important in our numerical evaluations of the SUSY model discussed later. 



Hadron 

Hadronic particles are captured in the HCAL and produce hadronic showers to be observed as "jets." In 
the calorimeter, which consists of many cells, measured is how much energy is deposited in the respective 
cells. Thus we have to reconstruct jets from such information. Many algorithms, called jet algorithms, 
are invented to this end. Recent ATLAS SUSY searches usually utilized the anti-£, algorithm B4l with a 
distance parameter of 0.4. 

In the actual experiment the energy of the jet is calibrated to match the result from their full detector 
simulations, but it is beyond the scope of this dissertation. 



fe-tagging 

The b-quark always appears in new physics searches at the collider experiment in these days; the top quark 
always produces a b-quark, and the dominant decay branch of the Higgs boson is expected as h — > bb. 
Therefore, it is very important to identify hadronic particles which originate in /?-quarks. This is called 
"/7-tagging," and several algorithms are invented to this end. 
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One famous method for ^-tagging is the secondary-vertex method, which exploits a characteristic feature 
of the Z?-mesons that their lifetime is slightly longer than that of other mesons. It is about 1.5 ps, and thus the 
flight is a few millimeters, which is observable with the ID [45 1. 

^-tagging algorithms are applied to the reconstructed jets, and decisions are made that they are fo-jets or 
not. Here the efficiency and the fake rate appear. It is important that the fake rate for jets containing (or, 
originating in) c-quark is remarkably worse than that for the light-jets (jets without c- and fo-quarks). 

Recent analyses by the ATLAS collaboration utilize the MV1 algorithm, a neural-network-based al- 
gorithm whose input is taken from several stand-alone ^-tagging algorithms. As is expected, this neural- 
network-based algorithm generally yields better efficiency |46| and better jet rejection Il46ll47l . 



Electron 

Electrons produced at the collision point leave tracks in the trackers, and provide shower to stop in the ECAL. 
Thus identification and reconstruction of electrons utilizes the trackers and ECAL. Since similar signatures 
are provided by charged hadrons, it is important to reject jets faking electrons, but then we face a difficult 
trade-off between identification efficiency and rejection power. As a solution several criteria are defined for 
electron identification. 

The recent techniques for electron reconstruction are summarized in Ref. ||48l . In the document three 
criteria for the central region and two for the forward region are defined. 

For electron identification in the central region with the tracker coverage, \r/\ < 2.47, "loose", "medium" 
and "tight" selections are provided with an expected jet rejection of 500, 5000 and 50000, respectively. 
The "loose" selection employs only the calorimeter information. The "medium" selection further requires 
a track-cluster matching, and track quality provided by the pixel and the SCT In the "tight" selection the 
criteria of the "medium" are tightened, the information from the TRT detector is used, and discrimination 
against photon conversion is employed. 

For the forward region of 2.5 < I77] < 4.9, the region without trackers, "forward loose" and "forward 
tight" selections are defined, which utilize only the calorimeter information Electrons in this region are 
usually not used in searches for models beyond the Standard model. 

The electron reconstruction efficiency and the momentum resolution are publicly reported in Ref. 1481 * 1 
with the tag-and-probe method based on - 7 TeV collision data corresponding to the integrated luminos- 
ity of 40pb~', which were taken in the year 2010. Here one should note that the efficiency for an electron 
f eiec j s (j ecom p 0sec i j n to two factors as 

e dec = areco • em- (3.1) 

The first term ar reco is the efficiency for the cluster in the ECAL from an electron to be reconstructed well, 
and the second term enj is that for an electron to pass the selection criteria discussed above. a reC o is measured 
with Z — > ee events, and en) is with Z — > ee for pi e [20,50] GeV, W — > ev for pi e [15,50] GeV and 
J /if/ -> ee for p T s [4, 20] GeV. 

The Bremsstrahlung effect with materials before the ECAL is not negligible for electrons. To correct 
the energy loss, a technology called the Gaussian sum filter is utilized in some of recent analyses, e.g. in 
Ref. 11491 . which provides better performance in the electron detection [50 ] . 

The ATLAS collaboration seems to have re-optimized the selections and defined "loose++", 
"medium++", and "tight++" 1511 , whose definitions are hardly found in public documents. We do 
not discriminate the re-optimized selection with the original ones for simplicity. 

These reports are summarized on an web page* 2 , where preliminary plots of the efficiency and the reso- 
lution measurements are also available. 



1 This is the latest public report (conference note) at the date of November 20 1 2. 

"^https : //twiki . cern . ch/twiki/bin/view/AtlasPublic/ElectronGammaPublicCollisionResults for electrons, and 
https : //twiki .cern. ch/twiki/bin/view/AtlasPublic/MuonPer£ormancePublicPlots for muons, as of November 2012. 
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Muon 

Muons in the ATLAS detector, produced at the center, fly through the ID, pass through the calorimeters 
with a little energy deposits, and reach the muon spectrometer (MS). Muons are reconstructed with the MS, 
which is essential for muon identification and reconstruction. The ID are also used to increase resolution 
and efficiency. 

Recently the ATLAS collaboration reported three methods of the muon identification [52, 53 1. Stand- 
alone (SA) muon is the muon which is reconstructed with utilizing only the MS, i.e., without using the 
ID. Their direction and momenta are determined with the trajectory in the MS, taking energy losses in 
the calorimeters into account. Combined (CB) muon and segment tagged (ST) muon are, meanwhile, 
reconstructed with utilizing both the ID and the MS. In the former case the identification is respectively 
done in the MS and the ID, and a combination of an MS track and an ID track is employed. A successful 
combination results in a CB muon. In the latter case the identification is done in the ID. The ID track is 
extrapolated to the MS, and a successful association with MS track segments is reconstructed as a ST muon. 

CB muons are with the highest purity, and thus usually utilized. ST muons are sometimes used supple- 
mentarily to compensate the inefficiency of that for CB muons. 

The muon reconstruction efficiency is publicly reported in Ref. 1 53 1 and Ref. Il54l * 3 . In Ref. l53l the 
efficiency for pj > 20GeV is measured with Z — > pp events. The efficiency for pi e [2, 10] GeV is reported 
in Ref. If54ll . where J/if/ — > pp events are used. Both measurements utilize the tag-and-probe method, and 
are based on the yfs = 7TeV collision data corresponding to the integrated luminosity of 40pb~\ which 
were taken in the year 2010. 

The muon momentum resolution is reported in Ref. |52]* 4 , which is based on Z — > pp andW — > pv 
events in the 20 1 data of 7 Te V and 40 pb~ 1 . 

These reports are summarized on an web page* 2 , where preliminary plots of the efficiency and the reso- 
lution measurements are also available. Especially several plots for the efficiency measurement at the 8 TeV 
LHC can be found as a rapid communication |55 1. 



Section 3.3 Concluding Remark 

With the detector and the strategy we have seen in this chapter, the discovery of the Higgs boson was 
achieved. However, here another factor which helped the discovery should be emphasized. It is provided by 
the effort of the collaboration for the LHC accelerator. 

The LHC ran at the energy of Eqm = 7 TeV in 2010 and 2011, and data corresponding to an integrated 
luminosity of 5.3 ftr 1 were recorded at the ATLAS experiment. Then in 2012 the energy was increased to 
Ecm = 8 TeV, and data of 21.7 ftr 1 were obtained. Especially in 2012 run, the collision rate was increased 
to 20 MHz, and the peak luminosity of 7.73 x 10 33 cm _2 s _1 is recorded with fully utilizing the accelerator. 
With these great utilizations of the LHC, the discovery of the Higgs boson was achieved. 

On 17 December 2012, the LHC was shut down to prepare for collisions with Ecm = 13 TeV and 14 TeV. 
In the 13 TeV run expected to start in 2015, the collision rate will be increased to the nominal 40 MHz, which 
will provide larger instantaneous luminosity of ~ 1 x 10 34 cirT 2 s~' to allow us to expect data corresponding 
to (9(100) ftr 1 in the 13-14 TeV runs.* 5 

Now we have to wait for three years, but now, contrary to that after the accident in 2009, we have data, 
which are enough to allow us to find the "tail" of the physics beyond the Standard Model buried inside them. 

In the next chapter, we will review the MSSM as a promising candidate for such theories, and discuss its 
current status with examining the data. 



* 3 These are the latest public reports (conference notes) at the date of November 2012. 
4 This is the latest public report (conference note) at the date of November 2012. 

5 Also in order to obtain much more data the HL-LHC (High Luminosity Large Hadron Collider) is proposed, where an instantaneous 
luminosity of ~ 3 X 1 34 cirT 2 s~ 1 and an integrated luminosity of ~ 3000 ftr 1 are expected. 
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Chapter 4 

The MSSM and Its Current Status 



Section 4.1 The MSSM 

The minimal supersymmetric standard model (MSSM) |[TTl[T2l[T3l is the minimal supersymmetric extension 
of the Standard Model. It is characterized by imposed gauge symmetries and field content, as well as the 
Standard Model. The gauge symmetry of the Standard Model is the same as the Standard Model: 

Gmssm = SU(3) strong x SU(2) weak x U(l)y. (4.1) 

The field content is shown in Table |4.l| Here, unlike the Standard Model, we need two Higgs doublets H u 
and H&. This is because the SUSY severely constrains the Lagrangian to forbid a single Higgs doublet to 
have Yukawa interactions both with the up-type quarks (Q and U) and down-type quarks (Q and D). Note 
that the anomaly cancellation condition, in the consideration of which fermionic partners of the Higgs boson 
must be taken into account, is satisfied with these two doublets. 

The symmetry and the field content lead us to the following superpotential of the MSSM: 

W m = fiH u H 6 - (YJijHnQiUj + (Y&jHiQtDj + (X^jH^Ej 

1 _ , - 1 „ ( 4 - 2 ) 

+ KiH u Lj + -A ijk LiLjE k + A uk LiQjD k + -A ij k UiDjD k . 

In the terms in the second line of the expression, however, the baryon number B and the lepton number 
L are not conserved. This is not only inconsistent with the fact that B- and L-violations are not observed 
yet in experiments, but also makes protons decay in a very rapid rate via, for example, the Feynman di- 



agram described in Fig. 4.1 The decay rate of the proton is severely constrained by Super-Kamiokande 



experiments |56|. For this channel, the constraint is 

m 5 \A' A" '" 



FY +\ h' 3" |2 m pr°t°n rm'MU /lTeV\ / c ,, , A 33 Y l (A n\ 

T(p^ne)~\A m A n2 \ -^ = 29xlQ . 20yr (— ) <(8-2xl0 yr) , (4.3) 

which results in the constraint on the coupling constants of U' 112 ^j' 12 | % 10~ 27 . 

In order to solve this unnaturalness, we usually install the conservation of the R -parity |[T3l into the 
MSSM, which is a discrete Z 2 symmetry defined as Pr := (_i) 3B L + 2j . Here s is the spin of the particle. 
The exact conservation of the 7?-parity restricts the superpotential as 

W = nH»H A - (Y u ) u H u QiUj + (Y d ) u H^Dj + (Y e ) u H 6 L,Ej, (4.4) 

and now the baryon number B and the lepton number L are conserved at the classical level. (Note that the 



sphaleron process, mentioned in Sec. 2.1 ) Conservation of B and L at the tree level is achieved accidentally 
in the Standard Model, but in the MSSM we have to impose it by hand. We in this dissertation assume the 
/^-parity conservation, and use this superpotential (|4.4[). 
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Table 4.1: The field content of the MSSM. In the leftmost column the superfield notation is used, and the 
rightmost two columns describe the included fields of the superfield: complex scalar fields for the spin 
particles, Weyl spinors for the spin 1/2 ones, and vector fields for the spin l's. The gauge indices are 
omitted, while indices for the three generations are denoted as subscripts ;'. 



Matter and Higgs fields (chiral multiplet) 
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spin 1/2 


Qi 


3 


2 
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Figure 4.1: Feynman diagram of the proton decay caused by the interactions U\DiD 2 and L l Q l D 2 - To 
forbid this interaction we have to introduce the /^-parity. 
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This "MSSM with R -parity" provides a very nice explanation of the Dark Matter problem. Under the 
^-parity conservation the lightest particle among those having odd parity, actually which is the lightest 
supersymmetric particle (LSP), becomes stable. The LSP is considered as a promising candidate for the 
Dark Matter 1571. 

Strictly speaking, the ^-parity does not make protons stable. Besides the sphaleron process, the 
B and L are violated in higher dimensional interactions even in the presence of the ^-parity con- 
servation. Such terms can be forbidden by imposing, e.g., the baryon triality 6 3 |58 1 or the proton 
hexality [59 1 instead of the /^-parity, where these symmetries conserve B and L up to dimension 
five operators with holding the stability of the LSP. 

On the other hand, just for avoiding the proton decay problem, we do not have to forbid both the 
B- and L- violations. Protons remain stable if the baryon number is conserved, or the lepton number 
is conserved and the LSP is heavier than protons. In both choices the LSP becomes unstable, and 
difficult to be served as a Dark Matter candidate. 

Now we can write down the full Lagrangian of the MSSM, but it is still insufficient because we know our 
Universe is not supersymmetric. We have to introduce additional "SUSY-breaking (SUSY)" terms, which 
do not respect the SUSY, into the Lagrangian. However, we must be careful not to disgrace the MSSM. 
The hierarchy problem is solved in the MSSM because the SUSY guarantees the condition |fe| 2 = k s in 
Eqs. p,13[ ) and p,14| i. For this reason we usually introduce only the "soft" SUSY-breaking terms, or terms 
whose coupling constants have positive mass dimension. Then the condition \k\ 2 - k s does not violated even 
in quantum corrections because soft couplings do not appear in the renormalization group equations of the 
dimensionless couplings. 

The soft SUSY-breaking Lagrangian of the MSSM is given as 



(4.5) 



-Xsusr= \ {Migg + M 2 WW + M X BB + H.c.) 

-(a u )ijH u QiUj + {a d )ijHiQidj + {a e ) i jH i L i e j + H.c. 

+ |(4) fj Q*Qj + (m 2 L )..L*Lj + (m 2 D )SiUj + {m 2 ^..d i dj + {rr^,'X i l j 
+ [m 2 H H* u H a + m\H\H A + (bH u H d + H.c.)] , 

where M a and m 2 , are the (quadratic) masses for the gauginos and the scalar bosons (Higgs bosons, squarks 
and sleptons), a lud , e are called trilinear scalar couplings, and b is the off-diagonal mass term of the Higgs 
bosons. For later convenience, we here define the following parameters as usual: 



(A u )u :-- 



(a u )ij 



(A d ) u :-- 



(Y d )i 



(A e ) 



(Ye)i 



B 



b 



(4.6) 



In the MSSM we have redefined the phases so that the coefficient b, the vacuum expectation values of the 
Higgs bosons, and fermion masses are to be positive. Then the CKM matrix has one phase as is in the case of 
the Standard Model, and the phases are thrust to the SUSY-breaking parameters M a and a,. Nevertheless, we 
ignore these phases in this dissertation; in particular, M a are assumed positive throughout our discussions. 



Breaking a symmetry by hand is always accompanied by a new problem, and this is the case, too. Not 
only M a but also m 2 and a in Eq. ( |4.5| > have complex phases in general and induces CP-violation. Also 
those for quarks and leptons can have off-diagonal component, and large flavor violation may be yielded. 
These CP- and flavor violation may easily contradict with current observations; these are called the SUSY 
CP- and flavor problems. These problems originate in a fundamental question that what kind of mechanism 
generates the SUSY-breaking effects. Various models have been built, or are being built, to achieve a realistic 
SUSY-breaking. 

A promising model for SUSY-breaking is the gauge-mediated SUSY-breaking (GMSB) scenario l60l 
|6T1 . In this scenario the source of SUSY-breaking is "hidden" from the MSSM particles, i.e. there are no 
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direct interactions between the MSSM fields and the source of SUSY-breaking. Instead the existence of 
the source, or the hidden sector, is mediated with extra fields, called messengers, which interact with the 
MSSM particles via the Standard Model gauge interactions. Since the gauge interactions do not "know" 
CP-violating phases or the flavor indices, we can circumvent the SUSY CP- and flavor problems. This 
GMSB framework is adopted to the V-MSSM, introduced in Sec. [5] which we call "V-GMSB scenario." 



Section 4.2 Higgs Mass in the MSSM 

The Higgs sector of the MSSM is completely different from that of the Standard Model. We thus begin with 
the discussion on the SUSY from reviewing the Higgs sector and the electroweak symmetry breaking inside 
it. 



4.2.1 A brief review of tree-level result 

We start from the tree-level discussion The tree-level scalar potential related to the Higgs fields is given 
directly from the superpotential (|4.4|> and the SUSY-breaking Lagrangian (|4.5|l as 



(4.7) 



<L = (n 2 + <) (Kl 2 + K| 2 ) + (n 2 + <) (|< + K) + [b (fljflj - flSflj) + h.c] 
+ ^ (Kl 2 + K - Kl 2 - 1"/) 2 + f \khT + H°unr\ 2 ■ 

We choose the SU(2) basis not to break the electromagnetic symmetry, which leads us to 

Vg* = (N 2 + <) W + (N 2 + <) K - + H.c.) + (|// u f - |< ) 2 . (4.8) 



As discussed below Eq. (4.6) , we take as a convention > by redefining the phases of H^. Then at minima 
of V® H°H° d > 0, and thus we can take the convention under which the vacuum expectation values of 

and //9 are real and positive. 

We would like to derive the vacuum expectation values 

(H a ) = v n = ysin/3, (H d ) = v a = ocosjS, u»174GeV, (4.9) 

from this potential. Ideologically (v,/3) are extracted from the condition 

i} ^J=(iej=o, (4.io) 



flflO Higgs/ ^0 Higg 

but, since we actually know the value of v as 174 GeV, we use these conditions not to determine (u, tan/3) 
but to constrain the parameters of the MSSM (in 2 H ,m^j ,fi,b\ with treating (v, tan/3) as input values. Then 
we obtain 

m 2 Hu + \fi\ 2 -b cot/3 - ^m 2 z cos 2/3 = 0, m 2 ^ + |yu| 2 - fotan/3 + -/nfcos2/3 = 0, (4.11) 

where the Z-boson mass m| = + g^) v 2 /2is used, or equivalently 

^2j8=- —j, m 2 = I B, r - OTg - mjj - 1 |/x| , (4.12) 

<+<+2W ^1- sin 2 30 

as the constraints. 

1 Here we just incorporate a brief review. More detailed discussions can be found in, e.g., Ref. 1621 . 
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With the determined values of (m 2 H , rn^fj., £>), we can evaluate the mass terms of the Higgs sector. The 
mass matrix for the CP-even Higgs bosons, (h, H), is obtained to be 



where 



V 3 



Uh H)R a 



At™ Al^ 1 a 



eld 



m; 



(0) 



\jj\ 2 + m 2 Hu + -m|(l - 2 cos 2/3), 



At: 



(0) 



-Z? - -m| sin 2/?, 



and the tree-level masses are as 



(4.13) 



(4.14) 





2b 


m A 








sin 2/? 




1 




m l,H 


~ 2 


m|-i 






V 


m 2 H± 


= m A + m 



2 M 2 + + m 2 Hd 



+ m 2 A + J (m A - mj) + 4m|m^ sin 2 2/? 



(4.15) 

(4.16) 
(4.17) 



Here we have defined the well-known Higgs bosons* , h, H, A, H ± , and the Nambu-Goldstone bosons, G 
G ± , as 



(4.18) 



with the four rotation matrices 
/ cos a- sin a 

K„ = 

l-sina- cosff 



sin /Jo cos Bo 
-coSjSq sin/?o 



sin/? ± cosy6 ± 
-cosyS ± sinyS ± 



(4.19) 



The mixing angle a, which is traditionally chosen to be negative, is related to tan B at the tree level as 

tan 2a m 2 A +m 2 z 



tan 2/3 



m 2 A -m 2 z 



(4.20) 



The case where m A » m| is called "decoupling limit". In this limit the angles a and B are related as 



tan 2a 
tan 2/3 



(4.21) 



Let us check resulting features of the electroweak symmetry breaking. It is straightforward to check that 
the SU(2) gauge bosons successfully obtain masses: 



£ d |(d p - ig 2 W M - \ig Y B,) ff u | + \(d M - ig 2 W M + {igyB^H^ 

, 2i 2 \ (4.22) 

^-{d,h) 2 + V -[g 2 2 W^W u+ 9 4-Z^ 



which is completely the same as Eq. \2.9) , the Standard Model version. Meanwhile the fermion terms are 
slightly changed as 

-l d (Xu)ij Hy Li u ] Rj + (Y d ) u H° d du4ij + {Y ^i H W> e b 4 23 

(uY„y sin/?) ,7 cos/?) d^jd^j + (vY e ij cos/?) eue^./. 



2 Throughout this dissertation, the term "Higgs boson" in the context of the MSSM (and the V-MSSM) refers to the lighter CP-even 
Higgs boson 
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4.2.2 With one-loop level effective potential 

The above discussion is based on the tree level result. As will be mentioned in the discussion section 
(Sec. |4.2.4] >, quantum corrections are crucial to the mass of the Higgs boson in the MSSM, which currently 
is calculated up to the three-loop level accuracy [63, 64). Here, to examine the corrections analytically, we 
discuss the one-loop level correction with utilizing the effective potential method. 

We begin with the discussion how we should treat the effective potential in the calculation. What we 
should consider carefully is that the vacuum expectation value v ~ 174 GeV should not be changed even if 
the potential receives higher order corrections. That is, even the potential is modified as 
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must hold with the same value ofv » 174 GeV. Thus the constraints ( |4.11[ ) are modified to be 
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* + 2^U°/ + M / -^an^ + -^cos2^0, 



and the resulting parameters (m 2 H , m 2 ^, //, b} are different from those at the tree level. Here we have already 
imposed that tan f3 is not modified, or in other words, treated it as an input value. 

Here let us assume that (//, b) are kept unchanged, and the change of a, the mixing angle ofh and H, can 
be neglected. Then we can write down the parameters as 
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and the components of the mass matrix for the CP-even Higgs bosons are obtained as 
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et cetera, or simply evaluated as 
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Especially the mass of the lighter Higgs boson is modified to be 
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and if we take the decoupling limit, we obtain 
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Note that the tree level mass [m^] aee should be evaluated with absence of, and the correction terms should 
be with presence of, the correction AV. 



Here we set (fi, b) unchanged and neglected the change of a to obtain the equation l |4.30| >, and 
took the decoupling limit to Eq. {43T}. However, even if we instead set (m Hu ,m Hil ) unchanged, 
Eq. < |4.3 1 [ > can still be obtained at the decoupling limit. We should note that there lie several 
approximations, but since, anyway, corrections at higher levels are more important, we do not care 
these issues. 



4.2.3 MSSM Higgs mass at the one-loop level 

The effective potential of the MSSM at the one-loop level is known as, in the DR' scheme [65ll.* 3 
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(4.33) 



with renormalization scale Q, which is irrelevant in calculation of the Higgs mass. Note that the summation 
should be taken for all particles; i.e., one must not forget the color factor A^ c = 3. 

The dominant contributions come from the top-stop sector. The masses expressed as functions of (v a , t>d) 
are obtained from the mass matrices 
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Then through a bit tough calculation we can obtain 
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M top = Y t v a . (4.34) 

(4.35) 
(4.36) 
(4.37) 
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+ 0(g 2 2 ,g 2 Y ) (4.39) 
(4.40) 



where 



/)?: 



X, = A, - yUCOt/?, 



v ~ 174 GeV. 



(4.41) 



3 The DR' scheme is a modified version of the DR scheme |66 | in which effects of the mass of the e-scalars, appearing as the extra 
components of the vector bosons corresponding to the extra 2e dimension in the dimensional reduction methodology, are properly 
removed. More detailed description can be found in Ref. | 67 1. 
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4.2.4 Discussion — The 126 GeV Higgs 



The most important difference on the Higgs mass between the Standard Model and the MSSM is the Higgs 
quartic coupling. In the Standard Model the coupling, denoted by A in Eq. ( |2.6[ ), is an unknown parameter, 
and we can tune it to obtain the Higgs mass of ~ 126 GeV. On the contrary, under the MSSM, the quartic 
coupling is fixed by the SUSY, as is shown in Eq. ( |4.8| l. Therefore, the Higgs boson mass is restricted, and 
actually the mass has an upper bound at the tree level of [68] 



m\ cos 2 2/3 < m\. 



(4.42) 



As a result, in 1980's, the Higgs boson under the MSSM was considered lighter than Z-boson. Therefore, 
people considered that the Large Electron-Positron Collider (LEP), which was designed in early 1980's with 
£cm = 100 GeV and started operation in 1989, would discover the Higgs boson or reject the MSSM. In those 
days, the most promising channel for the Higgs search was therefore e + e~ — > Z — > he + e~ for the early stage 
and e + e~ — > Zh for Ecm <: 180 GeV, and this upper bound was referred in many reports from the LEP 
experiments. Ref. l69l from the DELPHI collaboration is an example, which was published in August 1990 
to report that the MSSM Higgs boson must be heavier than 28 GeV for all values of tan/3. 

However, Nature was not so simple. In the year 1990, several groups found that radiative corrections 
could raise the Higgs boson mass 117011711 . In 1991, before the exclusion limit reached to the Z-boson mass, 
those papers were published and the upper bound on the MSSM Higgs boson was loosened 1701171117211731 
1741 . For example, if one substitutes tan/? = 10, Ma = Ms — 1 TeV, X, = and A 2 = in Eqs. ( |4.16| l and 
( |4.40| i, the Higgs mass is evaluated to be 



/;?/, 



V(89.4GeV) 2 + (88.8 GeV) 2 = 126 GeV (X, = : no-mixing), 



(4.43) 



which could not be reached even after the upgrade of the LEP collider, i.e. with the LEP-II experiment. 
Higgs mass increase is enhanced at X, — + V6Ms due to the stop mixing, which is called to/, -max scenario 
or maximal mixing scenario. The Higgs mass becomes, with the same parameters and simply assuming 
A = 0, 



TO/, 



V(126 GeV) 2 + (82.1 GeV) 2 = 150 GeV (X, = + V6M 5 : maximal-mixing). 



(4.44) 



The above calculation was based on the approximation that only the top-stop sector was included. Other 
one-loop level contributions, such as threshold effects on the top-quark Yukawa coupling and corrections 
from bottom-sbottom and tau-stau sector, are important as well as higher-level loop corrections. The 
bottom-sbottom and the tau-stau one-loop level contributions are approximately given as l75l l76l 
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(4.45) 



(4.46) 



respectively. They give negative contributions of Ato/, 10 GeV to the Higgs mass. 

Conclusively, to realize the 126 GeV Higgs boson within the MSSM, the stop mass Ms must be 
(9(10) TeV in the no-mixing case and ~ 1-2 TeV in the maximal-mixing scenario (See, e.g., Ref. Il77l ). 



4.2.5 Discussion — The little hierarchy problem 

Another interesting feature of the Higgs sector in the MSSM is the so-called little hierarchy problem [78 1. 
This problem lies in Eq. ( |4.1 1[ ). The equation must be hold for the Higgs fields to have vacuum expectation 

values of v — Ji% + = 174 GeV, but it includes a cancellation among m 2 H , to 2 ^ and jj 2 . There we actually 
need a "little" tuning of order 
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where the MSSM parameters should be evaluated at the EWSB scale mz- 

This tuning is related to the mass of the stop. The parameter m 2 H receives the radiative correction of 



Am 



9 31? 9 A 

2 L^w (4.48) 



which can be seen from the renormalization group equations (RGEs) summarized in Appendix 4. A or 
more simply, from the effective potential approach as Am 2 H * (dAV/dH 1 ^ /(2v u ). Thus the tuning can be 
interpreted as 

tuning* 0.3% x(p^) ■ ( 4 - 49 ) 

We hope a lighter stop, or a lighter SUSY scale Msusy, to eliminate this little tuning. 

However, as we saw in the last section, the SUSY scale seems to be heavier than 1 TeV. How can we 
settle this collision? Is this inevitable price to pay for eliminating the fine-tuning? This is one of the hottest 
topics after the Higgs discovery, and currently under discussions. 



Section 4.3 Muon g - 2 in the MSSM 



As we saw in Sec. 2.2 the muon g - 2 has 3<x level discrepancy between its experimental result and the 
Standard Model theoretical value. The MSSM has extra contributions to the g - 2, and might be a solution 
to the g - 2 problem if the sign of the correction falls as the desired one, i.e., positive ll79ll . 

As we will saw later in Eqs. ( |4.53| l-( |4.57| i, all the extra contributions for a lepton are proportional to the 
mass squared of the lepton. Therefore, the SUSY contribution to the electron g - 2 can be neglected, which, 
thus, we do not consider afterwards. 

The dominant contributions come from the neutralino-smuon and the chargino-muon-sneutrino loop 



diagrams shown in Fig. 4.2 The calculation is straightforward, but summarized here as a reference. 



4.3.1 Formulae in mass eigenstates 

Relevant interactions of the MSSM are included in the Lagrangian as 
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Here one should impose proper phase-shifts to the fields so that the masses become positive; thus N x /s etc. 
are complex in general. Also the small mixing between /II and /7r must not be neglected, which actually 
does contribute to the final result. 

With these obtained couplings, employing the usual calculation of Feynman diagrams, we can obtain the 
MSSM contribution Aa^ to the muon anomalous magnetic moment :- (g^ - 2)/2 as 
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4.3.2 Formulae in gauge eigenstates 

The above formula ( |4.52| i is accurate and stiff. However, from a theoretical viewpoint formulae in the 
gauge eigenstates are more interesting and convenient even at the expense of accuracy. The formulae can 
be obtained from the Feynman diagrams shown in Fig. 4.3 with the mass insertion method, which results 
in [80] 
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(e) A fljU C7? R (bino)") = - ^-m 2 • M ' ■ • F b (m x , //; m^) , (4.57) 
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where 
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Since dA^ 2 (x)/dx < and dCi(x)/dx > for x > 0, these functions F a and Fb return positive values for any 
mass parameters. 
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Figure 4.3: The same as Fig. 4.2 but expressed in terms of the gauge eigenstates. 



4.3.3 Interpretation 



As is discussed in Sec. 2.2 we need a positive shift of the muon g - 2 to explain the observed value, and 
are happy if the shift is Aa^ — hl0~ 9 . From this viewpoint, here we will briefly interpret the obtained result 
(|433)-(|437]>. 

First let us take the limit where all the relevant mass parameters are common: Mi = M2 = = m^ L = 
m^ R = m? =: Msusy- Then, as the functions return F a — > 1/4 and Ff, — > 1/12, the results are approximated 
to be 
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In this case the chargino diagram (a) gives the dominant contribution, and from the expression we realize 
A^susy ~ (5(100) GeV is required to explain the discrepancy in this degenerated scenario. 

What will happen if one of the smuons is extremely heavy? The answer is simple: corresponding 
neutralino contributions just vanish. In particular, the case where mjj L » mj; R is very interesting. The 
neutralino contribution becomes negative with fi > 0, and moreover, since ~ mji L is expected, the total 
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contribution is given as 

[m^ ~ mj; L » ntfa] Aa^ = ~^ m l ' ' ^ — ■ ^fe (^i./t; »%) ■ (4-61) 

Thus yU must be, different from the usual scenario, negative to explain the discrepancy. 

The other extreme cases are Mi ~ M2 <s yu and Mi ~ M2 » yu. In the former only the diagram (a) 
contributes. Both of the smuons must be light, and /i > should be hold as usual. The latter case has no 
diagrams to contribute: Aa^ - 0. 

The above results are summarized in the next table. 



Table 4.2: Summary of the SUSY explanation of the muon g - 2 in several extreme scenarios. 
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Section 4.4 The Gauge Coupling Unification 



Remember the discussion in Sec. |2.3| physicists desire the gauge coupling unification to realize the grand 
unification, a unified explanation of the three forces; under the renormalization group running they approach 
to each other at a high-energy scale of ~ 10 15 GeV, but the trajectories have a slight mismatch, which is 



summarized in Fig. 2.3 



The unification of the gauge coupling is a virtue of the SUSY theories. In the MSSM the RGEs are 
modified with the presence of the SUSY particles as well as two Higgs doublets; they are summarized in the 



appendix of this chapter, Appendix 4. A 

In particular, for the gauge couplings, (g\,g2, ^3) := ( 
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where C a (G) is the quadratic Casimir invariant for the adjoint representation of the group, and IJj) is the 



Dynkin index of the chiral supermultiplets appearing in the model, 
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For the MSSM these equations are 



(4.63) 
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(4.65) 



These RGEs result in the well-known result displayed in Fig. 4.4 where the unification is much im- 
proved. These figures are drawn with the two-loop level RGEs, and three sets of tan/3 and the SUSY scale 
Msusy are taken for illustration: the top figure is with (Msusy, tan/3) — (ITeV, 10), the middle one is 
(10TeV,3), and for the bottom (100TeV,40). In any choices the unification is much improved from the 
Standard Model case. 
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Figure 4.4: The renormalization group evolution of the gauge couplings g^, g2 and g\ (black lines, from top 
to bottom) under the MSSM, together with the evolution of Y, (a red line), Yt, (a blue line), and Y T (a green 
line). For comparison the evolutions under the Standard Model are drawn with dotted lines. In the three 
figures tan/3 and the SUSY scale M SUSY are differently set: (M SU sy, tan/3) = (1 TeV, 10), (10TeV,3) and 
(lOOTeV, 40) from top to bottom. Yf, and Y T appear only in the bottom figure: they are below the plotted 



range in the other figures. The input values and the schemes of the RGEs are summarized in Appendix 4.B 
together with a note on approximations employed in drawing these figures. 
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Section 4.5 Current Status of the MSSM 

— Higgs, g - 2, and LHC SUSY Searches 

Let us summarize what we discussed in this chapter. 

First we saw the mass of the Higgs boson under the MSSM greatly depends on the stop sector, especially 
on the stop mass and the stop mixing parameter X, = A, - p cot/3. To realize the mass 126 GeV under the 
MSSM, the stop mass is required to be (9(10) TeV without large stop mixing, and even with the maximal- 
mixing of X, ~ V±6»if, it should be ~ 1-2 TeV. This fact prefers heavy SUSY. 

On the other hand, the argument of the little hierarchy, embedded in the Higgs potential, prefers light 
SUSY; in particular, light m^. This is obviously in collision with the argument of the 126 GeV Higgs mass. 



We then, in Sec. 4.3 saw that the muon g - 2 discrepancy is explained with the SUSY contributions 
if neutralinos, charginos, and smuons are of order 100 GeV. Here light SUSY is preferred; note that this 
requirement is not for the colored sector as the previous two, but for the non-colored sector of the MSSM. 

Now let us move on to the discussion on results from the LHC experiments. In this dissertation we 
concentrate on the R -parity conserving SUSY; the other possibility, R -parity violating SUSY, is beyond the 
scope of this dissertation and not discussed. 

The most promising channel in LHC SUSY searches is pair-production of colored superparticles, that is, 
squarks <7 and gluinos ~g~. Especially, the events from pp — > ~gg, ~gcf*^ and Tf-^Tf*^ channels can be detected well 
in the detectors; multiple hard jets plus large missing energy signature, where jets come from cascade decays 
of superparticles, and large missing energy is provided by the LSP. However, we have not detected such 
signatures. Now colored superparticles ~g and q-*' are approximately constrained as rrtg > 900 GeV and nig > 
1400 GeV, respectively, which are relatively heavier than we optimists expected. These constraints were 
obtained by the ATLAS collaboration [81 1, where they analyzed their data corresponding to an integrated 
luminosity of 5.8 fbr 1 collected at the 8 TeV LHC. The CMS collaboration also analyzed their ~ 5 fbr 1 data 
obtained at Ecm - 7 TeV, and obtained similar results I82, l83l . 

Searches focusing on the third generation squarks are also employed. These searches, especially direct 
stop searches, are important for the following two reasons. First, the lighter stop is expected to be lighter 
due to the mixing, and thus be produced more, than the other squarks. In addition, stop mass is crucial for 
the Higgs boson mass and the discussion of the little hierarchy. Such searches are performed with requiring 



leptons and/or fr-jets (cf. Sec. 3.2 1, which are expected from the decay of top quarks. However, we have 
observed no excess, and the stop mass is constrained as > 500 GeV as long as the separation between 
the stop mass and the LSP mass is larger than ~ 100 GeV l84ll85l . 

These two LHC results indicate that colored sector of the MSSM is heavy, which is also supported by 
the 126 GeV Higgs mass. 

Here one possibility arises: colored superparticles are heavy, but non-colored are light to keep the SUSY 
explanation of the muon g - 2 anomaly. The SUSY with such a scenario can be searched with focus- 
ing on pair-production of charginos, neutralinos, and sleptons via electroweak interactions, i.e., pp — > 
^^,^^1,7^/1*, etc. Generally it is difficult to be searched for at the LHC, a hadron collider, because 
hard jets which are useful to distinguish signal events from background ones are not expected; generally 
multiple leptons are instead required in relevant searches. Obtained bounds are ~ 300-500 GeV on masses 
of^J and^|, but these bounds depend to a large extent on the mass spectrum in non-colored sector, and thus 
further searches are required. Here it should be noted that these searches do not need higher energy, but a 
higher luminosity. As the instantaneous luminosity is to be increased and data of (9(100) fbr 1 are expected 
in the 13-14 TeV runs of the LHC, the runs are of great importance for this scenario. 

The following prognostication should also be emphasized: if the three forces are unified under the SU(5)- 
GUTs, gaugino masses M a are expected to be yielded from a common mechanism; then, if charginos and 
neutralinos have masses of order 100 GeV, gluino is expected to have a mass of the same order. For example, 
if they have the same masses M\ = = M3 at the GUT scale, their masses are expected to be 

Mi : M 2 : M 3 * g\ : g\ : g\* 1 : 2 : 7 (4.66) 



at the low-energy scale. Therefore, if masses of charginos and neutralinos are < 1 TeV, the gluino mass 
is expected to be < 3.5 TeV. This possibility, i.e., charginos, neutralinos, and sleptons have (9(100) GeV, 



Chapter4 The MSSM and Its Current Status 



35 



gaugino mass is a few TeV, and the squarks are too heavy to be produced at the LHC, is also expected to be 
focused on at the 14 TeV LHC. 

* * * 

However, from a theoretical viewpoint, this scenario with heavy-colored light-non-colored particles is 
not preferred, because it does not respect the SU(5)-GUTs. The masses of the squarks and the sleptons 
should be close, because they are considered to have a common origin, i.e. to originate from the same 
multiplets, under the SU(5)-GUTs. Actually, this SU(5)-GUTs hypothesis is fully respected in the simplest 
models for the SUSY-breaking, such as the CMSSM (constrained MSSM) scenario, which is sometimes 
called the mSUGRA (minimal supergravity), and the GMSB (gauge-mediated SUSY-breaking) scenario. 

Unfortunately, experimental results do not support this anticipation. Under the assumption that the 
masses of squarks and sleptons are at the same order, it is not so easy to realize the 126 GeV Higgs mass with 
explaining the muon g - 2 anomaly. Especially, it is known that the CMSSM scenario does not realize these 
two features of the SUSY [ 14, 86|, and the GMSB neither. Nature seems to force us to leave this dream. 

The GMSB scenario is, however, too beautiful to be abandoned; the SUSY CP- and flavor problems do 
not arise in this framework. Therefore, we consider an extension of the MSSM: the V-MSSM, in which a 
vector-like pair of the decuplets of the SU(5), 10 + 10, is appended as extra matters to the MSSM. In the 
model, which respects the SU(5)-GUTs, the extra matters contribute to raise the Higgs boson mass, and 
therefore the squarks do not need to be so heavy as in the MSSM case. Since the muon g - 2 is scarcely 
affected by the extra matters, both features, the 126 GeV Higgs boson and the SUSY explanation of the 
muon g — 2, can be simultaneously realized within the simplest SUSY-breaking frameworks, the GMSB and 
the CMSSM. 

In the next chapter, we will discuss the V-MSSM. Especially, we will examine the V-GMSB scenario, 
i.e., the V-MSSM with the GMSB framework, and to see that the 126 GeV mass and the explanation of the 
muon g - 2 problem can be simultaneously realized, even under the current LHC constraints. This is the 
main topic of this dissertation. 
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Appendix 4.A Renormalization Group Equations for the MSSM 

In Chapter|5] we will discuss an extension of the MSSM, called the V-MSSM, and introduce the renormal- 
ization group equations (RGEs) for the model. Here, paying the cost of verbosity for completeness, the 
RGEs for the MSSM parameters are displayed up to two-loop level. 

The y6-functions are calculated with Susyno 1.1 1 87 1 and checked by comparing with Refs. 



4.A.1 Restriction and notation 

Here we simply employ the following assumptions, that is, we ignore CP- and flavor violations thoroughly. 

• The R -parity is conserved. 

• The scalar soft mass terms m x are diagonal, 

• For the A-terms ax and the Yukawa coupling Yx, all the components but the (3, 3) are neglected. 

• The gaugino masses M a are real. 
The y6-function is defined as 



1 



dX(Q) 
dlogg 16tt 2 



[X] 



1 



(16tt 2 )- 



jS (2) [X] , 



(4.67) 



where X are one of the MSSM parameters. Their definitions are found in Eqs. ( |4.4| i and ( |4.5| l. and Q is the 
renormalization scale. The DR scheme (65) is chosen as the renormalization scheme. 
The following variables are used in the expressions of the /3-functions. 



X, 
X h 

X T 

a (l,b,T) 



= 2a] + 2Y 2 {m 2 Hu + (m^) 33 + (m|) 33 ) , 
= 2a 2 b + 2Y 2 h (m 2 Hd + (m 2 Q ) 33 + (m\ ) 33 ) , 
= 2a 2 + 2Y 2 (m 2 Hd + (m 2 ) 33 + (m|) 33 ) , 

= Y( t ,b,T)Cl(t,b,T), 

3 

= m 2 Hu - m 2 Hd + 2^ [(m 2 Q )a - 2{m 2 ) ii + {m\)n - (m£) a + . 



S 

S i2) := -Y 2 (3;4 u + (m 2 Q ) 33 - 4(m 2 & ) 33 ) + Y 2 b (3m 2 Hi - (m 2 ) 33 - 2(m 2 ) 33 ) + Y 2 T (m 2 Hs + (m£) 33 - 2(m|) 33 ) 



+ [To g * + 3 2 g2 



2 1 
m H u ~ m Hi 



1 2 , ^ j2 8 2 



30 yi 2 

3 



+ -</~2 + •; iX""^'" 



0"i := 



y 0? + y gi) Yjin&a + (y 9i + f »f) + §0? 

' i=l * ' 1=1 (=1 

3 

3m 2 Ha + 3m 2 Hi + ^ ((to 2 ,),,- + 8(m 2 ? ), 7 + 2(m\ ) u + 3(m 2 ) l7 + 6(m|) l7 ) 
(=i 

3 



CT 2 := 2 

°"3 := 3 X (2(»ig), 7 + (m^),-,- + (m 2 ,),,) . 



4 The chosen scheme of Ref. 1671 was originally noted as the DR scheme, but actually was the DR' scheme, as is mentioned in note 
added in the paper or, e.g., Ref. 1 89 1. 
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4.A.2 One-loop level ^-functions 



^ [9a 

[M a 

ft" [Or 

fF>\b 

p m K 
p (X) K 

[{m 2 Q )u 
P m [(m 2 D ) u 

P W [(ml). 



"a Va 



where (b^,B^,B^) 
Yl-^ g 2 -3 g l-f g j + 6Y?y t 



33 



-,1,-3 



6Y, 



2 7 „2 o 2 16 , 2 . v 2 



15 



gt-3ti-^-Sf 3 + Yt + Y;\Y b 



3Y£--gi-3g£+4Yt\Y T 

/2 , v 2 13 „2 o 2 1" 2 



26 



32 



1«7 + ^ - - 3gi - —gij a, + \2a h + -g\M x + 6£M 2 + — g\M 3 1 Y, 

Y 2 + 18F 2 + F 2 - ^ 2 - 3^2 - y 03) Ob 

(14 32 \ 

2a, + 2a T + y^!? 2 ^ + 6gjM 2 + —glM^Yb 

3Y 2 + UY 2 T - 9 -g 2 - 3^| a T + + y^ 2 M, + 6a 2 M 2 ) F T 

3^-^-3^ + 37 2 + y 2 L 

6a ( + 6a b + 2a T + ^a 2 M, + 6glM 2 Jfi + \3Y? + 3Y\ + Y 2 - ^g 2 - 3g^ b 
~g 2 M 2 -6g 2 2 Ml+^g 2 S +3X t 
--g\M\ - 6g 2 2 M 2 2 - ^g 2 S + 3X b + X T 



= -Y 5 9\M\ - 6g 2 2 M 2 2 - ^g\M\ + l -g\S + ((x, + X,,))^ , 



= -g^ 2 -^M 2 -^ + ((2X t )) for! . = . 



32 



-g\M\ - -g\M\ + -g 2 S + ((2**))^ 
^ 2 M 2 -%2^2 2 -^^+((x T )) 



for i = 3 



(4.68) 
(4.69) 
(4.70) 
(4.71) 
(4.72) 

(4.73) 

(4.74) 
(4.75) 
(4.76) 
(4.77) 
(4.78) 
(4.79) 
(4.80) 
(4.81) 
(4.82) 

(4.83) 
(4.84) 



4.A.3 Two-loop level ^-functions 



0- 2) [gi] = 

£ (2) fe] = 



26 



■Yf 



14 



18 



199 



27 



n-^Y: + -gi + -g^-g 2 3 g] 



-6Y 2 - 6Y 2 h - 2Y 2 + 9 -g\ + 25g 2 2 + 24^ 2 J g\ 



11 

y 



-AYf-Wi+'-gl+^l + UgWgl 



(4.85) 
(4.86) 
(4.87) 
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j8 (2) [M^ = 

j8 (2) [M 2 ] = 

/3 (2) [M 3 ] = 
/? <2) [Y t ] = 

P i2) [Y„] = 

/3 (2) [Y T ] = 

/? (2) [fl,] = 



j8 (2) [a*] 



52 (a, - if Mi) + 28 (a 6 - Y 2 h M x ) + 36 (a T - Y 2 M X ) 
796 



+ — g\M, + 54g z 2 (M, + M 2 ) + 176^ (Mi + M 3 ) 



0t 



12 (a, - if M 2 ) + 12 - YlM 2 ) + 4 (a T - F 2 M 2 ) 
18 



+ — (Mi + Mi) + \00g z 2 M 2 + 48g z 3 (M 2 + M 3 ) 
22 



02 



8 (a, - F 2 M 3 ) + 8 (a* - if M 3 ) + — g\ (M x + M 3 ) + 18g 2 (M 2 + M 3 ) + 56 9 2 M 3 
f J? + ^n 2 ) 0i + + 16lf ^ - 22Y? - 5Y 4 b - 5Y?Y 2 b - Y 2 b Y 2 



27 43 4 I64 99 136 ^ 1 9 9 

+ 4gQ-»i + y02 - y0 3 + 0102 + "^-0103 + 8 0203 



(4.88) 



(4.89) 



0? (4.90) 



(4.91) 



Yt 



\y 2 - 2 -Y 2 h + 6 -Y 2 }g\ + 6Y 2 h g 2 2 + \6Y 2 h g\ - 5F 4 - 22F 4 - 5Y 2 Y 2 b - 3Y 2 h Y 2 - 3Y 4 T 



287 4 15 4 16 4 22 ^22 r>22 

+ + y02 - y03 + 0102 + g0103 + 8 0203 



Y b 



2 -Y 2 h + 6 -Y 2 \g\ + 6Y 2 T g 2 2 + \6glY 2 h 



- 9Yt - 10Y 4 - 3Y 2 Y\ - 9Y 2 b Y 2 + f , 4 + ^ + \g\g\ 
^g 4 (a, - 4F ( M,) + y g 4 (a, - 4F ( M 2 ) + ^ g 4 (4F,M 3 - a,) 



(4.92) 



(4.93) 



+ ^ 2 (a, - 27,M, - 2F,M 2 ) + (fit ~ 2W " 2K ( Af 3 ) 



136 



.2 „2 . 



+ 89^2 (a, - 2Y,M 2 - 2Y,Mi) + g\ \^a,Y z + -Y,a„ + ^a,Y, - -Y,Y b M\ - —Y, 
+ g\ (-12M 2 F, 3 + 18a,F,) + g\ (48a,F, - 32M 3 lf ) 

- 1 10a ( lf - a,Y 2 Y 2 - 2Y t Y 2 a T - 2Y,a h Y 2 - 5a,Y 4 h - \5Z t Y t Y 2 b - 2QYmY 2 b - \QY^a b 

(4.94) 

287 15 16 

— g\ (a b - 4Y b M 1 ) + —g\ (a b - 4Y h M 2 ) - —g\ (a b - 4Y h M 3 ) 



,2 / ^ n v2 4 . 



18, 



12. 



; 3 mi) 



+ g\g\ {at - 2M, Y b - 2M 2 Y h ) + "-g\g\ (a b - 2Y b M l - 2Y b M 3 ) 



2„2 , 



+ g\g\ (%a b - \6M 2 Y h - 16M 3 Y b ) 

'4 



+ g\ 



12 _ 4 



6_ 



-Y 2 a b - % -M x Y 2 Y b + ^a,Y b + ^a b Y 2 - ^-M^Y 2 + ^Y b a T - ^M,F 3 + -a b Y h 



g 2 (lSa h Y b - \2M 2 Yl) + g 2 (48Z b Y b - 32M 3 F 3 ) - 1 10a b Y 3 h 

5Y 4 a b - 20a t Y 2 Y b - \5Y 2 a b Y b - 10a f if - 3a b Y 4 - 9a b Y h Y 2 - 6F 3 a T - 12Y h a T Y 2 



(4.95) 



Chapter4 The MSSM and Its Current Status 



39 



/? (2) [flr] = - AM\ Yj) + ^-g 4 2 (fir - 4M 2 Y T ) + \g\g\ (a T - 2Y T M\ - 2Y T M 2 ) 



P (2) [b] 



4_ 



12 



+ 9x\— 5 n^ + l-M^Yr - -a b Y T - -jM^ 



18_ 



-a T Y T 



+ g\ (l8a T F T - 12M 2 F T 3 ) + 93 (l6y£a T - 32M 3 F fc 2 F T + 32a fe y T ) - 3lf l^a T 

- 67i,y 2 y T - 6y, 2 fli,y T - 9y^ T - 367i fc y 2 y T - 27y 2 7i T y T - \u h Y\ - 50a T y T 3 



y?-^y 2 + ^y 2 ) 2 + i6(y 2 + y 2 ) 2 



- 9Yf - 9Yt - 3Y\ - 6Y 2 Y 2 h + + ^ + ^ 2 

I (Si - y 2 M!) g 2 - 1 (a, - y 2 Mj) + y (or - y 2 M x ) 9 2 

+ 32 (a, - M 3 Y? + a b - M 3 Y 2 ) g\ - 12 (3a,y 2 + 3a h Y 2 b + a T Y 2 + a,Y 2 + Y 2 a b ) 

414 18 

- ^9\Mx - 30g 4 2 M 2 - —g\g\ (M, + M 2 ) 



y 2 -^y, 2 + ^y 2 )g 2 + i6(y 2 + y 2 ) 9 2 



_ 97 4 _ 9y 4 _ 3y 4 _ 6Y 2 Y 2 + |^4 + + |^ 

621 18 6 3 

— + 33^M 2 + — g\g\ (m 2 + M 2 2 + M X M 2 ) + -g\S (2) + -g\(T X + 3g\a 2 

+ \o\ (*» - *M x a, + \M\Y 2 ) + 16<? 2 (x, - 4M 3 a t + \M 2 Y 2 ) 

- \M 2 X t - 36S 2 - 3Y 2 b X t - 12a,a b - 3Y 2 X b 

621 18 6 3 

— g\M\ + 33g A 2 M 2 2 + — g\g\ (m 2 + M 2 2 + M X M 2 ) - -g\S (2) + -g\(T X + 3g\a 2 

+ X -g\ {$M x a b - %>M 2 Y 2 h - 2X b - 24M \a T + 24M 2 Y 2 + 6X T ) 
+ 16g 2 3 (X b - 4M 3 a b + \M\Y 2 ) 

- 127J,a 6 - 3YfX b - \%Y 2 h X b - 3y^X ; - 36% - 12a 2 - 6Y?X T 

199 128 2 

g\M\ + 33g\M\ - —g\M\ + -g\g\ (m 2 + M X M 2 + M 2 2 ) 



75 



32 
45 



g\g\ {m\ + M 2 + Mi M3) + 32g\g\ (m\ + M\ + M 2 M 3 ) 
+ \9\S (2) + ^ 2 o"i + 3g\a- 2 + y^<r 3 
+ i^gi (4X, - \6M{a, + \6M\Y 2 + 2X b - t>M{d b + 8M 2 y 2 ) 

- lOYfX, - 20a 2 - \QY 2 h X b - 20a 2 b - Y 2 h X T - Y 2 X b - 4a b a T Jj 



(4.96) 



(4.97) 



(4.98) 



(4.99) 



(4.100) 



(4.101) 



for i = 3 



40 



Dissertation / Sho Iwamoto 



[imlh] = - fg\Ml + 5 -^ g \gl (m\ + M\ + M,M 3 ) - ^ + j|rf<r, + y ^3 

J?? + ^Mjfl, - \m\Y^ + g\ (6X t - 24M 2 a, + 2\M\Y 2 ) 

- 16Y 2 X, - 32a 2 - 2Y 2 b X t - Sa t a b - 2Y?X b \) 

''fori = 3 



(4.102) 



+ |s? - ^MiOj + ^M 2 F 2 J + 9 2 ( 6Xi - 24M 2 a b + 24M 2 F 2 ) 

- 8a b a T - 16Y\X b - 32a 2 - 8a,a b - 2Y 2 X h - 2Y 2 h X, - 2Y 2 h X T - 2Y 2 T X h \ 

//for i = 3 

(4.103) 

P i2) [(«!)»] = !?X + 33^m| + y a 2 9 2 ^ M 2 + M 2 + MiMz ) _ 6^(2) + | gr 2^ i + 

+ |<7? |^X T - yM,a T + yM 2 ^ 2 ) - 123^ - 3Yf,X T - 3Y 2 X b - 12a 2 - 6Y 2 X T Jj 



for i = 3 

(4.104) 



P (2) [(»£>«] = + y> 2 S< 2 > + fgfr 

+ l^gl (- 6 -X T + y M x a T - y M 2 F 2 J + g\ (6X T - 24M 2 a T + 24M 2 F 2 ) (4 _ 105) 

- 24^ - 6F^X T - 6Y 2 X h + -16a 1 . - 8F 2 X T \\ 

'/fori = 3 
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Appendix 4.B Note on Figures Showing Gauge Coupling Evolution 

Here are documented several remarks on the figures displaying the renormalization group running of cou- 
pling constants, i.e., Fig.|2.3|and Fig|4.4| 



Evaluation for the Standard Model (Fig. \2.3) The gauge coupling running under the Standard Model is 
employed in the MS scheme [39]. The RGEs for the gauge couplings g a are evaluated at the two-loop level. 
The Yukawa couplings appearing in the RGEs are run under the one-loop level RGEs, which is sufficient to 
obtain the gauge couplings at the two-loop level accuracy. 

Input for the coupling constants are evaluated at the electroweak scale, i.e. at the Z-boson mass m z . The 
values are 119011 

^(mz)*® =0.1184, sin 2 <9 w (m z )^ = 0.2312, a EW (m z f^ = 127.9, (4.106) 

which results in g\ = 0.462, gi — 0.652 and g^ = 1.22. Note that gi = V5/3gy. 

For precise calculation the Yukawa couplings should be also evaluated at mz with the MS scheme. 



However, in drawing Fig. 2.3 the constants are evaluated as 



m,(m t ) MS 160 m b (m b ) MS 4.18 m? ole 1.77682 

i t % % 1 1. % % /_ % % . (4.1U/) 

v 174 v 174 v 174 



Therefore, precisely speaking, the running of Fig. 2.3 should be understood as a rough estimate; but this 



simplification would be compensated with the large uncertainty of the MS top mass and the ineffectiveness 
of Y b and Y T in the evaluation. 



Evaluation for the MSSM (Fig. 4.4 1 As we are not interested in the precise values, several approximations 



for simplicity are employed in numerical evaluation for Fig. 4.4 First the procedure is summarized, followed 
by the discussion on the approximations. 

The ratio tan/? := u u /i>d and a SUSY scale Msusy are set as input values. The running below Msusy is 
employed with the Standard Model MS RGEs, and that above M SU sy is with the MSSM DR RGEs. The 
exchange at the SUSY scale is employed as 

gf ssu = gl M , F™ = J/ SM /sillA Tf SSM = j/fV cos/?, T T MSSM = y* M / cos p. (4.108) 

Input values at the electroweak scale are taken as the same as the Standard Model case discussed above. 

One should first note that the renormalization schemes are different above and below the SUSY scale. 
Nevertheless, the RGEs in the two schemes are the same at the one-loop level expression, and thus it is 
expected that this approximation does not cause serious problem. 

More important simplification is that no threshold corrections at Msusy are introduced. Precisely speak- 
ing, the matters decouple gradually during the running down of the mass scale Q, and thus the RGEs, or the 
coefficients in the RGEs, depend on the scale Q. However, in usual calculations the RGEs are treated as if 
they are independent of Q, and instead so-called "threshold corrections" are introduced at the SUSY scale 
to compensate the effect from the matter decoupling. Therefore, ideally, we have to introduce the threshold 
corrections into the conversion of Eq. ( |4.108| i; but as they depend on the SUSY mass spectrum, they are 



simply not introduced in the numerical evaluation for Fig. 4.4 



Consequently, the running of Fig. 4.4 should be understood as a rough illustration of the gauge coupling 
unification. 



Nevertheless, the numerical evaluations in Chapter[5] where we utilize SOFTSUSY 3 . 3 ll9TI to calculate 
the mass spectrum, the above subtle effects are introduced; i.e. the approximations discussed here are not 
performed. 
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Chapter 5 

The MSSM with Vector-like Matters 



Now we are ready to discuss the "V-MSSM," an extension of the MSSM with a vector-like pair of super- 
multiplets. This model ll92l is characterized as 

V-MSSM = MSSM + (10 + To) , (5.1) 

where 10 and 10 are the decuplet and the anti-decuplet of the SU(5) gauge group. 

In this model the Higgs mass can be raised by the extra vector-like quarks |92|. There very large squark 
masses of ~ lOTeV are not required to realize the 126 GeV Higgs boson even with a small mixing parameter 
X t as we will see in Sec. |4. 2.4] This feature further allows us to utilize the SUSY as a solution to the muon 
g-2 anomaly even in the simplest SUSY-breaking frameworks. One should also note that this model respects 
the underlying SU(5) symmetry. 

As the gauge-mediated SUSY-breaking (GMSB) scenario is very promising for its freedom from the 
SUSY CP- and flavor problems, the V-MSSM with the GMSB is particularly interesting. In the GMSB 
scenario, the slepton soft masses have the same origin as those of the squarks, and the A-terms are much 
smaller than the soft mass terms. Therefore, under the MSSM framework, the GMSB cannot realize the 
126 GeV Higgs mass with holding the SUSY explanation of the muon g-2 anomaly. Here, the "V" resurrects 
the GMSB scenario. As the extra matters raise the Higgs mass, the mass of 126 GeV is realized without 
exploiting heavy squark masses of order 1-10 TeV. It allows the lighter masses of the sleptons and gauginos 
even under the GMSB framework, and with those masses the muon g-2 discrepancy can be explained. 

In chis chapter we will examine the V-MSSM, especially focusing on the combination with the simple 
GMSB scenario, which we call "V-GMSB model." 

Historically, the V-MSSM was introduced in Ref. [92], where the increase of the Higgs boson mass was 
unveiled. In Ref. [93 1 the model was examined; especially discussed were the infrared fixed point behaviour 



of the couplings and decays of the extra fermions, which will in this dissertation be explained in Sec. 5.1 
and Sec. |5.7| respectively. 

Then author's works follow.* 1 In Ref. ifTBI the V-MSSM was firstly combined with the muon g-2 
problem. In Ref. l94l the LHC phenomenology is discussed, and Ref. |16| was devoted to the vacuum 



stability condition, which we will see in Sec. 5.4 Ref. [95 1 was on constraints from SUSY searches at the 



LHC, which corresponds to Sec. 5.6 of this dissertation 



In summary, notable features of this model are: 

• It respects the underlying SU(5). 

• It goes perturbatively up to the GUT scale. 

• It allows us to solve the (g - 2) M anomaly even with realizing m/, ^ 126 GeV. 

After defining the model, we will see these advantages, and then discuss current constraints on this model. 



'This chapter is based on the works by Author, completed in collaboration with Dr. M. Endo, Prof. K. Hamaguchi, Mr. K. Ishikawa 
and Dr. N. Yokozaki 1 15 16 94 95); a conference article by Author was published in proceedings [96|. 
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Table 5.1: The field content of the V-MSSM. The gauge indices are omitted, and subscripts i denote indices 
for the three generations running from one to three. The gauge group is the same as the MSSM, and thus the 
gauge fields are the same. 



Matter fields (chiral multiplets) 


field 


SU(3) color SU(2) weak U(l) y 


/^-parity 


Q 




3 2 1/6 









3 1 -2/3 




E 




1 1 1 




D 




3 1 1/3 




L 




1 2 -1/2 




H u 


1 2 1/2 


+ 


H d 


1 2 -1/2 


+ 


Q' 


3 2 1/6 




U' 


3 1 -2/3 




E' 


1 1 1 




Q' 


3 2 -1/6 




U' 


3 1 2/3 




E' 


1 1-1 





Gauge fields (vector multiplets) 


gauge group 


field 


SU(3) co l or SU(2) weak U(l)y 


SU(3) co i or 


G 


8 1 


SU(2) weak 


W 


1 3 


U(l)y 


B 


1 1 
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Section 5.1 The V-MSSM Model 
5.1.1 Definition 

The V-MSSM is an extension of the MSSM, where a vector-like pair of the SU(5) decuplets is introduced. 
That is, it is defined with the same symmetry as the MSSM, 

symmetry: [SU(3) strong x SU(2) weak x U(l)y] x Z\ , (5.2) 

and the extended field content, 

fields: MSSM + 10 + 10 = MSSM + (q, U',E') + (q' , U',E') . (5.3) 

The field content is summarized in Table 15. ll 

These ingredients yield the following generic superpotential, 

Wmssm = fiH u H d - (Y u ) u H u QiUj + (Yj^HtQiDj + (Y e ) u H^Ej, (5.4) 
W ex , ra = -Y'H U Q'U' + Y"H d Q'U' + MqQ'Q' + M^U'U' + M E ,E'E', (5.5) 
^mixing = - fe); H a QiU' - «),. H u Q'Ui + H d Q'Dj + {e e \ H d Lj£', (5.6) 

and SUSY-breaking terms 

--Csoft = --C™l SM + m 2 Q , \Qf + m\, |ef + m%, |[7'| 2 + <, \uf + m% \E'\ 2 + m% \E'\ 2 

+ {-a'H u Q'U' + a"H d Q'U' + H.c.) + (b Q ,Q'Q' + b^U'U' + b E >E'E' + H.c.) (5.7) 
+ [- OU H u QjU' - (A' u \ H u Q'Ui + (A d )i H d Q'D t + {A e \ H d L,E' + H.c] . 

Note that extra yu-terms, such as [i,Q,Q ', can be absorbed into the e-terms, and similar absorption can be 
operated for the SUSY-breaking terms. 

However, this general Lagrangian has several problems. The principal one is on the mixing terms e's be- 
tween the Standard Model fermions and the extra fermions. These mixings must be small enough, especially 
for i = 1,2, not to evade current experimental bounds on flavor violating processes. However, on the other 
hand, absence of these mixings causes a problem that some of the extra fermions, especially the lightest 
vector-like quark (t\ , which will be introduced soon), become stable. This is obviously disfavored because 
we know our Universe has no heavy stable colored particles. Quantitative evaluation of these mixing param- 
eters is interesting, and has great importance in search for the extra fermions at the LHC, which is discussed 



in Sec. 5.7 because the decay branch of the fermions is determined by these mixings. Nevertheless, the 
quantitative discussions are left as future works. In this dissertation we simply restrict the model to have no 
mixings between the vector-like matters and the Standard Model fermions in the first and second genera- 
tions, and to have small mixings between the vector-like matters and the third generation.* 2 Furthermore the 
SUSY-breaking mixing parameters A's are all ignored. These terms just govern the mixing between MSSM 
sfermions and the extra sfermions, and are thus not important in phenomenological discussions. 

Later we will impose one more condition that Y" = 0; this is because the term Y"H d Q' U' tends to reduce 



the Higgs boson mass if tan/? is as large as ~ <9(10). This mechanism is discussed in Sec. 5.1.3 For the 
SUSY-breaking sector, similarly, a" = is imposed. 

Therefore, in the following discussions, we use the superpotential and the SUSY-breaking terms of 

(5.8) 



Wvmssm = Wmssm - Y'H U Q'U' + Mq.Q'Q' + M W U'U' + M E ,E'E' 
- e u H u Q 3 U' - e' u H n Q'Ui + e d H d Q'D 3 + e e // d L 3 £\ 

_ X VMSSM = _ X MSSM + ^ |g,|2 + ^ + ^ (f? ,|2 + ^ ^2 + ^ ^2 + ^ ^2 

+ [-a'H a Q'LJ' + b Q ,Q'Q' + b^U'U' + b E >E'E' + H.c.) , 

2 Discussions on the mixings with the first and second generation Standard Model quarks are found in Ref. |:97) etc. 



(5.9) 
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where e's are assumed to be tiny; more precisely, they are small enough to avoid disfavored flavor mix- 
ings, but not so feeble as the longevity of the extra matters does not cause phenomenological, especially 
cosmological, problems. 

We also define the following variables for later convenience: 



A' := a'/Y', 



a"IY", B & :=b a IM Q ,, B L 



bu'/Mo-. 



B E , := b E 'IM E , 



(5.10) 



Actually there remains one subtlety: CP-violating phases in the parameters. First, in the supersymmetric 
sector, one phase cannot be removed from (Y' , Y" ,Mq>,Mu>) because we already fixed the phases of H u and 
H& to set u u and v A positive. This does not matter under the simplification Y" — 0, but anyway we just assume 
all of (1", Y" ,Mq,Mu') are positive for simplicity. In addition, the SUSY-breaking terms yield many phases 
as is the case of the MSSM, but we simply neglect all the complex phases in these terms. 



5.1.2 Extra particles and their masses 



In this model we have eight extra scalar particles and four extra fermions. Let us discuss the masses of the 
particles at first. 

The starting point of our discussion is Eqs. ( |5.8| l-( |5T9] >, but we restore Y" for completeness. After the 
electroweak symmetry breaking, where and acquire vacuum expectation values of u u = usin/3 and 
i>d = vcos/3, the fermionic mass terms in the Lagrangian are expressed as, with mixing terms, 



(M Q , Y" Vi \ 



Y'v u M v , e u v a 
v e,X m, 



Q'u 
v 



-£=>(& U' t R ) 



(5.11) 



where t E R, b E ,R and t E r are the Standard Model top, bottom, and tau, and their masses are denoted by m t , 
nib and m T , respectively. Similarly, the scalar mass terms are obtained as 



-£d(g';c/" Q' u u')m u 



(OS 

U' 

Q«* 

U'*) 



• (e;; & d ) M d (g) + (et e>) m c (5.12) 



where the mass matrices are 



Y'v u M* Q , + Y"*v d M Q , 
b Q> 

a'v u — fi'Y'vd 



-b a M d 



b* , 

S) < 



Y"v u Mu, + Y"v A M* 

a"vd — /j*Y"v u 
b v . 



M, = 



Q 

a"*Vd - fiY"*v u 
M™ 
Y'v u M* Q , + Y"* V dM v , 



b* 

(2,2) 



a'*v u - fiY'*Vd ^ 
Y'*v u M Q , + Y"v A M\j, 

Mr ] 
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and the diagonal components are 



\M Q >\ 2 + \Y'v u \ 2 + m 2 Q , + (~ - ^ sin 2 6> w j m 2 z cos 2/3, 
M {2 ' 2) = \M V > | 2 + \Y"v A \ 2 + ml, +\~ sin 2 <9 W j m| cos 2/3, 
M[, 3 ' 3) = |M G ,| 2 + \Y"v A \ 2 + m|, + |-i + | sin 2 W J m| cos 2y8, 
7V(( 4 ' 4 ) = IMf/-! 2 + |i"f u | 2 +«4, +^sin 2 6» w j;74cos2yS, 



|M e < | 2 + m 2 Q , + (- ~ + i sin 2 6» w ) w| cos 2/j. 
|M fi - 1 2 + m 2 , + I ^ - i sin 2 # w ) m| cos 2/j. 



M^ 1 ' = |M £ -| 2 + w|, + (sin 2 6> w )m|cos2/?, 
M (2 ' 2) = \M E , | 2 + m\, - (sin 2 W ) m\ cos 2/?. 



Note that we use the same symbols for both the scalars and the fermions, but it should not be confusing. 
Finally, it is obvious from the above expressions that we have following particles as the extra matters: 



scalars: (^ l ,7 1 7^,7 A ,b\,b' 1 ,'T' l ,~T'^, 



fermions : \t[ , f 2 , b', t / J, 



(5.13) 



with the definition that t' x < t' 2 < t' 3 < t' 4 etc. in terms of the masses. 



5.1.3 Higgs mass increase 

Let us move on to the most important feature, the Higgs mass in the V-MSSM. We have discussed much 



about this topic for the MSSM in Sec. 4.2.3 and the discussion given here is completely parallel to that. 

As the effective potential has the same form as Eq. ( |4.32| i, the increase of the Higgs mass specific to the 
V-MSSM at the one-loop level is given as, in the decoupling limit, 



Am 



2 |l-loop 



IVMSSM 



sin z /? f d 2 

fou 2 



v a dv a 



cos 2 /? / d 2 

to 



i d \ . o 2 

— - — + sin/icos/j- — — 
va dv A j ov u dv d 



(AV« SSM ), (5.14) 



where 



av: 



<i) 

VMSSM 



16tt 2 



F(x) 



x 2 ( x 3 



(5.15) 



Here one should note that the other extra particles, b' n b', t', have no contribution to the Higgs mass in the 
absence of Y". Therefore, with the masses we calculate in the last section, the increase of the Higgs boson 
mass can be easily calculated numerically. 



Let us assume M Q > - Mu>(=: M F ), m 2 , = m% = m^, - m^X-'- M 2 - M 2 ), and b Q > = b v > = 0, in 
order to obtain an analytic expression of the increase. With a straightforward computation, one can obtain 
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the following result: 

"Ivmssm An 2 



log 



Ml 1 



5- 



M 



2 \ 



(r sinj0) 3 (Y" cos/3) 



-(F'sin/3) 2 (y" cos/3) 2 



- ~{Y' sin y8)(y" cos/?) 3 



M 2 ^ 



1 - 



■ + — 

M 2 



M 2 J 



1 - 



Ml 



3M 2 S 12M 2 ) 



M 2 ^ 



Ml)\ Ml) 



M 2 < ~' 2 



Mj {Mj 2M\ t 



1 - 



mV\ m 2 (»' + h") 2 



M 2 J 



3M* 



2- 



M 



2 \ 



+ (y" COS/3) 



M 2 1 
l0g ^"6 



2 
S 

5 T 



1 - 



M 



M 2 j 



Mj 
M? 



v M 2 2M 2 J 



2 \ 



1 - 



Ml) E" 2 
M 2 



Ml) 



1 - 



Mi 



3M 2 12M 2 J 



(5.16) 



with 



S' := A' - /j cot/3, 



A" - /itan/3. 



(5.17) 



Especially, if we set Y" = as is mentioned above and employed in the following phenomenological dis- 
cussions, the increase becomes 



Am: 



,i-ioop 3 V 2 Y' 4 sin/3 4 



Ivmssm 



4tt 2 



M 2 1 
l ° g Mj-6 



5- 



M 2 / 
M 2 , 



1 - 



M 2 F ) E' 2 



Mj) 



1 - 



M 2 
3M 2 



12M 2 J 



(5.18) 



This is similar to the MSSM top-stop contribution ( |4.40| >; the difference simply comes from the fact that the 
extra top-like quarks are vector-like. Therefore, the extra contribution is considered to be comparable to that 
from the MSSM, and thus the Higgs mass is to be increased considerably. 

It is interesting that, similarly to the MSSM case, a large ratio of Ms /Mr yields a larger increase in 
the Higgs mass; the Higgs mass does not raised under Ms = Mf. Therefore, since the parameter Ms is 
characterized by Msusy, a lighter M F gives a large contribution to the Higgs mass for a fixed Msusy- We 
will see this feature in Sec. 



5.1.4 Renormalization group flow 



The unification of the gauge coupling constants is a virtue of the SUSY theories, which are gy = V3/5{7i f° r 
U(l)r, gi for SU(2) wea k and g s = g^ for SU(3) co i or , as is discussed in Sec. |4.4| Here the vector-like matters 
modify the renormalization group equations (RGEs). When we extend the MSSM with n 5 copies of (5 + 5) 
and «io copies of (10 + 10) supermultiplets, the one-loop level RGEs are modified to be 



d#3 


g\ 


ding 


16n 2 


dg2 


g\ 


ding 


16tt 2 


dg\ 


g\ 


ding 


16tt 2 



[-3 + (3nio + "s)], 
[l + (3nio + "5)], 



y^ + ( 3«io+«5)]. 



(5.19) 
(5.20) 
(5.21) 



In the V-MSSM, where (« 10 ,n 5 ) = (1,0), it is interesting that g^ does not run at the one-loop level. 
Including the two-loop level RGEs, we obtain the renormalization group evolution as Fig. 5.1 where the 
V-MSSM parameters are set as: (F, Y", tan/3) = (1.0,0,20) at Q = M SUSY = 1 TeV. In the figure the 
evolutions of Y, and Y' are drawn, but here one should be careful: the evolution of Y' significantly depends 
on its value at the SUSY scale Msusy, which we call y suSY - It is illustrated in Fig. 5.2 In the left figure 
the renormalization group flow for Y' is drawn; we can see that a wide range of Y' G{]T , the value of Y' at the 
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Q [GeV] 



Figure 5.1: The renormalization group evolutions of the coupling constants in the V-MSSM. Here Y' is set 
to be Y' = 1 at the SUSY scale M SUSY = 1 TeV. Other parameters are: Y" = 0, tan/3 = 20. The drawing 
procedure is similar to that summarized in the latter half of Appendix 4.B The RGEs are summarized in 
Appendix 5. A The evolutions of g3, g2, and g\ are drawn with black lines (from top to bottom); that for Y, 
is with a red line, and for Y' is with a turquoise blue line. One should be careful that, as we will observe 
in Fig. 5.2 the running of Y' is sensitive to its value at the SUSY scale, Ts USY , for its infrared fixed point 
behaviour. Moreover, the value of Y' SI]SY affects the evolutions of the other Yukawa coupling constants; their 
flows also significantly depend on Tg USy . However, as the Yukawa couplings appear in the RGEs of the 
gauge coupling constants only at the two-loop level, the renormalization flows of the gauge couplings are 
less affected by Ts USY - 



GUT scale Mquj, converges at the low-energy scale to Y' SVSY ~ 1. This is called "quasi infrared fixed point 
behaviour," and the fixed point is Y' ~ 1 ||93l . In the right figure Y' SUSY is drawn as a function of Y' Gl]T . 

The above discussion assumes Y" = 0, but actually, in the presence of Y" , it behaves in similar manner; 



moreover, Y" also has a quasi infrared fixed point. This is shown in Fig. 5.3 Tg USY in the left figure and 
^s'usy m t ^ le right fig ure are plotted as a function of (Y^ m , Y^ m ). We can see that the both have quasi 
infrared fixed points, and respective behaviour is almost independent of the value of the other. 

It is also known lfT31 l93l that A,, A' and A" among the scalar trilinear couplings perform the infrared 
fixed point behaviour. To illustrate this feature, the CMSSM scenario is utilized, which has the following 



five parameters: (otq, M\ /2, tan /3, Aq, sgn fi). Fig. 5.4 shows the infrared fixed point behaviour of A, and A' 



in the V-MSSM with Y" = A" = 0. Here (m ,Mi/ 2 ) = (0, 1.5 TeV) is set at the GUT scale, and A'(= A ) 
at the GUT scale is used as an input (horizontal axis). The dependence is drawn for three values of Y'. For 
any input value of A, and A' at the GUT scale, A,/Mq3 and A' /Mq> converge to 0.5 at the SUSY scale 



Msusy '■- 1 TeV as long as Y' > 1, which is nothing but the infrared fixed point behaviour. Fig. 5.5 is that 
for A" in the V-MSSM; here Y' — A' — 0, and A" = Aq are the input values. We can easily see infrared fixed 
point behaviour, and it is similar to that of A'. 

5.1.5 Contribution to the muon g — 2 

The MSSM parameters relevant to the muon g-2 are affected by the extra matters; e.g. from the characteristic 
behaviour of g^ during the renormalization group running. Thus the SUSY contribution to the muon g-2, 
which is reviewed in Sec. |4.3| is quantitatively modified. 

However, it does not change qualitatively. Because the extra particles are assumed to have no direct cou- 
plings with muons, their contribution appears only as higher-order loop level corrections with the assumed- 
to-be-small mixing parameters. Moreover, the particles are expected to be heavy for the SUSY-invariant 



mass Mq' u\e' and the soft masses. Therefore, we can use the same formulae as in Sec. 4.3 for the V-MSSM 
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Q [GeV] Y'(GUT) 



Figure 5.2: (left) Renormalization group flow of Y' in the V-MSSM. We can see that, for a wide range of Y' 
at the GUT scale (which let us say YL m ), Y' falls into ~ 1 at the low-energy scale (which let us say Tg USY ). 
In other words, the renormalization group flow of Y' severely depends on }^ USY . This feature, or the quasi 
infrared fixed point behaviour, can be obvious in the right figure, (right) The value of Tg USY as a function of 
7q UT . Here, the GUT scale is defined with the condition ^i(Mgut) = ^(A^gut)- In both figures, the SUSY 
scale is set to be Msusy = 1 TeV, and tan/? = 20 and Y" = are used. Details of the evaluation procedure 



are summarized in Appendix 4.B 




Figure 5.3: Contour plots of (left) Y' S{]SY and (right) Tg' USY as functions of (Yq VJ , Y^ vj ). The parameters are 
the same as Fig. 5.2 Both Y' and Y" have quasi infrared fixed points ~ 1, and they are nearly independent 
in the plotted region. 
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Figure 5.4: (left) A,/rag 3 evaluated at the SUSY scale as a function of A,/Mi/2, which is an input at the 
GUT scale, (right) A' jrtiQ evaluated at the SUSY scale as a function of A' /M1/2, which is an input at the 
GUT scale. For both figures, the CMSSM scenario is utilized. (m , M1/2) = (0, 1.5 TeV) at the GUT scale, 
and tan/? = 20, sgn/v = +1; A' is set as the same as A ( ), as well as A,, at the GUT scale. Y' is set at the GUT 
scale as Y' = 3 (solid), Y' = 1 (dashed), and Y' — (dot-dashed). For simplicity, F" = A" = is used. 
Note that mg 3 and rag* is the square-root of the SUSY-breaking soft masses squared. From these figure, it 
is obvious that both A' and A, perform the quasi infrared fixed point behaviour if Y' works effectively. The 

infrared fixed points sit at A, 0.5wq 3 and A' 0.5»iq<, which means the stop mixing is rather small in 

the V-MSSM framework (cf. Sec.033). 
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Figure 5.5: A" /m e < evaluated at the SUSY scale as a function of A" jM\ji, which is an input at the GUT 
scale. Similarly to Fig. [54] the CMSSM scenario is utilized with (m , M\ /2 ) = (0, 1 .5 TeV) at the GUT scale, 
and tan/? = 20, sgn/i = +1. Here we set Y' — A' — 0; A" is set as the same as Aq at the GUT scale. Y" is set 
at the GUT scale as Y" = 3 (solid), Y" - 1 (dashed), and Y" = (dot-dashed). Here the infrared fixed point 
behaviour can be seen as Fig. |5.4| 
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Section 5.2 On the Choices We Have Done 

Now we have introduced the V-MSSM. In the above introduction one might wonder why we have chosen 
(10 + 10) pair, and why we will set Y" = 0. For such curious readers, before going into the V-GMSB model, 
the answers for these questions are provided here. 



5.2.1 Why did we choose 10 + 10? 

Extensions of the MSSM with complete SU(5) multiplets are preferred, because there the SU(5)-GUTs are 
respected and the gauge coupling unification is realized. The vector-like insertions, such as 5 + 5 and 10 + 10 
have a virtue for vanishment of gauge anomalies. 

Then one may consider that we can introduce a vector-like pair 5 + 5 instead as extra matters. This does 
not work, because no Yukawa interactions can be additionally introduced. 

Actually an alternate is (5 + 5) + (l + l) model, or LND model 



S = (L'+D'), 



5 = (r + d') , 



1 = N', 



l = N'\ 



(5.22) 



W LND = JZuL'L' + MjyD'D' + Jt N ,N'N' + &"H U L'N' [+^"H i L'N'). 



(5.23) 



As is obvious, the structure of the Lagrangian is almost the same as the V-MSSM. This model corresponds 
to («io,«5) = (0, 1) in Eqs. ( |5.19| >-( |5.21| i, and the gauge coupling unification does realize. One weakness 
of this model is that the Higgs mass increase is smaller than the V-MSSM. The effective Higgs potential is 
almost the same as that in the V-MSSM, but as the leptons, not quarks, form the Yukawa coupling, the color 
factor 3 in the superpotential lacks in this case, i.e., 



Ami 



i-ioop 1 2 ii-ioop v 2 ^' 4 sin/? 4 



= -Am 



4n 2 



log 



Ml 
M 2 F 



5- 



Ml j 



1 - 



M] j 



(5.24) 



Furthermore, it is known that the coupling W has also a quasi infrared fixed point as well as Y', which 
however is slightly smaller: W 0.765 Il93l . As the mass increase is proportional to £^' 4 , this is crucial. 
Therefore, although it is possible to realize = 126 GeV with a large Ms/Mp, we do not consider this 
model in this dissertation. 

On the other hand, we can consider adding more vector-like pairs: 3«io + «s > 3. In this case, however, 
the gauge coupling gj, might blow up before the grand unification; as is shown in Fig. |5.1| g^ increases 
during the running up to the GUT scale in the V-MSSM. If we add more matters, the RGE for g^ becomes 
positive even at the one-loop level (See: Eq. ( |5.19| l). Therefore, this direction is not so promising as is 
naively expected. Note that the upper bound on (3«io + «s) is dependent on other MSSM parameters, in 
particular the SUSY scale M SUSY (cf. Fig.|4~4i. 



It should be also mentioned here that the MSSM with four fermion generations, which can be written as 
"MSSM+10 + 5", is excluded by the Higgs boson searches ED (cf. Ref. ED)- 



5.2.2 Why is Y" = required? 



From the expression ( |5.16[ ), one may notice that the contribution from Y" to the Higgs boson mass is similar 
to that from Y', and thus may wonder why we set Y" = 0. This is due to a large tan/?, which is preferred to 
keep the SUSY explanation of the muon g - 2 problem, and a large //-parameter; as we will see in Sec. 5.3.3 



the yU-parameter is raised by the vector-like quarks during the renormalization group running and thus tends 
to be large in the V-MSSM. 

Under a large tan/?, the mass increase is expanded as 



Ami 



1 1 -loop 
1VMSSM 



3^ 

4n 2 



(5.25) 
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where we defined r :- Mf/Ms (which is usually smaller than one) and £ := /j/Ms ; we set A' — A" — here 
for simplicity. From this expression it is understood that, especially with a large yu-parameter, the interaction 
Y"HdQ'U' just decreases the Higgs mass to disgrace the model. 

This is the reason we simply employ an unnatural simplification of Y" = 0; studies with Y" + are left 
as future works. 



Section 5.3 V-GMSB Model 

Among various SUSY-breaking scenarios, the gauge-mediated SUSY-breaking (GMSB) scenario [ 100 1 is 
notably attractive because of natural suppressions of dangerous flavor-changing processes and CP-violations 
which can appear in soft SUSY-breaking terms. However, the GMSB under the MSSM faces the difficulty 
of explaining the mass of the Higgs boson in the range of 125-126 GeV. This is because the scalar trilinear 
couplings, A Ut d e , arise in two-loop level quantum correction and thus are generally small in the GMSB 
scenario; consequently the mixing in the stop sector, characterized by X, in Eq. ( |4.40| i, is suppressed, and we 
are forced to set the SUSY scale Msusy to be considerably large. Then the SUSY explanation for the muon 
g - 2 problem is spoiled, and the SUSY loses one sales point. 

In the V-MSSM, the extra vector-like matters efficiently increase the Higgs boson mass to 126 GeV with- 
out exploiting large Msusy- This was qualitatively discussed in Sec. |5. 1.3] and will be confirmed quantita- 
tively soon. Thus the GMSB scenario with the V-MSSM is very attractive and phenomenologically viable; 
the V-MSSM resurrects the GMSB scenario. We will see that the parameter space has a region where the 
Higgs boson mass 125-126 GeV is achieved and the muon g - 2 is consistent with the experimental value at 
the 1-2 cr level. 

* * * 



We adopt the simplest GMSB scenario as the SUSY-breaking framework of the V-GMSB model. This 
framework is parameterized by the messenger scale M mess , the soft mass scale A, the messenger number N5, 
the ratio of the Higgs vacuum expectation values tan/? = o u /i>d, and the sign of the yu-parameter sgn//. 



As we saw in Sec. 5.1.4 the gauge coupling constant g^ is marginal to blowing up; as the messengers also 



worsen the situation, we set Ns = 1 in order to preserve the perturbativity of the gauge coupling constants 
up to the GUT scale. We also set fi > 0; note that we have taken the convention under which the gaugino 
masses are positive: M a > 0; thus this choice results in positive contribution to the muon g - 2 (cf. Sec. 4.3 1. 
Here are summarized the parameters of the V-GMSB model and our setting: 

(A, M mess , tan/3, sgn// = +l; M Q = Mv = M E ,{=: My), Y' - 1.0, Y" = o), (5.26) 

where Y' is set at the messenger scale M mess , while My is an input at the low-energy scale Msusy, which is 
defined as Msusy '■— y mj\ 11% ■ Note that the SUSY-breaking parameters specific to the V-GMSB models, A' 
and Bq^ u',E',< are n °t input parameters but yielded by the messengers. 



5.3.1 SUSY-breaking 

In the V-GMSB model, the soft SUSY-breaking terms are induced radiatively via the messenger fields as in 
the ordinary GMSB scenario. The superpotential of the messenger sector is 

W mess = (M D + Td&WdVd + (M L + TlO 2 )^ l$ l, (5.27) 

where Md (Ml) is the supersymmetric mass, and Td CFl) is me SUSY-breaking F-term of the colored 
(non-colored) messengers. The messenger scale is taken as M mess := Md = Ml, and the SUSY-breaking 
F-terms are as To — Tl- Thus, the soft SUSY-breaking parameters are characterized by the two parameters, 
M mess and A := T D I 'Md - TlI Ml, at the messenger scale. Precisely speaking, the relations Md = Ml 
and Td — Tl do not hold generally at the messenger scale due to the renormalization group evolution, even 
if we set Md - Ml and Td = Tl at the GUT scale. Even though, A is almost unchanged during the 
renormalization group running, and hence, the soft mass parameters are not significantly changed under this 
effect. Finally, these messenger fields set the soft SUSY-breaking parameters at the messenger scale M mess , 



and the parameters are evolved down to the weak scale with the RGEs summarized in Appendix 5. A 
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Figure 5.6: The mass spectra of the GMSB (left) and the V-GMSB (right). The GMSB parameters are 
(A,M mess ,tan/?,iV m ess) = (140 TeV, 10 6 GeV,20, 1) in both cases. The SUSY-invariant masses of vector-like 
fields are set as Mq> = Mu> = Me< = 1 TeV for the V-GMSB model. The sfermions of the first and the 
second generations are labeled with #l,<?r, II, and Zr, where the left-right mixings are ignored. 



5.3.2 Characteristics of the V-GMSB model 



As is discussed in Sec. 5.1.4 the couplings Y', A, and A' show quasi infrared fixed point behaviour in the 
V-MSSM scenario. However, in the V-GMSB model, the situation is slightly different; since the parameter 
is input at the messenger M mess , the renormalization group evolution does not work effectively enough to 
drive them to the fixed points. This is important for the Higgs boson mass; the Higgs boson mass decreases 
as Y' is reduced away from Y' = 1 at the messenger scale. Similarly, A, and A' are much smaller than the 
fixed point value especially for smaller M mess as the usual GMSB models. This disables the Higgs boson 
mass to be enhanced by the trilinear couplings, but the constraint from the branching ratio of the b — > sy 
decay is safely satisfied, which is serious in so-called the m/,-max scenario of the CMSSM scenarios 04). 

Another characteristic feature of the V-GMSB model is the size of the Higgsino mass parameter fi at the 
SUSY scale. It becomes very large compared to the usual GMSB models. This is because the vector-like 
matters contribute to the RGE of the up-type Higgs mass squared, which is approximately given as 



dm 2 . 



d\nQ d\nQ 



MSSM 



Q' 



+ m,- 



+ m Hi + 



|Af), 



(5.28) 



where Q is the renormalization scale (cf. Appendix |5.A| >. This contribution, which originates in the vector- 
like matters, gives a large negative contribution to the up-type Higgs mass squared m 2 H . It is as large as that 
from the stop, since Y' - 1 . Hence, the electroweak symmetry breaking conditions require a large value of 
fi. Consequently, the next-to-lightest neutralino^ and the lightest chargino^ tend to consist of the wino, 
and the mass of the lighter stau is likely to be small, because the left-right mixing of the stau mass matrix, 
which is proportional to yUtan/3, becomes large. Especially when tan/3 is large, the lighter stau becomes the 
NLSP. These features can be found in the mass spectrum discussed below. 
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5.3.3 Mass spectrum 

In this subsection, we discuss the mass spectrum of the V-GMSB. In the numerical analysis we utilized 
SOFTSUSY 3.3 [91 1 and FeynHiggs 2 . 9 111 1 II to calculate the mass spectrum of the SUSY particles and 
the Higgs sector. They were modified to take the effects from the extra vector-like matters into account. 



A typical mass spectrum is displayed in Fig. 5.6 where the V-GMSB result is compared to that of the 



ordinary GMSB. Here and hereafter, we assume a common SUSY-invariant mass for the vector-like fields, 
M v = Mqi = Mjj> (= M E i). As is explained just above, the heavier chargino xt an d the heavier two 
neutralinos ^ 4 mainly consist of the Higgsino, and are much heavier than the ordinary case. In addition, 
it is found that the squarks become heavier; this is because the gauge coupling gi stays large during the 
renormalization group evolution. On the contrary, a ratio of the gaugino masses is less affected by the 
extra matters. This is because the ratios of M a /a a , where M a are the gaugino masses, are fixed during the 
renormalization group evolution at the one-loop level, and the gauge coupling constants g a are (should be) 
set at the low-energy scale. The masses of the vector-like squarks are close to those of stops and sbottoms 
up to the SUSY-invariant mass, because they have the same quantum numbers as the corresponding squarks. 



Let us explain how the mass spectrum is evaluated in SOFTSUSY and what kinds of modifications are ap- 
plied. First of all, the program estimates the Standard Model gauge couplings a a (m z ) and Yukawa couplings 
T,(;« z ) at the scale of the Z-boson mass m z 1 102 1 . They are evolved upwards from m z to the messenger scale 
A^mess with solving the V-GMSB RGEs. In the numerical calculations, the SOFTSUSY package was modified 
to include the two-loop level RGEs of the V-GMSB model. It is also important to include threshold correc- 
tions and self-energy corrections of vector-like matters to the gauge coupling constants at the m z scale, since 
they can give ~ 10% contributions especially to the coupling constant gj. Those corrections also affect the 
gaugino masses and the scalar masses squared generated at M mess through the gauge coupling constants. 

Next, the soft SUSY-breaking parameters are provided by the messenger loops at M mess . The extra 
Yukawa couplings are also set at the scale, which are T'(M mess ) = 1 and Y"(M mf , ss ) = 0. The gauge, 
Yukawa and soft parameters are evolved with the V-GMSB RGE down to the SUSY scale, Msusy, which is 
determined by the stop masses as Msusy = 

The masses of the superparticles are evaluated including the whole one-loop level corrections to the self- 
energies within the MSSM. The pole mass of the gluino also receives a correction from vector-like matter 
loops. However, this turns out to be around a few GeV, since the masses of the vector-like matters are 
close to Msusy, and their contributions to the self-energy are relatively small. Similarly, contributions of the 
vector-like fields to the electroweak gaugino masses are safely neglected. 

Then, the Higgs potential is investigated, which determines the Higgs boson mass and the /i-parameter. 



The MSSM part of the Higgs boson mass is evaluated at the two-loop level by the FeynHiggs package 1 101 1. 
The contribution from the vector-like matters is estimated with the one-loop level effective potential dis- 
cussed in Sec. 5.1.3 where one-loop level contributions from the vector-like matters are taken into account, 
and is added to the MSSM Higgs boson mass in quadrature to obtain the V-MSSM Higgs boson mass. The 
two-loop level contribution of the vector-like matters can shift the Higgs boson mass by ~ 1-10 GeV, study 
for which is reserved for future works. On the other hand, the yu-parameter is less affected by the vector- like 
matters except for the renormalization group evolution ( |5.28| l, because tree-level contribution to /j is fairly 
large in the V-MSSM. 



Section 5.4 Vacuum Stability Bound 

One of the most severe constraints on the V-GMSB model is the vacuum stability bound ifTBI . When tan/? 
is large, the large trilinear coupling of the staus and the Higgs boson can generate charge breaking global 
minima and destabilize the electroweak symmetry breaking vacuum [ 103 ]. Therefore, tan/? has a tight upper 
bound for given M mess and A in the GMSB framework from the condition that lifetime of the proper vacuum 
must be longer than the age of our Universe. 
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Let us see this condition in detail [104]. The MSSM (and the V-MSSM) Higgs potential together with 
the stau sector is expressed as (cf. Eq. (|4.8|l) 



(5.29) 
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Note that here the one-loop level correction to the quartic coupling, which actually can be found in Eq. ( |4.32| >: 
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(Avt 1 ))^— Llog-i, (5.30) 

is included as Sh, whose definition is 
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The terms with are also ignored because \H^\ is tiny in our concerning large tan f3 region, although they 
give considerable corrections [ 105 1. 

Hereafter we simply assume fi e R as well as Y T > 0. Expanding the expression with H® = v + h, we 
obtain the following mass terms for the stau: 
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(5.32) 



The vacuum stability is at classical level understood that the eigenvalues of the stau mass matrix must be 
positive 111041 ; that is, 



(m T ■ fi tan/?) 2 < i (2m 2 Li + 2m 2 w - m§) (m 2 ^ + m| - m 2 ,,) . 



(5.33) 



This condition is certainly not sufficient because it just asserts that "our" vacuum, the EWSB vacuum, is a 
local minimum. The condition which must hold is that the EWSB vacuum is the global minimum, or, as 
a more weaker condition, it is just a local minimum but the transition rate to other minima, especially the 
global minimum, is longer than the age of our Universe. For the latter case, the transition rate is estimated 
by a semi-classical method, searching for so-called bounce solutions B 1 061 . In Ref. II104L an approximate 
formula for the bound on //tan/? is obtained by using multi-dimensional bounce configurations, including 
top-stop radiative corrections to the Higgs potential. 



In a work of Author" 3 , Ref. Ifl6l . this vacuum stability condition on the V-GMSB model is discussed, 
where a fitting formula given in Ref. 11041 vl on arXiv| was utilized: 

H tanyS < 76.9 y/m L} m^ + 38.7 [m Li + m^) - 1.04 x 10 4 GeV. (5.34) 
Very recently, however, the fitting formula was revised as follows Ml 041 v2 on arXiv | : 

/ztanjS < 213.5 ^m L ,m Ei - 17.0 (m L , + m E] ) + 4.52 x 10~ 2 GeV" 1 (m Ly - m E3 f - 1.30 x 10 4 GeV. (5.35) 

In these expressions, m El and m Ei are the square-root of the SUSY-breaking soft mass squared for, re- 
spectively, the left- and the right-handed stau. Moreover, also very recently, the vacuum stability bound 
was reanalyzed with including effects of a radiatively-corrected tau Yukawa coupling [ 105 1 . The results in 
Ref. 1 105], however, cannot be directly applied to the V-GMSB model, because the masses and the mixing 
angle of the staus are different from those in Ref. 11051 . 

3 Done in collaboration with Dr. M. Endo, Prof. K. Hamaguchi, and Dr. N. Yokozaki. 
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The above result ( |5.35| > is obtained in the limit of zero temperature, but thermal effects may tighten the 
bound 1 1071 11081 . The thermal decay rate of a false vacuum is usually estimated with following the method 
in Ref. II 1 091 . Evaluating the Higgs potential at the one-loop level with including the thermal effect coming 
from the top quark and the electroweak gauge bosons, one can find that the stability bound can become more 
severe by ~ 10% than the zero temperature result for a small stau mass region [107]. In addition to the 
thermal corrections, one should note that the proper vacuum may be required to be the global minimum if 
the vacuum expectation values of the scalar fields stayed away from the ordinary ones in the early universe. 

Considering all these various factors together, this dissertation uses the bound of Eq. ( |5.35| > to draw the 
upper-bound, and also the bound but weakened by 10% is drawn for reference use. 

* * * 

Finally, let us touch on another possibility of the vacuum instability. In a class of SUSY models, the 
trilinear coupling of the top squark is predicted to be large in order to raise the Higgs mass. Such a large 
trilinear coupling may spoil the vacuum stability in the stop-Higgs plane similarly to the stau-Higgs plane 
discussed above. The meta-stability bound was studied in Ref. 1 1 1 1 and obtained as 

A] + 3ii 2 < 7.5 (m| + m|) . (5.36) 

However, as the A-terms are not so large in the V-GMSB model as ordinary GMSB models, this condition 
does not appear in our discussion. 



Section 5.5 Higgs mass, muon g - 2 and vacuum stability in V-GMSB 

Now we are ready to discuss numerical results: we are about to see that the V-GMSB model explains the 
126 GeV Higgs boson mass together with explaining the muon g - 2 anomaly. 



The current status of the V-GMSB model is summarized in Fig. 5.7 together with the current LHC 



bounds discussed in Sec. |5.6| (See also Fi g.|5.12| ) The mass spectrum including the Higgs mass is calculated 
with the procedure summarized in Sec. |5. 1.2] where SOFTSUSY3.3 ll9D and FeynHiggs2.9 II 1 01 1 are 
utilized. The muon g - 2 is calculated with FeynHiggs. Y' = 1.0 is set at the messenger scales, Mq = 
Mu> = Me> =: My is set at the low-energy scale, and Y" = A" = is assumed. The values of important 
Standard Model parameters are set as a s (mz) = 0.1 184 and m t = 173.5 GeV. 



In Fig. 5.7 contours of the Higgs boson mass and the muon g -2 are drawn in the (trig, tan/?)-plane for 
the messenger scales of (a) M mess = 10 ll) GeV, (b) 10 8 GeV, and (c) 10 6 GeV. Here, rrtg is the gluino pole 
mass, which is mainly determined by A, and tanyS is an input evaluated at the electroweak scale as usual. 
As the manner discussed in Sec. 5.4 the vacuum stability bound is imposed with Eq. ( |5.35[ ). In each of the 



figures, it is drawn with the blue dot-dashed line; the region above the line is disfavored. Above the line, the 
blue double-dotted long-dashed line is drawn as the vacuum stability condition which is weakened by 10%. 
The black solid and dashed lines show the LHC constraints, which we will discuss in detail in Sec. |5.6| The 
NLSP is the lightest neutralino below the light blue line, and the lighter stau above the line; the LSP is the 
gravitino in this model. In the gray shaded region, SOFTSUSY fails to calculate the stau mass. Note that such 
a parameter region is experimentally excluded as long as the NLSP (the lighter stau) is long-lived. 



In the green bands in Figs. 5.7 'a)-(c), the Higgs boson mass takes a value of 125-126 GeV for reference 



values of SUSY-invariant mass of the vector-like matter of My — 1 TeV and 1.2 TeV. The Higgs boson 



mass increases as the vector-like matter becomes lighter (cf. Sec. 5.1.3 1. Actually, the Higgs boson mass of 
125-126 GeV can be realized in the whole parameter region of Fig. |5.7| by changing the vector-like matter 
mass. For fixed My, the Higgs boson becomes lighter when the gluino mass is smaller, because the radiative 
corrections to the Higgs potential from stops and vector-like stops decrease. On the other hand, for higher 
messenger scale, the stop mass becomes larger during the renormalization group evolution, and thus the 
Higgs boson mass of 125-126 GeV is realized with a smaller gluino mass. 

On the contrary to the mass of the Higgs boson, masses of the MSSM superparticles and the SUSY 
contributions to the muon g - 2 are not sensitive to My. We have checked that the MSSM superparticle 
masses change only by < 2% when My is varied from 500 GeV to 1 TeV, and the muon g - 2 changes by 
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Figure 5.7: (a)-(c) Contours of the Higgs boson mass, the muon g - 2 and the LHC constraints in the V- 
GMSB model are shown for (a) M mess = 10 10 GeV, (b) 10 8 GeV, and (c) 10 6 GeV. In the green bands, 
the Higgs boson mass is 125-126 GeV for My - 1 GeV and 1.2 GeV. The yellow (orange) regions show 
the parameter spaces where the muon g - 2 discrepancy is explained within 2<x (lcr)-level by the SUSY 
contribution. The light blue lines indicate m^o = ntf l ; the NLSP is the lightest neutralino below the lines, 
and the lighter stau above them. The regions above the blue dot-dashed lines are constrained by the vacuum 
stability condition ( |5.35| i; as a reference the bound which is weakened by 10% is drawn with the blue double- 
dotted long-dashed lines. The LHC constraints, discussed in Sec. 5.6 are drawn with the black solid lines 
without theoretical uncertainties, where the regions left to the lines are excluded. The black dashed lines 
indicate the LHC bounds with theoretical uncertainties: 35% error of the production cross section is adopted 
for the neutralino NLSP, and 2% error in terms of the stau mass is for the stau NLSP. The gray shaded 
regions are excluded by the failure of the SOFTSUSY calculation, (d) Masses of the SUSY particles under 
A^mess = 10 6 GeV. The green contours are for the lighter stau, and the red lines for the lightest squark among 
the first two generations. The light blue line indicates m^o = ntf l . 
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less than 1%, correspondingly. In the figures, we fix My — 1 TeV in the calculations of all the quantities 
except for the Higgs boson mass. 

In the yellow (orange) regions in Figs. 5.7 'a)-(c), the muon g-2 discrepancy, Eq. ( |2.24[ ), is explained at 
the 2cr (lcr) level, i.e., 



10.1 xl(T 10 < [Aa^lsusY < 42A xl0~ 10 in the yellow regions, 
18.1 xlCT 10 < [Aa^JsusY < 34.1 xlCT 10 in the orange regions. 



These regions are sensitive to both the gluino mass, i.e., the soft mass scale Msusy, and tan/3 (cf. Sec. 4.3 1 



It should be stressed that the gluino mass, and thus Msusy, has an upper bound for successful explanation 



of the muon g - 2 under the vacuum stability condition, which we will discuss again in Sec. 5.6.4 (see 



Fig. 5.12| i. This bound becomes tighter for a larger M mess . The principal reason is that the renormalization 



group evolution, which raises the soft scalar masses during the running down from the messenger scale to 
the low -energy scale, works more significantly under a larger M mess . In addition, as the gauge couplings 
are larger at higher energy scales in this model (cf. Fig. |5.1| l, a larger M mess results in slightly heavier soft 
scalar masses at M mess . Thus heavier squarks and sleptons are provided for a larger M mess , and the SUSY 
contributions to the muon g - 2 are suppressed, in a fixed value of the gluino mass. 

As we shall see in the next section, LHC searches for the SUSY particles depend on the species of the 



NLSP. In Figs. 5.7 'a)-(c), the NLSP is the lightest neutralino (the lighter stau) below (above) the light blue 
lines. The stau NLSP region becomes smaller as M mess increases, for the renormalization group running, 
which increases the sfermion masses relative to those of the gauginos, works more significantly. For SUSY 
searches at the LHC, the masses of the gluino, squarks, and the stau are particularly important. For illus- 



tration, we show in Fig. 5.7 d) the masses of the lighter stau T\ and the lightest squark among the first two 
generations with green and red lines, respectively, for M mess = 10 6 GeV. The squark mass is almost inde- 
pendent of tan/3, since it is governed by the strong interaction. On the other hand, the stau mass, which 
determines the LHC constraint in the parameter region with the stau NLSP, depends on the gluino mass (the 
soft mass scale), and on tan/3 through the tau Yukawa coupling Y T » (m T /v) tan/3. 

An interesting feature of the V-GMSB model is that the mass parameter My should be < 1 .2 TeV in 
order to explain the muon g — 2 discrepancy at the 2<x-level together with the 126 GeV Higgs boson. This 
behaviour can be understood from Eq. ( |5.18| i; for larger Mf(= My) in the expression the SUSY-breaking 
mass terms, M s must be larger to realize enough increase, which means the masses of the sleptons and the 
gauginos get larger, and the shift of the muon g - 2 decreases. Therefore, in this context the vector-like 
quark mass should be < 1 .2 TeV, and this fact makes searches for vector-like quarks very interesting. We 



will discuss this topic in Sec. 5.7 after the discussion on the LHC bounds drawn in the figures. 



Section 5.6 LHC Bounds on the V-GMSB Model 

In this section the constraints on the V-GMSB model from LHC SUSY searches are discussed. Actually 



the bounds are already shown in Fig. 5.7 here the procedure to obtain those bounds and interpretations of 
them are provided. The V-GMSB model is also constrained by searches for the vector-like quarks, but the 
constraints are discussed in Sec. 15.71 

5.6.1 Overview 

In the V-GMSB model, the gravitino is the LSP, and collider signatures are determined by species and the 
lifetime of the NLSP. As discussed in the previous section, the NLSP is either the lightest neutralino or the 
lighter stau in the V-GMSB model, which are respectively realized below and above the light blue lines in 



Fig. 5.7 Let us summarize the cases. 



Long-lived Neutralino NLSP: The promising signature is multi-jets plus large missing energy (Z?x) yielded 
from pair-production of colored particles, where the multi-jets are provided by cascade decays of the 
colored particles and the missing energy is from the NLSPs. 
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The ATLAS collaboration reported a search for this channel using the data corresponding to an inte- 
grated luminosity of 5.8 fbr 1 obtained at Ecm = 8TeV in Ref. (ST|, and obtained a lower bound on 
the gluino mass for the case where the squarks are heavy, which is expected in the V-GMSB model, as 
ntg > 900 GeV. The CMS also analyzed their (7 TeV, ~ 5 fbr 1 ) data to obtain similar results 



Neutralino NLSP with prompt decay: In this case the above signature, multi-jets plus large missing en- 
ergy, is still expected. Nevertheless, as the NLSP can decay as ^ — > y + G, di -photon signature with 
large missing energy is promising. 

For the decoupled squark scenario, the CMS collaboration excluded ntg < 1.1 TeV for the bino-like 
NLSP case, and ntg < 750 GeV for the wino-like NLSP case, with (8 TeV, 4 fbr 1 ) data ifTTTl . and also 
obtained similar results from (7 TeV, 4.3 fbr 1 ) data Ml 121 : both analyses are based on jet(s) + photon(s) 
+ Xt- Note that the NLSP is bino-like in the V-GMSB model for the large yu-term. 

The ATLAS collaboration excluded gluinos with ntg < 1.1 TeV (< 950 GeV) for the bino- (wino-) like 
NLSP case in the limit where the squarks are decoupled; these bounds are based on jet(s) + di-photon 
+ Xt signature with (7 TeV, 4.8 fbr 1 ) data [l 13 1. They also employed the di-photon search without jet 
requirement 111 131 ; this search focuses the electroweak production of the SUSY particles: pp — » xx 
etc., and is interesting because it covers the case where the colored particles are extremely heavy and 
not accessible with the 14 TeV. 

If the neutralino is quasi-long-lived and decays during the flight in the detector, typical signature at 
the LHC would be a non-pointing photon 1 1 14 1 and a neutralino in-flight decay into a Z-boson 1 1 15 1. 
Since the V-GMSB has an upper bound on the gluino mass ntg < 1.2 TeV (1.8 TeV) in the neutralino 
NLSP case once a 125-126 GeV Higgs and the muon g - 2 constraints at lcr {2(f) are imposed, it is 
expected that a large part of the parameter space can be judged also in this case. 

Long-lived stau NLSP: In this case the NLSP does not yield large Xt, but is observed as a heavy stable 
charged particle (HSCP), where the term "stable" is used from the experimental viewpoint, i.e. it 
means long-lived enough to escape from the detectors. 

The CMS collaboration reported a lower bound on the mass of the lighter stau as 

rrtf, > 223 GeV (5.37) 

with analyzing their data of (7 TeV, 5.0 fbr 1 ) lfTT6 l. where the (lighter) stau is assumed to be produced 
only with the direct production, pp — > t{t*. It should be emphasized that this bound is generic, i.e., is 
most conservative and can be adopted to any models, because it targets the lighter stau pair-production; 
one should however note that the production rate depends on the stau mixing. 

The CMS collaboration also published constraints on the stau mass for typical GMSB models, where 
charginos and neutralinos involve the main production channels HI 161 . However, it cannot be directly 



applied to the V-GMSB case, because the mass spectra are different (cf. Fig. 5.6 1. The difference 
results in a deviation of the velocity distribution of the NLSP staus. Once the reconstruction efficiency 
of the HSCPs is published, the production channels of the charginos and neutralinos can be considered, 
which would provide much more strict constraints. 

The ATLAS collaboration also published a result on searches for the HSCPs 11 171 . but their search 
seems to assume direct productions of the whole sleptons, i.e., the production is not limited to the 
lighter stau, and the reconstruction efficiency of the particle is not provided either. Therefore, their 
results cannot be applied to the V-GMSB. 

Stau NLSP with prompt decay: In this case the promising channel is multi-taus + Xt- The ATLAS col- 
laboration, utilizing their (7 TeV, 4.7 fbr 1 ) data, searched for the events with > jet + > 1 tau + 
0-1 lepton + Xt, and constrained the parameter space of GMSB scenario with (M mess , N mess , sgn/z) = 
(250 TeV, 3, +) to obtain the bound of ntg > 1 .0 TeV, which corresponds to ntf > 160 GeV-70 GeV de- 
pending on tan/? (of 2-60). Nevertheless, this result cannot be applied to the V-GMSB model because, 
especially, the stau is relatively heavier in the V-GMSB model as is already discussed. 

An in-flight-decay of a stau leaves kink signature in the TRT detector, and is expected to be observed 
if the decay length is ~ 1 m HI 181 : this feature is the same as the ordinary GMSB scenario. 
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In this dissertation we concentrate on the long-lived NLSP scenario. The LHC exclusion limits in Fig. [577] 
are drawn under this assumption, i.e., they are based on the analysis for the first and the third case in the 
above listing. In the following of this section the detailed procedure in the evaluation of the bounds is 
explained. 

5.6.2 LHC constraint for the V-GMSB with long-lived neutralino NLSP 

For the region with a neutralino NLSP, we interpret the ATLAS result of the search for the superparticles 
in events with no lepton, 2-6 jets and missing energy, with data obtained at the LHC with Eqm = 8 TeV 
corresponding to an integrated luminosity of 5.8 fbr 1 lISTI . as a constraint for the V-GMSB model. The 
detailed explanation of the Monte Carlo analysis is given here. 

Data sample and triggering 

The ATLAS collaboration utilized their data taken in 2012 corresponding to a total integrated luminosity of 
5.8 fb obtained at the 8 TeV LHC. Their trigger requirement is a jet with pi > 80GeV and the missing 
energy / Ej > 100 GeV.* 4 Since the event selection, discussed later, includes requirements of a jet with 
pi > 130 GeV and missing energy of^x > 160 GeV, the trigger efficiency is at the efficiency plateau as the 
ATLAS collaboration wrote "full efficiency" in the report. Therefore, in our analysis the trigger efficiency is 
not considered. 

Our data are obtained with MadGraph 5 [ 1 19 1 package. It includes Pythia 6.4 1 120 1 to simulate parton 
shower and initial- and final state radiation (ISR and FSR), and Delphes 2 . [ 121 1 for detector simulation; 
both of them are utilized in our analysis. The parton distribution functions (PDFs) are obtained from the 
CTEQ6L1 set 01221 . Pythia setting is based on the MadGraph default. 

Simulated events 

The events in pp — > 'gg,'gq ( *\'^* ) q ( * i channels, which are relevant for the SUSY signals, are generated, but 
among them the production channels with t and b are neglected, whose production cross sections are less 
than a few % of the total cross section. 

The V-GMSB mass spectrum is generated with the customized SOFTSUSY 3 . 3 l9T| as is explained in 
Sec. 5.1.2 and passed to SUSY-HIT 1.3 |123| to calculate the decay pattern of the SUSY particles and 



the Higgs bosons. The generated events are normalized to the NLO cross section obtained with Prospino 
2 . 1 Ifl24l [12511 . where the CTEQ6L1 and the CTEQ6 . 6M PDFs lfl22l are used. 

Object reconstruction 

Our detector simulation is based on the Delphes 2 . package. The default parameter set for the ATLAS 
experiment, which can be found in DetectorCard_ATLAS.dat in Delphes package, is used for energy 
resolutions and calorimeter tile configurations. 

Jets are reconstructed using the anti-£ f jet clustering algorithm |44] with the distance parameter of 0.4, 
where Delphes utilizes FASTJET 2 . 3 [ 126]. For the missing energy, the Delphes vanilla outputs are used. 

In the electron reconstruction, the ATLAS collaboration uses the "medium" criterion. The detection 



efficiency is reported in Ref. [48 1, which are already introduced in Sec. 3.2 However, the report is based on 
the data taken at the 7 TeV LHC, and no efficiencies at the 8 TeV LHC are public. Therefore, assuming that 
the efficiencies at the 8 TeV LHC are the same as that at the 7 TeV, the efficiency used in the analysis is taken 
from Ref. [48 1. In principle, the efficiency is, and should be expressed as, a function of (p T ,j]). However, 
such functions are not reported yet due to limitation of the data; they instead reported the efficiency as two 



4 The transverse momentum pj and the transverse missing energy are defined as 

P7 :={px,Pv), Pt:=IW; fa:=(tf x ,fa), fa = (5.38) 
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separated functions of pj and of 77 with unignorable uncertainties. Thus, for simplicity, following efficiency 
is adopted for the electron reconstruction in our analysis: 

• The reconstruction efficiency as a constant: ff reco = 0.943 [48 Fig. 22(B)]. 

• The identification efficiency for "medium" as a constant for pi > 20GeV: 6j™ ed = 0.942 11481 Table 7]. 

The identification efficiency significantly drops below p T = 20GeV; for such electrons the 
efficiency is set as (a rcco , e^ & ) = (0.943,0.8) in our simulation set, although q- 1cc0 is not reported 
for this region. Note that, however, this efficiency is not used in this analysis, because electrons 
below p T = 20 GeV are discarded. 

The muon in this analysis seems to be reconstructed with the "combined" method, which is explained 



in Sec. 3.2 The corresponding efficiency is found in Ref. [53 1 for muons with pj > 20 GeV, which is also 
based on the 7 TeV data. Our analysis adopts the following simplified efficiency for the muon detection: 

• The efficiency for inner detectors as a constant: e lD = 0.88 ll53l Fig. 4]. 

• The efficiency for muon spectrometer in combined muon method as: = 0.973 for \rj\ > 0.25 and 
0.828 for < 0.25 [53, Fig. 6]. This difference comes from the fact that the muon spectrometer is 
only partially equipped in this region. 

Consequently, the lepton reconstruction efficiencies are defined as 

e(e) = areco • * 89%, e(p) = e ro • * 96% (82%) for \tj\ > 0.25 (< 0.25). (5.39) 

in our analysis. 

Actually these efficiencies are prepared for another analysis; in the analysis the events with three 
leptons are simulated and analyzed, and this set of efficiencies reproduces the ATLAS analysis in 
Ref. fT27l very well. 



Only the jets with p T > 20 GeV and \rj\ < 2.8, the electrons with p T > 20 GeV and I77I < 2.47, and the 
muons with pj > 10 GeV and \rj\ < 2.4 are considered. For overlap removal, jets are rejected if electrons are 
found within a distance of AR = 0.2 from the jet, and then leptons within AR = 0.4 of any jets are discarded. 
Here, pi is a missing transverse momentum, AR := ^{Aij) 2 + (A<f>) 2 is a distance parameter between two 
objects, and r\ and <p are pseudo-rapidity and azimuthal angle around the beam direction, respectively. 

Signal regions 

Our analysis is employed with the same definition of the 12 inclusive signal regions (SRs) as the original 
analysis I8TI . which is summarized in Table |5T2"| However, the selections based on jet quality selection 
criteria and primary vertex reconstruction are not employed due to the simplified detector simulation. 

The events with any reconstructed electrons or muons are vetoed. At least two jets with pj > 60 GeV 
are required, and the leading jet must be with pj > 130 GeV. The SRs are classified with the number of jets 
having pj > 60 GeV. The definition is inclusive; for example, if an event has four jets with pj > 60 GeV, it 
can be a candidate for the events in the SRs A, B, and C. 

Two requirements designed to reduce the background events from multi-jet processes are impo sed. One 
is on the azimuthal separations between }Zi and jets, which is expressed as A</>(j, v # T ) m j n in Table 



5.2 



For 

the SR A (B), leading two (three) jets are considered for this selection. For the others, leading three jets 
must have a separation of at least 0.4, and all the other jets with pj > 40 GeV must at least 0.2. The other 
requirement is on the ratio between Jtj and trr"2, which is defined as 

(») 



+ \\p T of ('-th leading jet| I 



1=1 



Finally, the inclusive effective mass is utilized, which is defined as m 1 ^ :- m^, where Af is the number 
of jets with pj > 40 GeV. 
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Table 5.2: The definition of the 12 inclusive signal regions (SRs), which is the same as the ATLAS original 
analysis lISTl . The effective mass m "^ (m 1 ^) is defined as the scalar sum of^x and pj's of the leading «-jets 
(all jets with pj > 40GeV). Some categories have two or three SRs, which are shown in the final two rows; 
they are called 'tight', 'medium', and 'loose', respectively. 





Signal Regions 




A 


B 


C 


D 


E 




(> 2-jets) 


(> 3-jets) 


(> 4-jets) 


(> 5-jets) 


(> 6-jets) 


# leptons 


= 


Xr [GeV] > 


160 


M/i)[GeV] > 


130 


p T (j 2 )[GeV] > 


60 


M/'sMGeV] > 




60 


60 


60 


60 


M/ 4 )[GeV] > 








60 


60 


M/sHGeV] > 








60 


60 


p T (j 6 )[GeV] > 










60 


A0(;' ; ,XT)min > 


0.4(/ = 1,2) 


0.4(j =1,2,3) 


0.6 (i < 3); 0.4 (pr > 40 GeV jets) 


AM« > 


0.3/0.4/0.4 
(n = 2) 


0.25/0.3/ — 
(n = 3) 


0.25/0.3/0.3 

(B = 4) 


0.15/ — / — 

(B = 5) 


0.15/0.25/0.3 
(n = 6) 


m* c [TeV] > 


1.9/1.3/1.0 


1.9/1.3/ — 


1.9/1.3/1.0 


1.7/ — / — 


1.4/1.3/1.0 
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Table 5.3: The result of the ATLAS analysis in Ref. ||8T| . The last two columns show the 95% upper limits 
(UL) on the excess number of events, Nbsm, and that on the cross section of the new physics, <x BSM . One 
should note that the experimental uncertainties are already considered in the calculation of the upper limits. 



SR 


Background 


Observed 


UL on jVocM 


UL on ctrcm [fbl 


A-loose 


650 + 130 


643 


224.8 


38.8 


A-medium 


140 ± 33 


111 


33.9 


5.84 


B-medium 


115 ± 30 


106 


43.8 


7.55 


C-loose 


155 ± 31 


156 


65.7 


11.3 


C-medium 


33 ± 8 


31 


17.9 


3.09 


E-loose 


5.7 ± 1.7 


9 


10.4 


1.79 


E-medium 


3.5+ 1.7 


7 


9.9 


1.71 


A-tight 


14 ± 5 


10 


8.9 


1.53 


B-tight 


8.7 ± 3.4 


7 


7.3 


1.26 


C-tight 


2.8 ± 1.2 


1 


3.3 


0.57 


D-tight 


6.3 ± 2.1 


5 


6.0 


1.03 


E-tight 


10+ 4 


9 


9.3 


1.60 



Analysis procedure 



The ATLAS collaboration found no significant excesses; therefore the cross section of models beyond the 
Standard Model, ct B sm, receives an upper bound. The expected number of the background events, the 



observed number, and the upper limits are summarized in Table 5.3 This result is interpreted to obtain the 



constraints (shown in Fig. 5.7 1 with the following procedure, so-called CL s -method. 
First, the expected sensitivity CLs Xp is defined for each of the 12 SRs. It is defined as 



CLP 



cC p b 

CL^ xp ' 



(5.40) 



CL^ xp 
CC b 
N h 
AU 

0" s+ b 

mo-) 



1 - (the probability that a random variable which obeys f{Nb, cr b ) exceeds Nb) , 
1 - (the probability that a random variable which obeys f(N s+ b, c s+ b) exceeds Nb) . 
(expected number of background events) , 
Nb + (expected number of SUSY events) , 
(uncertainty of Nb) , 



^jcr^ + erf, where cr s is the uncertainty of the number of SUSY events, 
Poisson (Normal (ju, a)) , 



where Poisson(x) is the Poisson distribution with the mean x, and Normal(//, cr) is the normal distribution 
with the mean fi and the variance cr 2 . As the uncertainties in the cross section of the SUSY events are 
not included in the procedure to obtain the upper limits (shown in Table |5.3j l but are considered afterwards 
(cf. Fig. 5.9 etc.), our procedure uses <x s = 0. Here one should note that CL S p is determined for each model 
point, and that the number of background events is assumed for simplicity to distribute with the normal 
distribution. 
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CL" P and CLf? can be interpreted as the expected confidence level of the background-only and 
the signal-plus-background hypothesis, respectively, or the expected probability that the respective 
hypothesis describes the data. Here, since the "expected data" are nothing but the expected number 
of background events, CLjj xp and CL^ P are defined as the probability that the respective hypothesis 
describes the number of background. 

The SR which gives the smallest Cl4 xp is expected to give the most stringent constraints on the model 
point, and therefore adopted as "the SR for the point?' In Figs. 5.9 and 5.10 which will be introduced 
later, which SR is selected for each model point is displayed. Then, the observed confidence level CL S is 
calculated for the selected SR, which is defined as 



CL S := 

CL b 

CL s+b 

Ndata 



CL s+b 
CL b 



(5.41) 



= 1 - (the probability that a random variable which obeys /(M>, o"b) exceeds iVdata) , 
= 1 - (the probability that a random variable which obeys /(iV s+ b, c s+b ) exceeds iVdata) > 
= (observed number of events) . 



If the CL S , calculated for the selected SR, satisfies CL S < 0.05, the point is excluded at 95% confidence 
level (CL). Otherwise, the point is not excluded. It should be important that, although we have 12 SRs, the 
analysis based on the data is employed for a sole SR, and thus the condition CL S < 0.05 results in 95% CL 
limit. 

It is known that the CL s -method generally gives a conservative limit, i.e., the obtained constraint 
is looser than the "true" one. 



Verification of our analysis 

To verify our analysis, the distributions of the inclusive effective mass, m m S, are compared. The distributions 
for the SRs C-tight, D-tight and E-tight are displayed in Fig. 5.8 where the results from ATLAS Monte Carlo 



simulations and our Monte Carlo simulation are shown, in which all the selections but on mS£ itself are em- 
ployed. As the ATLAS collaboration provides only the histograms which include both the Standard Model 
and the SUSY contributions, the comparison was done between the ATLAS Standard Model simulation plus 
the ATLAS SUSY simulation (displayed with gray plus blue histograms), and, the ATLAS Standard Model 
simulation plus our SUSY simulation (displayed with red lines). The results do agree well with slight deficit 
in our analysis. 



The exclusion limit on the CMSSM (mo, Mi/j) plane from our analysis is also shown in Fig. 5.9 This is 
based on the procedure, the CL s -method, described just above; the letters in the parameter spaces denote the 
expected-to-be-most-sensitive SR, that is, the SR which is selected for the point. Interpolation to obtain the 
curve is based on log CL S . The constraint is slightly looser than that reported by the ATLAS collaboration, 
but as ours is still conservative, we decided to accept this analysis. 



Result 

The results are already shown in Fig. 5.7 the black solid lines in the neutralino NLSP region, i.e., below 



the light blue lines. The exclusion limits with considering theoretical uncertainties are also drawn as the 
black dashed lines. The theoretical uncertainties mainly originate in the evaluation of renormalization and 
factorization scales, and choice of the PDFs 1128L For the uncertainties we adopt ±35% in the production 
cross section. 

The left sides of the black lines are excluded. Consequently, the gluino mass is required to be larger 
than approximately HOOGeV for M mess = 10 6 GeV, if the theoretical uncertainties are not included. The 
theoretical uncertainties can shift the mass bound by < 100 GeV. When M mess is larger, the exclusion 
becomes weaker for the fixed gluino mass. This is because the squarks become heavier, and the production 
cross section of the squark, especially that of pp — > ~gq, becomes smaller. When the messenger scale is as 
large as 10 10 GeV, the bound becomes nig > 1030 GeV. 
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C-tight D-tight 




500 1000 1500 2000 2500 3000 3500 4000 



meff(inc) [GeV] 

Figure 5.8: The rrtfg distributions before the selection on itself. Here the lepton efficiencies are taken 
into account (cf. Fig. |5.1 1] ). As a benchmark model, the CMSSM with (mo, M\j%, Aq, tan/3, sgn//) 
( 1 600 GeV, 400 GeV, 0, 10,+) is chosen. Only the distributions for the SRs C-tight, D-tight, and E-tight 
are shown here. The results from the ATLAS Monte Carlo analysis [81] are shown as the histograms; the 
gray are the events from Standard Model background, and the blue, stacked on the gray, are from the SUSY 
signal. The results based on our analysis are drawn with red lines, but since our analysis does not simulate 
background processes, the background results reported by the ATLAS collaboration are used as the back- 
ground histogram for our analysis. These two results are compared to certificate our analysis. Note that the 
observed data are not shown here. 
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mSUGRA: tan(3=1 0, A =0, (a>0 [ATL-CONF-201 2-1 09/Preliminary] 

800 n - 1 — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r~ 




m n [GeV] 



Figure 5.9: 95% CL exclusion limit for the CMSSM scenario with (Aq, tan sgn fi) - (0, 10, +) presented in 
the (mo, M^Vplane. The black line with a gray band is the result reported by the ATLAS collaboration in 
Ref. [81] with ±lcr uncertainty on the calculation of SUSY cross section. The red line is the result from our 
analysis, and the letters denote the SR which is adopted for the point. (See text for detail.) As a reference, 
the limits obtained in our analysis with +20% theoretical uncertainty are shown as the red dotted lines. 
Interpolation to obtain the curve is based on logCL s . The matching is not complete, but a decision was 
made to accept this result since it is conservative. The blue line describes the result from our analysis but 
with setting the lepton detection efficiencies as 100%. It tightens the lepton veto, and slightly loosens the 
bound. This figure is based on a figure produced by the ATLAS collaboration and used in Ref. |81|. 
(ATLAS Experiment ©2012 CERN) 
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Table 5.4: The cross section, the acceptance, and the expected number of events, on the V-GMSB model 
points (A,M mess ,tan/?,My) = (140GeV, 10 6 GeV,20, 1 GeV) and (130GeV, 10 I0 GeV,20, 1 GeV). The 
gluino mass, rrtg, and the mass of the lightest squark among the first and the second generations, rrtg, are 



also displayed. See Fig. 5.6 for the full mass spectrum of the former point. For both points, the "D-tight" SR 
is selected as the expected-to-be-most-promising SR, and thus adopted for the evaluation. For each process 



5000 events are generated. As can be seen in Fig. 5.7 the former model point is not excluded since N < 6.0 
(cf. Table|5.3]>, while the latter is excluded. 



(A, M mess , tan/3) 


(140TeV, 10 6 GeV, 20) 


(130TeV,10 10 GeV,20) 


rrtg 


1 1 16 GeV 


1002 GeV 




1813 GeV 


1887 GeV 


production channel 


99 


99 


99 


99 


99 


99 


cross section cr (fb) 


5.54 


2.75 


0.238 


16.2 


3.31 


0.119 


acceptance A 


0.0758 


0.184 


0.201 


0.0488 


0.155 


0.164 


N = o-xA x5.8fb"' 


2.44 


2.93 


0.28 


4.59 


2.98 


0.11 


N (total) 


5.65 


7.68 



The SRs used in the limit calculation, i.e., the expected-to-be-most-stringent SRs, are shown in Fig. 5. 10 
In the figure the red lines show the exclusion limits under assuming full lepton reconstruction, i.e., setting the 
lepton efficiencies as 100%. Discussion on this topic is performed later. What is important is that for almost 
all points the "D-tight" SR is selected. The "D-tight" SR requires five hard jets and large missing energy. It 
is easy to understand the origin of the four hard jets: the gluino is lightest among the colored superparticles, 
and it generates at least two hard jets in its cascade decay. The origin of the fifth jet can be understood as 
follows. In the channel of pp — > "gq, the hard fifth jet is produced in decays of the squarks into the gluino. In 
the pp — > "gg channel, the largest fraction of events which pass the "D-tight" cut have additional hard jet(s) 
from ISR and FSR (cf. Ref. [ 129 1). The ISR can yield additional hard jet(s) in the pp — > "gq channel as well. 
In the rest of the events the fifth jet comes from decays of W ± and r. Here W ± is generated in the decay of 
the top quark, and r comes from the cascade decay ofpj 1 ^) — > tt — > ttJ^, where x\(Xt) 15 provided from 
decays of colored superparticles. 

One might worry that the hardness of the ISR and FSR jets might not be simulated correctly in Pythia. 
In order to resolve this uncertainty, the analysis for the "D-tight" SR with the MLM-matching scheme 11301 . 
which is implemented in MadGraph, is performed at several model points. The analysis employs a shower 
k, clustering with avoiding double counting between the gluino and squarks 11311 . With this analysis it is 
checked that the results from Pythia agree with those with the MLM-matching. 

The main production channels of SUSY events are pp — > "gg,"gq,"qq. Although the squarks are rel- 
atively heavy compared to the gluino as shown in Fig. 5.6 and Fig. |5.7[ d), they are not decoupled from 
the productions. Since the squark masses are almost independent of tan/?, so are the LHC exclusion lines 
in Figs. 5.7 a)-(c). For illustration, in Table 5.4 we show the cross section of each channel at two model 
points, (A,M mes8 ,tanj8,My) = (140TeV, 10^GeV,20, 1 TeV) and (130TeV, 10 10 GeV,20, 1 TeV). These 



two points are, as shown in Fig. 5.10 analyzed with the "D-tight" SR, and close to the LHC exclusion limit. 



It is found that the channels of pp — > "gg and "gq comparably contribute to the SUSY searches for a low 
messenger scale (i.e., for a relatively light squark), while pp — > "gg dominates for a high messenger scale (a 
heavy squark). 



Appendix — On lepton efficiencies 

To check the effect of the lepton efficiencies, the above analysis is also performed with setting all the lepton 
efficiencies as 100%. Here the leptons are perfectly detected, and it is expected to tighten the lepton veto. 
As is expected, the 100% lepton efficiencies result in slightly looser limits, which are illustrated in 
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Figure 5.10: The same as Fig. 5.7 but here, for the cases where the neutralino is the NLSP, the SRs used 
for the limit calculation are shown instead of the upper limits with uncertainties. In addition, the limits 
obtained with assuming the lepton efficiencies as 100% are shown as red lines to see the effect of the lepton 
efficiencies. 
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Figure 5.11: The same as Fig. 5.8 but here the lepton efficiencies are set to be 100%, which tightens the 
lepton veto and makes event yields decrease by ~ 5%. 



Figs. 5.11 and 5. 10 and also in Fig. 5.9 In Fig. 5.11 the histogram of Fig. 5.8 is reproduced with assuming 
100% lepton detection. The event yields decreases by ~ 5%. In Fig. |5. 10| the exclusion limits with nominal 
cross section is drawn with red lines. It is observed that the bound is imperceptibly looser than the original 
one (the black solid lines). This comparison is also performed in Fig. 5.9 , there the bound with 100% 
efficiencies are drawn with a blue line. 



Now everything what Author wanted to explain is fully shown. Let us move on to the stau NLSP case. 



5.6.3 LHC constraint for the V-GMSB with long-lived stau NLSP 



The LHC constraints for this case are also already drawn in Figs. |5.7| a)-(c); the black solid lines above the 
light blue lines. The left-side regions of the lines are excluded. 

To obtain this limit we just utilize the CMS bound of 223 GeV, Eq. ( |5.37[ ), corresponding to the 95% 
CL exclusion, obtained from searches for the heavy long-lived staus via the direct production |116|. The 
CMS collaboration evaluated the theoretical uncertainties as 3-7% in the cross section calculations due 
to renormalization and factorization scales, a s , and choice of PDFs. In the figure, we assigned theoretical 
uncertainty of 8% for the cross section, which corresponds to 2% uncertainty for the stau mass bound, shown 
by black dashed lines. Consequently, the regions with the gluino mass lighter than 1200 GeV are excluded; 
the bound becomes more severe for a larger tan /3. 
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Figure 5.12: (left) The lower bound, obtained from the LHC SUSY searches, on the gluino mass when the 
NLSP is the neutralino. (right) That on the lighter stau mass for the cases where the NLSP is the stau. The 
yellow (orange) regions show the parameter spaces where the muon g - 2 discrepancy can be explained 
within 2<x- (lcr-) level. 



5.6.4 Discussion 

The LHC constraint can be interpreted as a function of the messenger scale in each category of the NLSP. 
The results are shown in Fig. |5.12| The left and right panels correspond to the cases with the neutralino NLSP 
and the stau NLSP, respectively. Note that, throughout this section, the NLSP is assumed to be long-lived. 
For a given M mess , the muon g - 2 discrepancy can be explained at the lcr (2(f) level when the gluino/stau 
mass is within the orange (yellow) region. The upper ends of these regions represent the upper bounds on 
the gluino/stau mass in order to explain the muon g - 2 at the lcr (2cr) level. 

In the case with the neutralino NLSP (left panel), the upper bound is determined by the requirement 
that the neutralino is lighter than the stau. In fact, for a fixed valu e of the muon g - 2, the gluino mass 



is maximized when tan/? comes on the light blue lines in Fig. 5.7 When M mess is as large as 10 GeV 



the vacuum stability bound can give a more severe bound, but the result is almost unchanged as can be 



noticed from Fig. 5.7 a). This upper bound should be compared with the lower bound on the gluino mass 
from the LHC SUSY search, drawn with the black solid line. Here, the exclusion limit at tan/? = 20 is 
used for illustration, since it does not depend much on tan/?. The black dashed lines show the theoretical 
uncertainties. (See discussion above.) It is found that the whole region where the muon g - 2 discrepancy is 
explained at the lcr level is already excluded by the direct searches for the superparticles at LHC. The region 
with the 2cr explanation is still viable, and expected to be covered with the 14TeV LHC. 



For the stau NLSP case, the upper bound on the stau mass is shown in the right panel of Fig. 5. 12 The 



upper ends of the yellow and orange regions are obtained from the requirement that the stau is the NLSP. 



According to Fig. 5.7 c) and (d), the stau mass is maximized for a fixed muon g — 2 when it is close to the 
neutralino mass. The black solid line is the lower bound on the stau mass from the LHC. It does not include 
the theoretical uncertainties; they are taken into account in the black dashed lines. Although the vacuum 



stability condition can give a tighter bound for M mess ~ 10 9 -10 10 GeV, such parameter regions are already 



constrained by the LHC (see Fig. 5.7 'a)). Consequently, it is found that the region where the muon g - 2 
is explained at lcr-level is fully excluded by the searches for the heavy long-lived charged particles at the 
LHC. It is even expected to become more severe, for instance by analyzing the data obtained at the 8 TeV 
LHC. 
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Section 5.7 Searches for the Vector-like Quark 

Finally, let us discuss searches for the extra vector-like quarks. The LHC experiments have not discovered 
any extra quarks, and thus the V-MSSM is constrained. What is important is that what is constrained is the 
SUSY-invariant mass terms Mq and Mw together with Y', and therefore the constraint is independent on the 
SUSY-breaking scenario. Therefore, the discussion performed in this section is not specific to the V-GMSB 
model, but generic for the V-MSSM scenario. 

The extra vector-like quarks are produced with considerable cross sections at the LHC, a hadron collider. 
Searches for these particles are of great interest and importance because the existence of these particles would 
be a direct evidence of the V-MSSM, and besides, because these particles should be near the TeV-scale as 



we have seen in Sec. 5.5 they are expected to be within the reach of the LHC. 

The vector-like quarks have characteristic decay channels, which are different from the ordinal fourth 
generation quarks. We first check their masses and decay modes, and then current experimental bounds on 
the masses of those particles are discussed. After that we will mention prospects of further searches. 



In the following discussions we use Y' = 1.05, the infrared fixed point value 1 931 , as a reference value. 
The approximation Mq = Mj/<(=: My) is also exploited; one should understand Mq - Mu> when we refer 

to My. 

As the Standard Model parameters relevant to the decay branch, the following are used: 

m f = 173.5GeV, m z = 91.2GeV, v = 174GeV, 02= 0.652, 

m h = 4.78 GeV, m w = 80.4 GeV, m h = 126 GeV, g z = 0.743. 



5.7.1 Masses and decay modes 



The masses of the vector-like quarks are already reviewed in Sec. 5.1.2 When we set Mq - Mjj> =: My, 
the masses are approximately given as 



M 2 Q , + Ml, + Y' 2 vl 



M 2 Q , + M 2 V , + Y' 2 vl 



-M 2 ,M 2 *M V \l±a + 



my - Mq = My, 

where one can see the mass hierarchy 

< my < m t i , 

and the mass splittings among them are characterized by 

Y'v&inp 91 GeV 



2M, 



Q 



My 



Y' 

L05 



(5.42) 
(5.43) 

(5.44) 
(5.45) 



with a large tan (5. As a reference, Fig. 5. 13 is provided, which shows the masses of the vector-like quarks as 
functions of (Mq>,Mu>). 

As we assume that the vector-like quarks are only mixed with the third generation Standard Model quarks 
(63, U3, and D3), the possible decay channels are summarized as 



t' 2 -> b'W, t\h, t\Z, 



b' 



t\W, 



bW, th, tZ, 



(5.46) 



where some of them may be kinematically forbidden if the mass separation is smaller. 

The lightest one, fj , is the most important for the vector-like quark search because of its large production 
cross section. As we assume very tiny mixing between the Standard Model quarks and the vector-like quarks, 
the pair-production pp — > t'^ is the most promising channel, and it leaves characteristic final state particles 
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Figure 5.13: Masses of the vector-like quarks as functions of (Mq>, Mjj>). The mixing term is taken as 
Y'v u = 1.05 x 174 GeV. 



Table 5.5: Benchmark points for the mixing parameters. The shown values as the branching ratios are 
calculated at m, ; = 400GeV (M Q = M w =* 483GeV), T = 1.05, tan/? = 30 and m h = 125GeV; they 
have nontrivial dependence on , but are almost stable under changes of the absolute values of the mixing 
parameters (e's) as long as the mixing parameters are much smaller than 0(1). 





eu ■ 
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<f 


Br(f; -> bW) 


Bi(t[ -> fZ) 


Br(t[ -> th) 


(A) 








1 


1 








(B) 


1 


1 


1 


0.51 


0.44 


0.05 


(C) 


1 








0.48 


0.13 


0.39 
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1 





0.15 


0.19 


0.65 
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1 
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0.01 


0.48 
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Table 5.6: Summary of current vector-like quark searches. The displayed limits are at 95% CL. Several 
obsoleted analyses are also shown with parentheses as a reference. For ?' — > bW search, "2/" denotes the 
di-lepton channel bbWW — > bbllvv, while "/ + /' does the lepton plus jets channel bbWW — > bblvjj. 



Decay channel 


EXP. 


Analyzed data 


Obtained limit 


References 


Br{t[ -> bW) = 1 


ATLAS 


7 TeV, 4.7 ft- 1 (/ + j) 


m f > 656 GeV 

'i 


ii32i rni33n 


CMS 


7TeV,~ 5 ft)- 1 (/+ /) 


»V > 570 GeV 

'i 


111341 (11135]) 


7 TeV, 5.0ft- 1 (20 


m,< > 557 GeV 

'i 


11361 


CDF 


1.96 TeV, 5.6 fb" 1 (21) 


m e > 358 GeV 

'i 




Br(f' t -> qW) = 1 


ATLAS 


7 TeV, 1.04 ft- 1 (20 


m f > 350 GeV 


cm 


CDF 


1.96 TeV, 5.6ft- 1 (/ + j) 


m f > 340 GeV 

'i 


EH 


DO 


1.96 TeV, 5.3 ft" 1 (/ + j) 


m e > 285 GeV 

'i 




Br(f'j -> fZ) = 1 


CMS 


7 TeV, 5.0 ft- 1 


me > 625 GeV 


11401 (111411) 


f'j -> feW, fZ, f/z 


ATLAS 


7 TeV, 4.7 ft-' 


See Fig. 5.17 


mm 



with multi-fe-jets plus leptons depending on the decay pattern. Thus the decay branching ratio of f. is crucial 
for the vector-like quark search. 



The decay widths 1931 are summarized in Appendix 5.B and the branching ratio can be calculated 
straightforwardly. Then one finds that the ratio has nontrivial dependence on the mass me ] , and the ratio 
among the mixing parameters ejj, e^, and eu (93] . Here note that the branching ratio is insensitive to the 
absolute value of the mixings as long as they are tiny, and also that the dependence on tan /J appears only 
with a form eol tan/? as we focus on the cases with a large tan/?. 

In order to proceed the discussion, we pick up several mixing patterns shown in Table |5.5| as benchmark 
points. At the benchmark point (A), t' { exclusively decays into bW, or Br(f'j — > bW) = 1, for any m,j. 
However, except for that point, the decay branching ratio of t\ has nontrivial dependence on the mass of t\ . 



(cf. Figs. 5.14-5.16 ) This is mainly because the t\ 
below "thresholds." Especially, if me < m, + m z - 



tZ and t\ 



th channels are closed if the mass of fj is 



264 GeV, only the t \ — > bW channel is open regardless 
of the mixing parameters; this case corresponds to My < 316 GeV under the Mq = My approximation. 



5.7.2 Current experimental bounds 

Searches for the vector-like quarks can be employed in a similar way to those for the fourth generation 
quarks. The current reports from the LHC experiments are summarized in Table [5~6] together with results 
from the Tevatron collider. If t\ decays exclusively into bW, i.e., Br(?j — > bW) = 1, results from searches 
for the fourth generation up-type quark can be applied, and the tightest bound is m,j > 656 GeV 11321 . The 
ATLAS collaboration also obtained a bound for the case in which Br(fj — > qW) = 1, where q is a generic 
down-type quark: q = (d,s,b), of me > 350 GeV 11381 . For the case where Br(fj -» tZ) = 1, the CMS 
obtained a limit of me > 475 GeV 11411 . 

Those bounds can be applied to the above model points (A)-(E). The result is shown in Fig. 5.14 5.16 
In Fig. |5.14| the 7 TeV LHC bound for the f -> tZ channel reported by the CMS collaboration [ 140 1 is 



displayed as a red solid line. They reported upper bounds on the production cross section in a mass range of 
400 GeV < me i < 650 GeV. The black solid line is the total cross section of pp — > tfj^, and the other (dotted 
and dashed) lines are the "effective" cross section of pp — > fjfj — > (bW)(bW). It is observed that the model 
points (B) and (E) have lower bounds of m,j > 400 GeV and 460 GeV, respectively. Note that the line for 
(A) does not appear in this figure because Br(fj — > tZ) = 0. 

The results for the t' { — > bW channel from the LHC experiments are shown in Figs. 5.15 here one 
result from the ATLAS collaboration |[T32l and two from the CMS llT34l[T36l are displayed (cf. Table [BTB). 
Benchmark points (A), (B), and (C) are constrained as, respectively, m^ > 650 GeV, 420 GeV, and 420 GeV. 
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Figure 5.14: The CMS experimental 95% CL upper limit (red line) on the f'jf'j pair-production cross section 
with 5.0 fhr 1 data, assuming f. exclusively decays via t' x — > tZ 1 140 1. Their limit is obtained for 400 GeV < 
m,j < 650 GeV. The black solid line is the NNLO total cross section of t'f production, calculated with 
HATHOR II 1421 . Considering the branching ratio Br(fj — > fZ), this limit can be applied to the vector-like 
quark with generic decay branch. We show the corresponding f'jfj cross section with the branching effect at 
the benchmark points (A)-(E) as dashed and dotted lines. Note that the line corresponding to the point (A) 
is not shown since Br(f j — > tZ) = at the point. 
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Figure 5.15: The same as Fig. 5.14 but here f. is assumed to decay exclusively via f. — > bW channel. The 
red line denotes the upper limit from the ATLAS B132I . while the blue solid and the blue dashed lines are 
from the CMS II 1 341 1 1 361 . The line for (A) overlaps the total cross section line since Br(f'j — > bW) = 1. 
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Figure 5.16: The same as Fig. 5.15 but from the CDF experiment at the Tevatron B137L t' x is assumed to 
decay via fj — > bW channel, and the cross section corresponds to 1.96 GeV pp collision. The red solid line 
denotes the upper bounds on the cross section. The other lines are the same as Figs. 5.14 -5.15 but the cross 



sections are for 1.96 GeV pp collision. Note again that the line for (A) overlaps the total cross section line. 



The point (D) is not constrained by the above results; one reason is that the LHC results are reported only 
for m,j > 350 GeV. To give constraint for this benchmark point the analysis by the CDF collaboration 1 1 37 1 
can be utilized. Their result obtained at the Tevatron collider with Eqm = 1 -96 Te V pp collision is displayed 



in Fig. 5. 16 and yields a limit of m,j > 300 GeV for the benchmark point (D). 

In this analysis the production cross sections are calculated with HATHOR 1 . 3 at the NNLO level. The 
CTlQnnlo PDFs lfl43l are used. 



An interesting analysis was reported by the ATLAS collaboration [ 132 1 . They employed an inclusive 
search for the vector-like quark t\ in the events with at least three jets, at least one of which should be 
tagged as a b-jet, and at least one lepton (e or p), and reported the excluded regions on the parameter space 
of [Br(f' — > bW),Br(t' — > f/i)]-plane with the mass m,< fixed. The result is shown in Fig. 



5.17 



with a 

modification that our benchmark points (A)-(E) are plotted instead of theirs. Interpreting the figure, we 
obtain tighter bounds for the benchmark point (B) as > 450 GeV and (C) as > 500 GeV. 

5.7.3 Prospects of further searches 



In order is discussion on future prospects to conclude this section. Fig. 5.18 shows the production cross 
sections of pp(pp) — > t' y 7 y as functions of the t' { quark mass, together with an extra axis of My(= Mq = Mu>). 
The cross sections are calculated with HATHOR 1 . 3 014211 and the CTlQnnlo PDFs 1 1431 at the NNLO level. 
The cross section doubles at the 8 TeV LHC, and is ten times larger at the 14TeV, than that at Eqm = 7TeV. 

As we have just seen, the inclusive search for pp — > t't' — > (bW,th,tZ) + (bW,th,tZ) by the ATLAS 
collaboration 11 1 3211 is of great importance and interest. Especially, because the pair-production cross section 



Chapter5 The MSSM with Vector-like Matters 



77 



A 



0.8 



0.6 



CC 0.4 
00 

0.2 




: V 

- % 

(D) . 


m t ,= 


400 GeV 








r 






i , , , i , 


(B) 


**X(A) 


0.2 


0.4 0.6 


0.8 1 



ATLAS 

7 TeV, 4.7 fb" 



95% CL expected exclusion 



95% CL observed exclusion 




: V 

: (D) ° 


m t - = 500 GeV 


-_ ★ 




1 (E) 




r* 


(C) i \ 


i , , , i . . , 


(B) ' 



0.2 0.4 0.6 0.1 




0.2 0.4 0.6 0.8 

BR(t' -> Wb) 



Figure 5.17: The result of the inclusive search for the vector- like quark fj reported by the ATLAS collabo- 



ration 111321 . In this analysis the decay branch of t' { is not specified, but the assumption that the vector- like 
quarks only mix with the third generation Standard Model quarks is employed. Therefore, this search gives 
an inclusive study for the three decay patterns: fj — > bW, tZ, th. The result is shown in the plane of the 
decay branch; here one should note that Br(f'j — > fZ), which is not displayed explicitly in this figure, is equal 
to 1 - Br(f'j — > bW) - Br(f'j — > th), and thus implicitly shown. This report gives tighter bounds for our 
benchmark points (B) and (C) than the analyses focusing on a specific decay branch, i.e., of Figs. |5.14| and 
|5.15| Original version of this figure is produced by the ATLAS collaboration and used in Ref. 1 132 1. 
(ATLAS Experiment ©2012 CERN) 
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Figure 5.18: The production cross sections of pp(pp) — * tft as functions of the t' { quark mass, calculated 
with HATHOR 1 . 3 lfl42l at the NNLO level. 



drastically improves at the 14TeV LHC, we can expect that the t' x quark with ~ 1-1.2 TeV will be in our 
experimental reach. As the V-GMSB model requires My < 1 .2 TeV to explain the muon g - 2 anomaly, the 
14 TeV LHC would be the court for the V-GMSB model; that is, the fate of the V-GMSB scenario, and thus 
the GMSB framework, with the explanation of the muon g - 2 problem will be determined there. 

Here it should be emphasized that the ATLAS analysis utilizes the event with at least one fo-jet. As the 
fj decay generally yields 1-3 b-quarks, its pair-production would be searched more efficiently with multiple 
Z?-jets requirement ll94l I144L In Ref. |144| the prospects of this strategy are discussed with Monte Carlo 
fast simulations, and it is found that the requirement of 1 lepton+ > 3 b-jets would yield complementary 



exclusion regions on the parameter spaces of Fig. 5.17 This analysis would also be an important key for 
searches for the heavier vector-like quarks, V and t' 2 . 

On the experimental side, we now realize that ^-tagging, discussed briefly in Chapter [3] is acquiring 
importance from the fact that the Higgs boson is expected to decay mainly into fe-quarks. Improvements 
on the efficiency and the mis-tag rate, and deep understanding of the algorithms, are now much awaited for 
searches of new physics beyond the Standard Model. 



Section 5.8 Summary 



Now the main chapter of this dissertation is finishing. Let us summarize what we learned. 

We in this section considered the V-MSSM, an extension of the MSSM with vector-like matters. It has 
three virtues. First, it respects the daydream of the SU(5)-GUTs. Also, the Higgs boson mass of 126 GeV is 
realized with help of the vector-like quarks without exploiting the heavy stops of ~ 1-10 TeV. Then squarks 
need not to be so heavy, and thus sleptons either. It allows us to explain the muon g - 2 discrepancy with the 
SUSY contributions, even under the GMSB scenario. 



These features are summarized in Fig. 5.7 for the V-GMSB model. Even after imposing the LHC limits, 
the yellow regions, where the muon g - 2 anomaly is explained within 2<x-level, survive in all the three plots. 
With the LHC limits, the gluino mass is constrained as nig > 1.0 TeV under the V-GMSB model. 

Another important point in this figure is that the mass parameter My is also constrained for successful 
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explanation of the muon g - 2 anomaly. For (c)M mess = 10 6 GeV case, the mass parameter My should be 
~ 800-1200 GeV; for (a) and (b), it should be lighter as ~ 700-1000 GeV. Therefore, the vector-like quark 
should be lighter than < 1 . 1 Te V, and this is obviously within the reach of the 14 TeV LHC. 

As for the vector-like quark searches, the one by the ATLAS collaboration Ml 321 is very interesting and 



promising (cf. Fig. 5.17 1. This analysis is based on the requirement of at least one fe-jet, but since more 
/7-quarks are expected from the decays of t[, much more parameter regions is expected to be excluded with 
multi-b-jets analysis [ 144 1. The importance of ^-tagging should be, therefore, emphasized now. 

Several future works are remained just for the vector-like quark searches. First, in this dissertation 
we simply restrict the model to have no mixings between the vector-like matters and the Standard Model 
fermions in the first and second generations, and to have small mixings between the vector-like matters and 
the third generation. Actually the Z?-jet signature is provided by this assumption; if the vector-like quarks 
mix with the lighter quarks, other strategies are mandatory for vector-like quark searches. In this context, 
quantitative evaluation of current experimental bounds on the mixing parameters is awaited. 

Analyses for the cases with a short-lived NLSP are, although we briefly discussed, left as future works. 



From a theoretical viewpoint, the V-MSSM scenario has several interesting features. One of them that 
the large /i-term. This feature actually worsens the argument of the little hierarchy problem. Therefore, it is 
still difficult to settle the collision between the little hierarchy and the 126 GeV Higgs mass, and hence, with 
examining this argument, another window towards fundamental theories might arise, as we have seen that 
the muon g - 2 anomaly bore the V-GMSB scenario. 

The mechanism to suppress the coupling Y" is also worth investigated, for which a model with Peccei- 
Quinn symmetry was proposed II 1451 . 

Cosmological discussion, which is not performed in this dissertation at all, is also left as future works, 
such as the effect of the vector-like matters to the history of our Universe, and consequences of the large 
strong coupling gj, in the early universe (cf. Fig. |5. 1 [ >. 
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Appendix 5.A Renormalization Group Equations for the V-MSSM 

The two-loop level renormalization group equations of the V-MSSM are shown here, which are used in our 
numerical evaluation in this chapter. As emphasized in Ref. [93], the two-loop level effect is significant 
especially for the running of the gaugino masses. 

For completeness all the yS-functions up to two-loop level of the renormalization group equations are 
included. The /^-functions of the equations are defined as 

jt^ = i^/ 3 m + r^ (2) [x] ■ (5 - 47) 

dlogg I6n 2 (16tt 2 ) 2 

where X are one of the parameters and Q is the renormalization scale. Shown are not only the /3-functions 
for the extra parameters but also corrections to the MSSM /^-functions ySMSSM, i-e. 

^mssm=<ssm + ^ (,) - (5-48) 
Note that the definition of the parameters is shown in Eqs. (j5.4|i and (|5.7|), and the MSSM y6-functions are 



listed in Appendix 4. A 



The /J-functions are calculated with Susyno 1.1 |87|. 



5.A.1 Restriction and notation 



The following assumptions are employed, as is done for the MSSM /^-functions (in Appendix 4. A I. 

• The R -parity is conserved. 

• The scalar soft mass terms m 2 , are diagonal, 

• For the A-terms a x and the Yukawa coupling Y x , all the components but the (3, 3) are neglected. 

• The gaugino masses M a are real. 

The DR' scheme ll65l is chosen as the renormalization scheme. 



The following variables are used in the expressions of the /^-functions. 

X, := 2a) + 2T 2 (m^ + (m 2 Q h 3 + (m 2 D h 2 ) 
2a 2 b + 2Y 2 h (m 2 Hd + {m 2 Q ) 33 + (m 2 ) 33 ) 
2a 2 + 2Y 2 T (m 2 Hi + (m 2 ) 33 + (m|) 33 ) 
2a' 2 + 2T' 2 {m 2 Hn + m 2 Q , + m 2 ^) 



X T 
X' 
X" 



a (t,b,T) 



2a" 2 + 



2F" 2 ( 



m % + m % + m l 



>) 



'(')„'(') 



= Y' w a 



S := m 2 Hu - m 2 Hi + ^ [(»ig)« - 2(m 2 - J ) ii + (m^u - (m 2 L ) u + (m|) s ] 



S' 

S (2) 



i=\ 



(iii 2 q, - 2m 2 y, + 'w 2 ?,) — ( m2 Q' ~ Zm^j, + m 2 E ^ 

-if [3m 2 Hu + (m 2 ) 33 - Mm 2 ^) + Y 2 b (3m^ - (m 2 e ) 33 - 2(m 2 ) 33 ) + T 2 {m\ & + (m 2 ) 33 - 2(m|) 33 ) 



2 2 
m H a ~ m H d 



- YjP&ii + (-U 2 + ]-g\ + L 2 ] 2(m 2 e )«« 



16 
T5 ! 



15 



30' 

3 



i=1 



6 , 
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o-i := p0i 



0- 2 



0- to t;l 
O" tot;2 
CT tot;3 



02 



3m^ + 3m^ d + J] ((m 2 ,),, + 8(m^) s + 2(m 2 ),,- + 3(m£) rt + 6(m|) ff ) 
i=i 

3 



cr 3 := 



;(2) 



i=i 

5 (2) - 7' 2 (3ml, + m 2 , - 4m 2 ,) + 7" 2 (3m 2 fj + m 2 , - 4m 2 ,,) 



1 2 3 2 8 
30 9 ' 1 + 2 91 + 3 



g 2 J (m 2 ,, - m|,) - |^0 2 + yd/ 2 J (m 2 ,, - m 2 ,,) + ^g 2 (m 2 , - m%) 



0-1 + 



\g\ {m 2 Q , + m|, + %m 2 , + im\j, + 6m 2 , + 6m|,) 



cr 2 + g\ (3m 2 Q , + 3m|,) 

0-3 + g\ (2m 2 Q , + 2m|, + m 2 , + m^,) 



5.A.2 One-loop level ^-functions 

Corrections to the evolution of the MSSM parameters 





= 3g? 


(5.49) 


A/3 (1) [Mi] 


= 6g 2 M, 


(5.50) 


Aj6 (1) [7] 


= 3Y' 2 Y, 


(5.51) 


A/? (1) [Y b ] 


= 3Y" 2 Y h 


(5.52) 


A/3 m [Y T ] 


= 3Y" 2 Y T 


(5.53) 


A/3 (1) [a,] 


= 37' (7'a, + 2a' 7) 


(5.54) 


Ay6 (1) [fl 6 ] 


= 37" (7"a fo + 2a" 7,) 


(5.55) 


A/3 (1) M 


= 37" (7"a T + 2a"7 T ) 


(5.56) 


AyS (1) Lu] 


= 3//(7' 2 + 7" 2 ) 


(5.57) 


A^ (1) [ft] 


= 6//(a' + a") + 3ft(7' 2 + 7" 2 ) 


(5.58) 


¥ l) KJ 




(5.59) 


A/? (1) KJ 


= 3X'- 3 -g\S' 


(5.60) 


A/3 (1) [(m 2 G ),] 


- \v\s' 


(5.61) 


A/3 (1) [(m 2 & ),] 


= >' 


(5.62) 


A/? (1) [(m£) ff ] 


- f rfs' 


(5.63) 


A/? (1) [(m 2 ),] 




(5.64) 


A/3 (1) [(m|) ff ] 




(5.65) 
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Evolution of the extra parameters 

fF> [Y>] = [3Yf + 6F' 2 - f 5 g\ - 3g 2 2 - fg 2 } r (5 . 66) 

^ [Y"] = {sY 2 b + F 2 + 6F" 2 - gg 2 - 3gj - fgjj Y" (5.67) 

^ [Mq>] = (y' 2 + Y" 2 - lg 2 - 3g 2 - ^ 2 ) M G , (5.68) 

ffi» [M v ,] = {lY' 2 + 2F" 2 - ^ 2 - y^M^ (5.69) 

0- 1) W E ,] = ~g\M E , (5.70) 

[a'] = [a t + ^ 2 M, + 6;? 2 M 2 + y g 2 M 3 ) F' + ^3Y? + 18F' 2 - ^ 2 - 3g 2 2 - y g 2 J a' (5.71) 

/ 26 32 \ 

[a"] = 674 + 2a T + — 2 Mj + 6g 2 2 M 2 + — g 2 3 M 3 \ Y" 

+ (3y 2 + y 2 + i 8 y" 2 -^-3, 2 -f, 2 ) fl " 

l) [b & ] = {la' + 2a" + ^g\Mi + 6g 2 2 M 2 + y0 2 M 3 j A% 



(5.72) 



(5.73) 



ft» [b„.] = \*3 + 4a" + ^ 2 M, + |g 2 M 3 ) My. + (-j^ 2 - fgj + 2F' 2 + 2F" 2 ) V (5.74) 

(5.75) 

24 12 

y6 (1) [b P ] = -^g\M,M E> - —g\b E , (5.76) 

P m K'l = _ ^i M i " 6 ^ M 2 " f 03 M 3 + +S') + X' (5.77) 

H>] = "1^ M ? " f & M \ ~ + 5 ') + 2X ' ^ 

f?> [ OT |,j = ~g 2 M 2 + ^g 2 (S + S') (5.79) 

P m HI'] = "^i M i " 6g ^ " f^ M 3 - +S') + X' (5.80) 

K'] = -§ 9 i M i " y ? 3 M 3 2 + \&S +S') + 2X' (5.81) 

[ml] = -y 2 M 2 - ^ 2 (S + S') (5.82) 

5.A.3 Two-loop level ^-functions 
Corrections to the evolution of the MSSM parameters 

A/? (2) [ gi ] = (f + |^ + f ^ - |r 2 - |y" 2 )^ (5.83) 

A/? (2) [02] = (^ + 21g 2 + 16^1 - 6F' 2 - 6F" 2 J g\ (5.84) 



A/? (2) [03] 



^g 2 + 6g 2 + 34g 2 - 4F' 2 - 4F" 2 J g\ (5.85) 
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A/? (2) [M!] = 

A/3 (2) [M 2 ] = 
A/? (2) [M 3 ] = 

A^ (2) [FJ = 
A/3 (2) [y 6 ] = 
Ay6 (2) [F T ] = 
AyS (2) [a,] 



A/3 (2) [fc] 



^ (7 - F' 2 M! + a" - Y" 2 M X ) + ^-g\M x + -g\ (M, + M 2 ) + ^gj (M x + M 3 ) 



01 



(5.86) 



12 (a' - F' 2 M 2 + a" - Y" 2 M 2 ) + -g 2 (AT, + M 2 ) + 84g 2 M 2 + 32g 2 (M 2 + M 3 ) 



<7 2 (5.87) 



8 (a' - F' 2 M 3 + a" - y" 2 M 3 ) + ^-g\ (Afj + M 3 ) + 12^ (M 2 + M 3 ) + 1360 3 A/ 3 



y 0? + 9g 4 + ^4 + ^ + ^ _ 9yr2 _ 9y2 j Y>2 _ 3y2y „ 2 

^1 + 9 g \ + 16^ - 3 y 2 r 2 + (^ 2 + 16 2 - 9 y 2 - 9 y" 2 ) y" 2 



^ .4 , n.4 , |^ 2 



y 0i + 9^ + 1 5 + 16 9 2 - 9y T 2 - 9y" 2 ) y" 2 



Yr 



y <tf (a, - 4yM0 + 9^ (a, - 4yM 2 ) + 16g 4 ( fl , - 4yM 3 ) 

+ -0 2 y' (y'a, + 2a'y - 2M{Y'Y t ) + \6g\Y' (Y'a, + 2a'Y, - 2M 3 Y'Y t ) 



(5.88) 
(5.89) 

(5.90) 

(5.91) 

(5.92) 



+ y (-27y'a ( y 2 - 9y'V - 36a'y' 2 y - i8a'y 3 ) 

+ Y b Y" (-3Y"a,Y h - 6Y"Y t a b - 6a"Y t Y b ) 
A/? (2) [a b ] = n -g\ (a b - AY b M x ) + 9g 4 2 (a h - AY b M 2 ) + 16a 4 (a b - AY b M 3 ) 

+ A -g\Y" (Y"a b + 2a"Y b - 2M\ Y"Y b ) + \6g 2 3 Y" {Y"a b + 2a" Y b - 2M 3 Y"Y b ) (5 m) 

+ Y t Y' (-3Y'Y t a b - 6Y'a,Y b - 6a'Y,Y b ) 

+ y" (-27y"a fe y 2 - 9Y"\ h - 36a"Y" 2 Y h - 18a"y 3 ) 

A/3 (2) [Or] = y<? 4 (fh - 4Y T Mi) + 9g\ (a T - 4Y T M 2 ) 
4 

+ -g\Y" (Y"a T + 2a"Y T - 2M X Y"Y T ) + \6g\Y" (Y"a T + 2a"Y T - 2M 3 Y"Y T ) 



+ y" (-27y"a T y T 2 - 9y" 3 a T - 36a"y" 2 y T - i8a"y 3 ) 
\g\ + 9g* + \g\ (r 2 + y" 2 ) + i6a 2 (y' 2 + y" 2 ) - 9 (y 4 + y" 4 ) 



36 



g\M\ - 36g A 2 M 2 - 36a' Y' 2 - 36a 'Y 



+ -g\ {a' + a" - (y' 2 + Y" 2 ) M x ) + 32g\ (a' + a" - (y' 2 + Y" 2 ) M 3 ) 

%g\ + %\ + \g\ {y' 2 + y" 2 ) + i6g 2 (r 2 + y" 2 ) - 9 (y 4 + y" 4 ) b 



(5.94) 



(5.95) 



H (5.96) 



54 

—g\M 2 + 54g 4 2 M 2 +^gfS' + ^gfa^ + 3g 2 2 a J 2 



+ -g\ (X' - \M{3 + 4M 2 y' 2 ) + 16^ (x' - \M{a + 4M 2 y' 2 ) - 36^' 2 - 18y' 2 X' 



(5.97) 
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+ -g\ (X" - 4M,a" + AM\Y" 2 ) + I6gj (x" - 4M 3 a" + 4M 2 F" 2 ) - 36a' 



A/?< 2 ' [(m 2 e ),] = 6 -g\M\ + 54g$M 2 2 + 96g 4 3 Mj + g\ ||s' + + 3g 2 2 (/ 2 + y ^ 
+ | - 12a'a, - 3Y 2 X' - \2a"a b - 3Y 2 X" - 3Y" 2 X h - 3Y' 2 X, 



for i = 3 



A/^ (2) [(«&)«] = + 96^Mf - + ^ + ^ 3 



16 



16 



15 1 



+ (( - 24a'a ( - 6F 2 X' - 6F' 2 X, 



for i = 3 



^ (2) [(«£)»] = -r9i M i + 96^M 2 + ^5' + + ^ g 2 ^ 



16 



5' 71 i 3 ' 1 3 
+ (( - 24a" a b - 6Y 2 h X" - 6Y" 2 X h 



¥ 2) [(mi) fi ] 



54 



-<? 4 M 2 + 54^M 2 - ^S' + ^cr'j + 3g 2 cr' 2 



+ - 12a"fl T - 37?X" - 3F" 2 X T 



for i = 3 



for i = 3 



A/?® [(«|) fi ] = ™g\M\ + |g 2 5' + + | _ 24a"« T - 6F 2 X" - 6Y" 2 X T jj 



18F" 2 X" 
(5.98) 

(5.99) 



(5.100) 



(5.101) 



(5.102) 



(5.103) 



for i = 3 



Evolution of the extra parameters 



/3 (2) [Y'] = 



3913 4 33 4 128 4 2 2 136 2 2 „ 2 2 

+ y 02 + -g-03 + #102 + -^-VlVl + 8 0203 
+ \t Y ' + 5 9 ' + 692i " 2 + 16 ^ (F ' 2 + 7 ' 2) ~ 3Y ? Y *> ~ 9Y > ~ 9Y ' Y ' 2 ~ 22Y " 



r 



02) p"] 



3913 4 33 4 128 4 2 2 136 2 2 2 2 



+ 0? + ^+ 6 -Y" 2) j + 6g\Y" 2 + 160 2 (7 2 + F" 2 ) - 91* - 3F T 4 

- 3Y 2 Y\ + F" 2 (-9^ - 3F 2 - 22F" 2 ) 



^[M Q> ] = 



289 4 33 



450 



0i 



"02 



128 



16 



03 + 50102 + ^0103 + 16 0203 



+ ^g\(Y' 2 + Y" 2 ) - 3Y' 2 Y 2 - 3Y" 2 Y 2 h - Y" 2 Y 2 - 5Y' 4 - 5Y" 4 



(5.104) 



(5.105) 



(5.106) 
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y6 (2) [M v .] 
ff* [a'] 



fP [a"] = 



2432 
~225 



+ ^3 + ^-M - U {r 2 + y" 2 ) + e g j (r 2 + y" 2 ) 



0i 



45 



6i" 2 if - 6r ,2 y 2 - 2y" 2 y 2 - sr 4 - sr 



7/4 



Mr 



648 4w 

25^"* 

3Q13 33 128 

— (a' - 47'MO + - g\ {a' - AY'M 2 ) + — g 4 (a' - 4F'M 3 ) 



(5.107) 



(5.108) 



+ %g\g\ (a - 2Y'M 2 - 2Y'M 3 ) + g\g\ (a - 2Y'M 1 - 2Y'M 2 ) 

136 to,, , , . i / 8 , 4,9 8 , i 18 , , 12 , 

+ -tffal ( fl " 2Y M i - 2y M 3) + 9? I 5 Y 'a, + -a'F 2 - -M.Y'Y 2 + —a'Y' - —M{Y 

+ g\ (lSa'y' - 12M 2 F' 3 ) + g\ (?>2Y'a t + \6a'Y 2 - 32M 3 Y'Y? + A&a'Y' - 32M 3 F' 3 ) 
- 6Y'a,Y? } - 6Y'Y 2 a b - 3a'Y 2 Y 2 h - lSY ,y a t - 36Y%Y? - 275 'Y'Y? - 9a' 'Y 3 - llOa'F' 4 



(5.109) 



3913 



33 



128 



m g\ {a" - AY" M\) + -g\ {a" - 4Y"M 2 ) + —g\ {a" - AY" M 3 ) 
+ %g\g\ (a" - 2Y"M 2 - 2Y"M 3 ) + g\g\ (a" - 2Y"M V - 2Y"M 2 ) 
+ ^g\gW ~2Y"M 1 -2Y"M 3 ) 

+ 9\ (I (Y 2 b ~ 3Y 2 - 3Y" 2 ) F"M, - 4 -Y"a b Y b + l j-Y"a T Y T + I^Y" - \y\ + ^F 2 ) a") 

+ g\ (l8F"a" - 12M 2 F" 3 ) + \<zg\ (3Y"a" + 2Y"a b Y h + a"Y 2 h - 2M 3 Y"Y 2 h - 2M 3 F" 3 ) 

- 27F"a"y 2 - 9F"a"y 2 - 110y"V - 6Y"a t Y,Y 2 h - 6Y"Y?a b Y b - 3a"Y 2 Y 2 b 

- \%Y" 3 a h Y b - 36Y"a b Y 3 h - 9a"Y A h - 6Y" 3 a T Y T - \2Y"a T Y 3 T - 3a"Y A T 



(5.110) 



578 512 2 32 

— g\M x - 66g 4 2 M 2 - —g$M 3 - -g\g\ {M\ + M 2 ) - —g\g\ (M\ + M 3 ) 



- 32g 2 2 g 2 3 (M 2 + M 3 ) + -g\ (a' + a" - M X Y' 2 - MiY" 2 ) 

- 6Y' 2 a, - 6a'Y 2 - 6Y" 2 a b - 65"Y 2 - 2Y" 2 a T - 2a"Y 2 - 20a'Y' 2 - 20a" Y" 2 



M Q (5.111) 



289 



33 



128 



4500? + y 02 + '-j- 9% + 5 g\g\ + ^gt&i + ^g L 2 g\ 

+ \g\\Y' 2 + Y" 2 ) - 3Y' 2 Y 2 - 3Y" 2 Y 2 h - Y" 2 Y 2 - 5Y' 4 - 5F" 4 



16 



,2„2 



,2 2 
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P i2) [b v >] 



0- 2) \b E ,] 



512 512 
g\M x - —g\M 3 - 4g-0i03 (Mi + 



+ 



9728 
225 

- ^g 2 (a' + a" - M\ Y' 2 - M X Y" 2 ) + Y2g\ (a' + a" - M 2 Y' 2 - M 2 Y" 2 ) 

- 32Y' 2 a' - 32Y" 2 a" - 12F' 2 a ( - \2Y" 2 a h - 4Y" 2 a T - US Y 2 - \2a"Y 2 b - 4a" Y 2 

+ l ~frt + ^ - \A ( y ' 2 + Y " 2 ) + ^ ( Y ' 2 + Y " 2 ) 

- 6Y' 2 Y 2 - 6Y" 2 Y 2 h - 2Y" 2 Y 2 - 8F' 4 - 8F 



-y/4 



2592 



648 4 

-25^ E '- 25 



g\M x M E , 



™g\M\ + %lg\M\ + ™g\M\ + \g\g\ (m 2 + M\ + M X M 2 ) 



32 
45 



g\g\ [m\ + M\ + M X M 3 ) + 32g\g\ [m\ + M\ + M 2 M 3 ) 



+ g^tot + Y501°tot;l + 3^2°"tot;2 + yg30" tot ; 3 

+ ^g 2 [X' - 4M{a + 4M 2 Y' 2 ) - ITSa, - 3Y 2 X' - 3Y' 2 X, - 20a' 2 - 10Y' 2 X' 
0* [m 2 d ] = ^g\M 2 + + ^,fo£ (m? + M 2 + M 1 M 3 ) 

- ^^tot + Jf^^tot;! + y^3^tot;3 " ^ " 4M^ + 4M 2 Y' 2 ) 

+ 6g\ [x' - 4M 2 a' + 4M 2 2 Y' 2 ) - 24Sa t - 6Y 2 X' - 6Y' 2 X t - 32a' 2 - 16Y' 2 X' 
^gtM 2 l + l j g l(s^ + ( r tot ,) 



16 



28 %X + %lg\M\ + ^-g\Ml + \g\g\ [m\ + M\ + M X M 2 ) 



75 



32 

' 45 
2 

' 5 

■ 4a"(3a h + a T ) - {3Y 2 h + Y 2 )X" - 3Y" 2 {X b + X T ) - 20a" 2 - 10Y" 2 X" 



g\g\ [M\ +M\ + MiM 3 ) + 32g\g\ [m\ + M\ + M 2 M 3 ) 
9i S m + A^°"tot;i + 3g\cr Xol . 2 + \-g\o- m $ + A -g\ (x" - 4M x a" + 4M 2 Y" 2 ) 



ff) = A ^-g\M\ + ™g\M\ + ^g\g\ {m\ + M\ + M X M 3 ) 

+ \<&S™ + ^<rtot;, + y g\<T^ - \g\ {X" - 4M x a" + 4M 2 Y" 2 ) 
+ 6g 2 2 (X" - 4M 2 a" + 4M 2 Y" 2 ) 

- 8a" (3a b + a T ) - 2(3Y 2 h + Y 2 )X" - 2Y" 2 (3X h + X T ) - 32a" 2 - 16Y" 2 X" 



15 
160 



512 



M„. 



(5.112) 
(5.113) 



(5.114) 



(5.115) 



(5.116) 



(5.117) 



(5.118) 



(5.119) 
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Appendix 5.B Decay Rates of t' 



Here are the decay rates of the lightest vector-like quark in the V-MSSM, f' p summarized, which are cited 
from Ref. [93 1. The relevant mixing parameters and the mass matrix are, as already referred, 



W = -Y'H U Q'U' + MqQ'Q + M V .U'U' - e u H a Q 3 U' - e' u H u Q'U 3 + e d H d Q'D 3 ; 

(Q'u) 



-£?(Q' U V t L )M, 



U' 



M t := 



'M Q , Y'v u e' u v^ 

My 
, £ u v u Y,v UJ 



Mi, : = 



-M Q e d u d 
Y b v d 



(5.120) 
(5.121) 

(5.122) 



Utilizing the singular value decomposition method, we can diagonalize these two matrices with unitary 
matrices R,R,L',R' as 



L*M,R 



m, ' 
m t > 
0' m. 



L'*M b R^ 



nib 
my 



(5.123) 



The gauge interactions are extracted to be 



+ 4 

n + Z 



lV2 V2 " 

ei (e^e;, + u*w + fa - q'Jq:, - u r w - r RfR ) - he\ (Q'jQ' d + b[b L - - j^h) 



2f|e| (c/'t^_^ R _c;'tc7')-^ R 



3 v / 3 

+ (-T^fRfL - Y b H° d b R b L - Y'H^Q' U U' - e u H° u t L U' - £ u H°uQ',M ~ edH°Q' u b R + H.c.) 

= ^ (L[\L 2l + L\\L 2i ) MTb^ + ^R\\R 2 i ■ W~b% 
V2 V2 



+ ^(^21^11 +^23^13)^1'?- ^-R* u Ri\ -Z^t\ 



2 v Zi ZJ ' 1 2 

- (T,L23^i3 + Y'L 2 iRi2 + e u L 23 R\ 2 + e' u L 2 \R\ 3 ) H®tt\ 

- (Y,L n R 23 + Y'L\\R 22 + e u L 13 R 22 + e u L n R 23 ) H^t[t + H.c. 

Therefore, we obtain the decay rates as, defining the following coupling: 



(5.124) 
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4 



02 



V2 



^'1*1^21, 



r- Wl^21 + ^12^23) 
V2 

y ( L 2i L n +£23^13)' 

cos a , , 

— i=r (Y,L 23 L\ 3 + Y L 2 iR x2 + e u L 23 Ri 2 + e L 2 \R\ 3 ) 

V2 
cos a 



V2 



(Y,L 13 R 23 + Y'L\\R 22 + e u L\ 3 R 22 + e' u L u R 23 ) , 



(5.125) 

(5.126) 
(5.127) 

(5.128) 
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rOi - w+b) = ^ (i - .,) 2 (2 + tf) (|^ t f + |,- t f) . 



(5.129) 



m,< 



\+r,-2r z + 



r z 



: )(l4t| 2 + |4,f)+i2v^Re(4 t 4,) 



r {t\ ^ht)=^ ylM\,r h ,r t ) [(1 + r, - r h ) ^ + |^ ( | ') + 4 V^Re (j^,) 



where 



:= x 2 + j/ 2 + z 2 - 2xi/ - 2j/z - 2zx, r x = — 

mi 



(5.130) 
(5.131) 

(5.132) 
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Chapter 6 

Coda 

In this coda, the theme of this dissertation is repeated. 

* * * 

On 4th July 2012, the ATLAS and the CMS collaborations claimed that they respectively observed a new 
boson with a mass approximately 126 GeV in the search for the Standard Model Higgs boson |T||2]], which 
completes the Standard Model J3HU. This model has the electroweak symmetry breaking as its heart, which 
is governed by the Higgs mechanism, and explains almost all of Nature. 

However, we already know that the Standard Model is not the ultimate theory. We have the Dark Matter 
problem, the fine-tuning problem, and the muon g - 2 discrepancy. Also the mechanism which generated 
current baryon asymmetry of our Universe is still unknown. Moreover, we have four forces, not one. To 
achieve the grand unification, the slight mismatch on gauge coupling unification should be resolved. These 
topics were discussed in Chapter|2] 

Now we have the LHC experiments as a powerful tool to investigate the next theory beyond the Standard 
Model. In this dissertation the ATLAS detector BT1I421 . a general-purpose detector for the LHC, is examined 
in Chapter|3] There we saw that the discovery of the Higgs boson was realized with help from the increases 
of the energy and the collision rate in the 2012 run, and also the importance of /?-tagging was emphasized. 

Then we reviewed the SUSY, a silver bullet for the problems in the Standard Model. It solves the 
hierarchy problem in a miraculous manner, and as we saw in Chapter [4] the muon g - 2 problem and 
the slight mismatch on gauge coupling unification can be solved with the SUSY. Moreover, a promising 
candidate for the Dark Matter is provided. 

However, we have not detected any hints of the SUSY at the LHC. Now the colored superparticles are 
highly constrained as m-g > 900 GeV and nig > 1400 GeV, which indicates that the SUSY is heavier than 
expected. We saw that this indication is supported also by the Higgs boson having a mass of 126 GeV. 
Under the no-mixing scenario, the stop mass is required to be (9(10) TeV. Even with the maximal-mixing of 
X, ~ + 'sfcmj, it should be ~ 1-2 TeV. There the argument of the little hierarchy is not fulfilled. 

Two possibilities were introduced there. One was the heavy-colored light-non-colored scenario, which 
can be investigated with searches focusing on pair-production of charginos, neutralinos, and sleptons via 
electroweak interactions. Such searches do not need higher energy, but require a higher luminosity; the 
data expected in the 13-14TeV run, corresponding to (9(100) fb _1 , would help this direction, and also the 
HL-LHC is of great importance (cf. Sec. |3.3) >. 

However, this scenario is not preferred from a theoretical viewpoint because it does not fully respect the 
SU(5)-GUTs. As a possibility which respects the SU(5)-GUTs daydream, there the V-MSSM scenario, an 
extension of the MSSM with vector-like matters, was introduced, and was investigated in Chapter|5]as the 
main dish of this dissertation. 

We saw that in this model the Higgs mass can be raised by the extra vector-like quarks. It allows us to 
explain the muon g - 2 anomaly even under the GMSB framework, which is a very promising framework 
for its freedom from the SUSY CP- and flavor problems, but disfavored because it cannot simultaneously 
realize the 126 GeV Higgs and the muon g - 2 explanation. There, the "V" resurrects the GMSB. 
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Fig. 5.7 was the conclusive figure of this dissertation. There it was obviously shown that the 126 GeV 
Higgs boson can be realized with keeping the SUSY explanation of the muon g — 2 at the 2cr-level, while 
the lcr-level explanation was excluded by the LHC SUSY searches. Also it was emphasized that the mass 
parameter My, which governs the masses of the vector-like quarks and the vector-like lepton, is constrained 
as My < 1.2 TeV. This means the lightest vector-like quark t' x should be lighter than < 1.1 TeV, which is 
within the reach at the 14 TeV LHC. 

Then collider searches for the vector-like quarks were investigated. We saw that the search performed 
by the ATLAS collaboration [132] was very interesting and promising, and that improvements and better 
understanding of ^-tagging algorithms is awaited for further progress. 



On 17 December 2012, the LHC was shut down to prepare for collisions with £cm = 14 TeV, leaving 
the following message on the LHC monitor: 

*** End of operation for 2012! *** 
See you again briefly for p-Pb in 2013. 
High energy proton proton physics 
will be resumed in 2015. 
So long and thanks for all the fish. 

In 2015, we will obtain a much more powerful tool to investigate new physics beyond the Standard Model. 
There SUSY searches would strikingly proceed, and we expect we will obtain many clues for the next 
theories. Especially, the fate of the V-GMSB scenario is expected to be determined. 

Now we have to wait for three years, but now, contrary to that after the accident in 2009, we have data, 
which are enough to allow us to find the "tail" of the physics beyond the Standard Model buried inside them. 
We should exhaust the obtained data to be ready for the 14 TeV run starting in, hopefully, 2015. 
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